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WELCOMING  ADD  RES 8 


Roman  R.  Birukoff 
Air  Force  Special  Weapons  Center 
Directorate  of  Research 
Structures  Division,  SUES 


In  reviewing  the  chronological  events  leading  to  this  First  Shock  Tube 
Symposium,  I  must  mention  that  about  two  years  ago  when  Wesley  Curtiss  of  BRL 
and  I  were  discussing  progress  in  shock  tube  research,  I  proposed  then  the 
idea  of  a  Shock  Tube  Symposium.  He  agreed  that  the  idea  was  sound  and  that 
there  was  a  need  by  all  to  integrate  experiences  and  exchange  Information 
among  the  numerous  shock  tube  laboratories  that  had  germinated  after  the  war. 
Tears  went  by.  Wesley  Curtiss  left  BRL  and  took  up  ceramics.  I  remained  in 
Air  Force  shock  tube  laboratory  work,  presently  under  AFSWC. 

Last  December  with  the  approval  and  cooperation  of  Col.  2.  Oilier, 
Chief  of  Directorate  of  Research  of  the  Air  Force  Special  Weapons  Center,  and 
Mr.  Eric  H.  Wang,  Chief  of  Structures  Division  of  the  same  center,  question¬ 
naires  were  mailed  for  an  AFSWC -sponsored  first  Shock  Tube  Symposium.  The 
prompt  enthusiastic  response  permitted  us  to  plan  this  symposium  at  such  an 
early  date.  Unfortunately,  neither  Col.  E.  Oilier  nor  Eric  H.  Wang  are  able 
to  attend  -  Col.  Glller  because  of  other  duties,  Mr.  Wang  because  he  became 
ill  just  two  weeks  ago.  Therefore,  in  behalf  of  the  Air  Force  Special  Weapons 
Center,  I  welcome  you  here.  I  can  assure  you  that  the  Center  deeply  appre¬ 
ciates  your  enthusiastic  response  and  participation  in  this  venture. 

I  shall  now  turn  the  meeting  over  to  Mr.  Jack  Kelso  of  AFSWP  who  will 
be  the  Moderator. 


Til 


ARMED  FORCES  SPECIAL  WEAPONS  PR0JBCT*3  INTEREST 
IN  SHOCK  TUBE  RESEARCH 

Jack  R.  Kelso 

Armed  Forces  Special  Weapons  Project 


I  vould  like  to  discuss  the  use  of  the  shock  tube  from  the  point  of 
view  of  research  In  the  weapons  effects  field.  Other  speakers  will  discuss 
the  design,  development,  and  operation  of  shock  tubes,  and  their  associated 
equipment,  however,  I  would  like  to  consider  their  use  In  terms  of  the  infor¬ 
mation  thus  developed.  The  use  of  the  shock  tube  in  weapons  effects  research 
began  with  studies  of  basic  phenomena  by  Mr.  L.  G.  Smith  at  Princeton  during 
World  War  II.  Information  on  shock  reflections  and  interactions  was  obtained 
at  that  time  which  supported  certain  other  theoretical  studies  being  conducted 
by  the  National  Defense  Research  Council. 

Of  course  most  of  this  work  was  applied  to  the  use  of  high  explosive 
bombs  and  demolition  charges.  Scwever,  some  of  the  basic  theory  was  carried 
over  Into  the  field  of  weapon  effect?  for  atomic  bombs.  An  example  of  such 
incorporation  Is  given  In  the  unclassified,  "Effects  of  Atomic  Weapons”  pub¬ 
lished  by  the  Atomic  Energy  Commission  and  the  Department  of  Defense  In  1950. 
About  thi3  same  time  work  on  the  diffraction  loading  of  rectangular  blocks  was 
Initiated  at  Pilnceton  by  Dr.  Bleakney  under  an  existing  Office  of  Naval  Research 
contract.  Shortly  afterward,  similar  work  commenced  at  the  University  of  Mich¬ 
igan  under  the  direction  of  Dr.  Otto  LaPorte,  which  was  subsequently  followed 
by  programs  of  a  similar  nature  at  the  Ballistic  Research  Laboratories  and  the 
Armour  Research  Foundation.  This  information  was  then  applied  to  the  analysis 
of  structures  exposed  to  full  scale  tests  in  the  Facific  during  Operation 
GREENHOUSE.  Since  that  time  a  number  of  shock  tubes  h&Ve  been  used  extensively 
In  the  weapons  effects  field  to  develop  Information  on  the  blast  loading  of 
various  types  of  structures  and  other  objects  as  well  as  certain  aspects  of 
basic  phenomena.  With  this  background,  I  should  like  to  discuss  the  general 
philosophy  governing  the  use  of  shock  tubes  In  research  and  development  pro¬ 
jects  that  are  supported  by  the  Armed  Forces  Special  Weapons  Project.  This 
headquarters  is  a  tri-service  organization  which  conducts  research  in  the 
weapons  effects  field.  Usually  these  contracts  are  administered  through 
various  groups  within  the  three  Services. 

Generally  speaking,  three  types  of  investigations  are  carried  on  in 
the  shock  tube.  The  studies  usually  concern  the  following! 

(1)  Basic  phenomena 

(2)  Loading  of  objects 

(3)  Response  of  objects. 

Most  of  tbe  effort  in  the  weapons  effects  field  is  concerned  with  basic 
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phenomena  and  loading  of  objects,,  since  only  limited  response  tests  can  be  con¬ 
ducted  in  the  shock  tube  because  of  size  and  c^her  considerations.  Testing  of 
scaled  ircdels  representing  buildings  and  similar  civil  engineering  structures 
to  obtain  response  data  is  usually  unsatisfactory  since  modelling  cannot  always 
account  for  mass  and  gravity  effects.  Fcr  this  reason  response  tests  sure 
usually  confined  to  prccf  tests  of  small  components.  In  the  latter  case,  such 
objects  are  exposed  In  a  shock  tube  to  determine  their  response  characteristics 
Insofar  as  damage  and  translation  are  concerned. 

The  stocks  themselves  may  he  either  flat-tcpped  waves  or  peaked  waves, 
depending  oc.  their  method  of  generation  (either  with  compressed  air  or  high 
explosives),  as  well  as  toe  petit  of  observation.  Basic  phenomena  are  usually 
studied  by  interferometer  techniques  and  Schlleren  photography .  These  methods 
may  also  be  used  to  obtain  information  the  loading  of  models;  but  a  great 
deal  of  analysis  Is  required  to  accurately  define  the  loading  pattern  in  this 
way.  Consequently,  in  recent  years,  miniature  pressure- time  gages  have  been 
positioned  within  solid  models  to  octain  loading  Information  at  various  points 
on  the  surface.  In  this  way  data  can  be  quickly  obtained  on  the  same  model 
for  various  s.nock  strengths  and  durations-.  This  Information  is  then  reduced 
and  used  to  develop  generalized  leading  schemes  vhicl  are  then  employed  to 
predicts  the  blast  loading  cf  a  tail  sized  s4ructure  under  various  conditions. 
With  a  predicted  loading,  the  response  of  an  actual  or  assumed  structure  can 
then  be  determined  analytically. 

The  Armed  Forces  Special  Weapons  Project  shock  tube  research  program 
is  primarily  designed  tc  obtain  lnforroaMcr.  on  basic  phenomena  and  the  leading 
of  structures.  The  specific  requirements  for  Information  are  concerned  with 
problems  common  to  the  three  Services.  In  many  areas  detailed  information  can 
be  obtained  by  a  snock  tube  pre-gram  at  a  very  reasonable  cost  while  it  could 
not  in  the  field  because  of  economic  or  practical  reasons.  For  example,  en- 
tranceways  and  ventilation  systems  can  be  studied  in  the  shock  tube  under  a 
variety  of  conditions.  In  this  way  almost  an  infinite  variation  of  parameter® 
can  be  obtained. 

In  discussing  the  areas  cf  interest  to  the  Armed  Forces  Special  Wea¬ 
pons  Project  in  shock  tube  research  I  would  like  to  refer  back  to  the  three 
main  categories  of  such  work  previously  mentioned.  Those  investigations  that 
are  primarily  concerned  with  the  effect  of  the  environment  on  the  6hock,  such 
as  surface  roughness  or  the  reflection  process,  are  considered  to  be  studies 
of  basic  phenomena.  Those  investigations  which  Involve  the  diffraction  of 
shocks  around  objects  are  considered  to  be  loading  studies.  An  example  of  the 
latter  is  a  test  program  presently  underway  which  involves  variously  shaped 
structures  such  as  arches,  cylinders  and  domes.  These  are  designed  to  improve 
our  knowledge  of  weapons  effects  by  determining  the  relation  between  free 
field  blast  phenomena  and  the  loading  of  representative  targets  which  can  then 
be  extended  to  full  scale  conditions  through  the  proper  scaling.  In  many 
cases  the  requirement  for  full  scale  tests  can  be  reduced  or  eliminated  through 
shock  tube  research.  In  other  cases,  information  can  be  obtained  which  will 
assist  in  planning  or  analyzing  the  results  of  a  full  scale  test.  An  example 
of  response  testing  would  be  the  exposure  of  dynamic  pressure  gages  and  other 
full  scale  instrumentation  in  the  shock  tube. 

These  examples  of  shock  tube  research  that  I  have  Just  described  art; 
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some  of  those  which ''are  in  the  current  Aimed  Forces  Special  Weapons  Project 
program.  I  should  like  to  mention  briefly  those  areas  in  which  further  work 
seems  necessary.  These  concern  primarily  the  development  of  shock  tube  capa¬ 
bility  to  provide  higher  shock  strengths  and  longer  durations  than  those  cur¬ 
rently  in  use.  With  regard  to  instrumentation,  it  appears  desirable  to  find 
a  more  direct  method  of  measuring  particle  velocity  and  dynamic  pressure  in 
the  shock  tube,  in  view  of  the  departure  of  shock  tube  flows  from  the  ideal. 

For  this  reason  we  would  like  to  know  what  the  flow  conditions  actually  are 
so  that  the  drag  coefficients  developed  from  shock  tube  data  for  various  models 
may  be  verified  and  correlated  with  wind  tunnel  data.  Of  course  more  infor¬ 
mation  would  be  obtained  on  transient  drag  coefficients  when  longer  durations 
are  obtained  in  the  shock  tube. 

In  conclusion,  I  should  like  to  say  that  a  great  deal  of  useful  infor¬ 
mation  has  been  obtained  in  the  weapons  effects  field  by  means  of  shock  tube 
research.  It  appears  that  present  operational  procedures  are  adequate  to 
obtain  data  on  current  problems  and  that  this  is  being  done  in  a  very  efficient 
manner.  However,  future  work  seems  necessary  in  those  areas  of  Interest  I 
have  Just  mentioned  in  order  to  meet  future  requirements. 


THE  AIR  FORCE  SHOCK  TOBE  AS  A  RESEARCH  TOOL 

T.  H.  Schiffman 
Armour  Research  Foundation 


The  main  purpose  of  the  Air  Force  shock  tube  as  a  research  tool  la  to 
simulate  large  explosions  under  laboratory  conditions  In  order  to  study  basic 
flows,  loading,  and  response  of  model  structures.  Auxiliary  applications  in¬ 
clude  its  use  as  a  short  duration  wind  tunnel  for  the  study  of  muzzle  blast 
effects  on  flying  aircraft,  and  as  an  explosive  decompression  chamber  simula¬ 
ting  a  pressurized  aircraft  cabin  suddenly  vented  as  a  result  of  enemy  flr«. 

In  view- of  tte  complexity  of  such  problems,  the  solutions  do  not  lend 
themselves  readily  to  an  analytical,  theoretical  approach.  Experimental  in¬ 
vestigations  have  been  utilized  in  order  to  obtain  *he  required  information. 

Initially,  loading  information  on  relatively  simple  geometric  obstacles 
was  obtained  from  tests  in  small  shock  tubes.  These  were  limited  to  two- 
dimensional  configurations  because  of  the  available  instrumentation  techniques. 

Later,  loading  information  resulted  from  two-and  three-dimensional 
shock  tube  studies,  and  full  scale  field  evaluations.  Shock  tube  tests  under 
controlled  laboratory  conditions  augment  and  help  design  full  scale  studies 
which  must  be  limited  in  number  and  are  fairly  expensive. 

Two-dimensional  shock  tube  studies  were  extremely  valuable  In  deter¬ 
mining  t.he  fhenomenology  of  flow  and  loading,  but  the  models  usually  are  not 
representative  of  the  majority  of  actual  structures.  Three-dimensional  shock 
tube  studies,  can  be  of  immense  value.  They  are  sufficiently  flexible  so 
that  the  most  representative  of  actual  loading  conditions  -  in  terms  of  both 
blast  and  structural  parameters  -  can  be  obtained. 

Three-dimensional  shock  tube  studies  can  be  used  to  maximum  advantage 
where  models  are  larger,  more  detailed,  and  complex.  Hence,  the  shock  tube 
itself  must  be  relatively  large  so  that  wave  reflection  from  the  *-ube  walls 
will  not  interfere  with  the  flow  field  being  observed.  This  concept  led  to 
the  Air  Force  sponsorship  of  a  6  foot  diameter  shock  tube  at  Gary,  Indiana. 
Feasibility  studies  for  methods  of  shock  generation  led  first  to  the  erection 
of  an  Air  Force  2-foot  shock  tube  at  Ft.  Sheridan,  Illinois.  Today  both 
3hock  tubes  are  operated  with  the  same  equipment  at  Gary,  Indiana. 

Loading  measurements  on  simple  shapes  have  been  augmented  by  theore¬ 
tical  analysis  and  engineering  Judgment  whenever  possible  in  order  to  provide 
ready  extrapolations  from  simpler  to  more  complex  objects. 


In  the  last  decade  the  shock  tube  has  become  an  Important  tool  for 
gas  dynamics  research  with  application  to  blast  effects  research,  to  transonic 
aerodynamics,  aerodynamics,  high  temperature  gas  dynamics,  and  hypersonic 
aerodynamics.  The  literature  concerning  shock  tube  theory,  design,  and  ap¬ 
plication  is  so  extensive  it  would  require  an  extremely  lengthy  discussion 
to  cover  the  subject  adequately.  However,  since  most  of  this  material  Is 
familiar,  we  shall  Just  triefly  define  the  flow  regions  of  interest  of  the 
Air  Force  shock  tube  and  acme  of  its  operational  characteristics. 

Figure  l  defines  the  shock  tube  as  a  piece  of  pipe  of  uniform  cross 
section,  closed  at  both  ends,  with  a  thin  diaphragm  separating  the  tube  Into 
low  pressure  and  high  pressure  regions.  In  the  Air  Force  shock  tube,  shock 
wave  generation  is  accomplished  by  one  of  three  methods:  compressed  air, 
explosive  gases,  or  Primacord.'  Figure  1  illustrates  the  basic,  compressed 
air  shock  tube  theory.  The  expansion  and  combustion  chambers  are  identified 
respectively,  as  (l)  and  (4).  At  time  (t?  zero,  regions  (4)  and  (l)  are 
filled  with  gas  at  a  predetermined  pressure  ratio  (Pj^),  are  rest  and 
In  thermal  equilibrium. 

On  rupture  of  the  diaphragm,  a  wave  cf  compression  and  another  of 
rarefaction  are  formed,  moving  in  opposite  directions.  A  contact  surface  is 
formed  between  these  two  waves  which  ac^s,  theoretically,  like  a  piston  in  a 
cylinder.  Ahead  of  the  piston  is  t.:ie  s bcc it  wave,  while  behind  It  is  the 
rarefaction  wave.  The  shock  wave  compresses  the  gas  in  the  expansion' chamber 
to  an  intermediate  pressure  and  accelerates  its  flow  in  the  same  direction  as 
the  shock  wave.  The  rarefaction  wave  reduces  the  press'.*’-*  in  the  combustion 
chamber  to  the  ssrae  intermediate  pressure,  and  accelerates  the  fluid  parti¬ 
cles  in  the  direction  opposite  to  the  advance  of  the  rarefaction  wave,  so 
that  the  latter  follows  and  eventually  overtakes  the  contact  surface  and  the 
shock  wave.  The  contact  surface  divides  two  regions  of  constant  pressure 
and  particle  velocity,  but  unequal  temperatures  and  hence,  Mach  numbers, 

(2)  formed  between  the  shock  wave  and  the  contact  surface  (hot  flow  region). 
And  (3)  formed  between  contact  surface  and  incident  rarefaction  wave  (cold 
flow  region).  Some  time  later,  conditions  exist  in  which  the  ref lec ted .rare¬ 
faction  wave  overtakes  the  incident  rarefaction  wave.  At  a  slightly  later 
time,  the  former  wave  overtakes  the  contact  surface.  Still  later,  i V  over¬ 
takes  the  shock  wave,  leading  to  the  customary  peaked  shock  profile,  and  the 
shock  front  is  followed  by  &  decaying  wave,  as  generated  in  atomic  blasts.. 

For  example,  compressed  air  up  to  an  initial  pressure  of  350  psi  1® 
presently  used  in  the  2-foct  tube,  giving  30  pBi  overpressure  at  the  shock 
front.  Compressed  air  also  seems  to  offer  the  best  solution  for  shocks  of 
long  durations.  In  the  6-foot  tube,  the  method  demands  the  provision  of  con¬ 
siderable  compressor  capacity.  *  Such  a  compressor  has  been  installed  and  la 
In  operation.  Vacua  to  1/2  psi  absolute  are  obtained  in  the  6-foot  tube  and 
to  5  mm  Hg  in  the  2-foot  tube.  The  timewise  variation  of  flow  In  the  tube 
is  well  known,  and  the  flow  variables  can  he  computed  for  a  given  pressure 
ratio  across  the  diaphragm. 

In  Figure  2  we  see  the  shock  pressure  ratio  plotted  as  a  function  of 
diaphragm  pressure  ratio  for  various  combinations  of  g&3es  on  the  high  and 
low  pressure  sides  of  the  diaphragm.  Such  a  chart  is  ueed  for  calibrating 


the  pressure  sensors  dynamically  in  a  4-inch  shock  tube,  .typical  pressure¬ 
time  variations  are  illustrated  in  Figure  3  wherein  the  pressure  distribution 
in  the  shock  tube  is  shown  at  three  different  times,  t  =  0,  1.7  milliseconds 
later,  and  finally  5*2  milliseconds  after  diaphragm  rupture.  Let  us  now  con¬ 
sider  the  shock  tube  flow  regions  and  indicate  some  of  its  applications. 

For  weapons  effects  work,  we  utilize  the  compressed  air  shock  wave  in 
the  expansion  section  at  a  distance  of  16  compression  chamber  lengths  which 
results  in  a  peaked  profile.  If  Prlmacord  is  used,  the  profile  is  peaked 
earlier.  The  shock  wave  is  made  to  diffract  around  models  to  be  loaded.  These 
might  be  simple  shapes  (blocks,  cylinders,  domes,  etc.)  or  complex  shapes  such 
as  multistory  buildings.  "Incidence"  refers  to  the  Mach  stem  or  regular  re¬ 
flection  region  created  by  the  tilting  plate  on  which  the  model  rests.  In 
structures  tests,  the  end  of  the  shock  tube  can  be  left  open  and  models  can  be 
observed  outside  of  the  shock  tube,  etc. 

For  muzzle  blast  tests,  we  utilize  the  flat  top  flow  region  as  a  short 
duration  wind  tum  el.  (The  duration  is  long  however  compared  to  the  duration 
of  the  muzzle  bla.it  wave.)  The  hot  flow  region  furnishes  Mach  numbers  up  to 
I.89,  and  the  cold  flow  region  up  to  much  higher  Mach  numbers.  The  shock  tube 
can  be  evacuated  to  simulate  any  desired  altitute.  Thus  by  firing  a  50~caliber 
gun  upstream  into  the  compression  chamber,  one  simulates  on  an  instrumeuted 
plate  the  muzzle  blast  effect  of  en  aircraft  firing  while  flying  both  subsonic 
and  supersonic.  It  was  found  by  emoirical  means  that  muzzle  blast  could  be 
correlated  to  the  detonation  of  a  spherical  charge  of  TNT.  In  this  way  the 
various  scaling  laws  of  peak  pressures,  impulse,  ambient  temperatures,  and 
ambient  pressure  changes  for  spherical  explosives  could  be  applied  to  muzzle  _ 
blast. 


For  explosive  decompression,  one  utilizes  the  compression  chamber. 
Tests  in  the  6-foot  shock  tube  were  conducted  on  equipment  furnished  by 
Sperry  Gyroscope  Company  to  evaluate  the  structural  response  of  the  Sperry 
Unit  and  its  mounting  when  subjected  to  explosive  decompression.  Such  explo¬ 
sive  decompression  world  be  experienced  during  high  altitude  flight  if  a 
pressurized  cabin  is  suddenly  vented  as  a  result  of  enemy  fire. 

Another  use  is  for  the  check  out  of  equipment,  both  for  shock  loading 
or  explosive  decompression  loading. 

One  of  the  earliest  uses  of  the  shock  tube  was  to  study  the  propaga¬ 
tion  and  interaction  of  shock  waves.  At  Princeton  University  considerable 
effort  was  expended  on  fundamental  studies  of  shock  propagation  and  reflec¬ 
tion,  as  well  as  blast  loading  on  structures.  Figure  4  shows  the  Princeton 
double  block.  For  this  application,  the  shock  wave  in  the  tube  simulates  the 
wave  resulting  from  an  explosion,  and  one  can  determine  how  the  pressure  on 
a  model  varies  with  time  as  the  shock  wave  passes  over  and  downstream  of  it. 
For  these  studies  a  tube  of  rectangular  cross  section  is  desirable,  and  use 
is  made  of  a  Mach-Zehnder  interferometer  in  order  to  obtain  a  quantitative 
measurement  of  the  density  field.  (The  shock  tube  used  in  the  Princeton  ex¬ 
periments  was  38  feet  long  and  had  a  4  inch  by  18  inch  cross  section.) 

An  instantaneous  picture  of  the  flow  field  in  the  test  section  can  be 
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obtained  with  a  short  duration  light  source.  In  order  to  get  a  time  history 
of  the  flow,  a  variable  time  delay  is  incorporated  in  the  spark  source  so 
that  on  repeated  runs,  a  sequence  of  pictures  taken  at  the  various  time  in¬ 
tervals  will  yield  the  desired  result.  This  is  made  possible  by  the  excellent 
reproducibility  of  shock  tube  runs.  The  interferometer  measures  the  density, 
and  the  pressure  can  be  computed  under  the  assumption  of  isentropic  flow  in 
regions  devoid  of  shock  waves.  Hence,  it  is  possible  to  obtain  the  pressures 
on  an  object  as  a  function  of  time  for  the  case  of  a  shock  wave  striking, 
passing  over,  and  going  downstream  of  the  object.  Inherent  in  this  method  is 
the  limitation  of  two-dimensionality.  It  is  possible,  under  suitable  condi¬ 
tions,  to  make  quantitative  three-dimensional  interferograms;  however,  for 
arbitrary  shapes  this  is  not  possible.  In  order  to  establish  three-dimensional 
effects  a  tube  of  large  cross  section  is  essential. 

The  Foundation  also  owns  an  8-inch  by  8-inch  square  sectional  shock 
tube,  30  feet  in  length,  located  at  Technology  Center.  This  shock  tube  is 
equipped  with  interferometer  and  Schlieren  systems  which  permit  the  explora¬ 
tion  of  scaled  pressure-time  boundary  conditions  for  the  Air  Force  shock  tubes. 
Whereas  in  interferometric  work  continuous  pressure-distance  plots  are  ob¬ 
served  at  discrete  times,  with  the  Air  Force  shock  tube  we  observe  continuous 
pressure-time  plots  at  discrete  points. 

Under  the  sponsorship  of  the  U.  S.  Air  Force,  the  Armour  Research 
Foundation  has  placed  into  operation  shock  tubes  which  will  permit  studies 
of  flows  past  three-dimensional  models.  Figure  5  is  an  artist's  sketch  of 
the  Air  Force  Shock  Tube  Laboratory  located  in  Gary,  Indiana.  Figure  6  shows 
the  outside  and  Figure  7  is  an  inside  view.  This  shock  tube  laboratory  con¬ 
tains  the  following: 

A  4-inch  diameter  tube,  15  feet  long  which  is  used  for  calibrating 
gages . 

A  2-foot  diameter  tube,  200  feet  long  composed  of  10  lengths,  each 
20  feet  long,  with  a  combustion  chamber  than  can  be  varied  in 
length  from  4  inches  to  20  feet  by  means  of  wood  blocks. 

A  6-foot  diameter  tube,  lr0  feet  long.  Thi^  length  comprises  a 

minimum  combustion  chamber  3  feet  long,  a  5“f°°t  test  section,  and 
ordinarily  a  40-foot  exhaust  section,  leaving  100  feet  maximum  or 
expansion  section.  The  tube  thickness  is  5/8  of  an  inch  in  the 
expansion  section,  and  1  inch  in  the  compression  chamber.  As  a 
result  of  using  interchangeable  tube  sections,  with  the  majority 
10  feet  long,  it  is  possible  to  s»t  up  a  vide  range  of  shock  wave 
pressures  (over  ambient  pressures)  and  durations,  the  former  from 
1  psi  to  50  psi  and  the  latter  variable  between  1  and  50  milli¬ 
seconds. 

Normal  test  duration  is  15  to  20  milliseconds.  Loading  of  the  combustion 
chamber  and  access  to  the  test  section  is  effected  by  mounting  these  tube 
sections  on  dollies  which  run  on  rails  extending  the  full  length  of  the  tube. 
The  tube  is  primarily  located  in  the  open,  with  a  building  to  protect  the 
test  section  and  adjoining  lengths,  as  veil  as  the  instrumentation.  Having 
both  ends  free  enables  the  tube  to  be  fired  with  the  expansion  chamber  open 
or  closed.  The  test  section  and  the  two  IG-foct  sections,  respectively, 
leading  and  trailing  it,  are  provided  with  10-foot  horizontal  stainless 
steel  floor  plates  to  accommodate  the  test  model  and  to  smooth  out  irregu- 


larities  In  the  shock  wave  before  and  after  striking  the  model.  The  model 
can  be  mounted  on  a  20-inch  turntable  and  rotated  about  a  vertical  axis  to 
present  any  face  to  the  Bhock  front.  Figure  8  shows  the  turntable  with  a 
thr  e-dimensional  double  block. 

The  work  is  facilitated  by  such  improvements  as  overhead  cranes  and  a 
puncturing  device  and  clamps  as  shown  in  Figure  9*  The  clamp  has  ten  hydrau¬ 
lic  Jacks,  each  with  a  80,000  psi  capacity.  Shock  generation  is  compressed 
air,  an  explosive  gas  mixture,  or  Prlmacord.  Compressed  air  is  used  with 
Mylar  up  to  15  pal  and  aluminum  up  to  65  psi  compression  chamber  pressure. 

Explosive  gas  is  used  in  the  form  of  a  hydrogen-air  mixture  of  1  to 
4  at  atmospheric  pressure.  A  6hort  length  of  Prlmacord  has  been  found  to  be 
the  most  effective  means  of  detonating  the  mixture  with  a  minimum  of  burning. 
As  such,  it  satisfactorily  fills  the  range  of  peak  shock  duration  up  to  70 
milliseconds  with  pressures  up  to  10  psi.  The  material  most  commonly  used  in 
the  compression  chamber  is  Prlmacord,  colled  in  a  flat  spiral  at  right  angles 
to  the  axis  of  the  tube.  This  method  offers  more  versatility  over  the  range 
of  practical  weapons  effects  pressures  and  durations. 

The  6-foot  Air  Force  shock  tube  was  placed  in  operation  on  September 
1954,  and  since  has  functioned  continuously  with  up  to  20  shots  a  day  as 
routine. 


Figure  10  shows  the  instrument  room.  The  instrumentation  records 
pressures  in  the  shock  wave  and  cn  the  model,  and  the  velocity  of  the  shock 
front  over  a  period  of  time  represented  by  steadyflow  conditions.  Therefore, 
all  electronic  and  photographic  equipment  available  at  the  shock  tube  labora¬ 
tories  is  controlled  by  a  time  delay  generator  furnishing  dexays  from  1  micro¬ 
second  to  1  second,  initiated  from  the  combustion  chamber.  This  enables  any 
selected  portion  of  the  shock  phenomena  to  be  recorded  on  the  appropriate 
instrument.  For  calibrating  the  gages  dynamically  in  the  4-inch  shock  tube, 
the  shock  wave  velocities  are  measured  electronically  by  means  of  light 
screens  viewing  through  windows  set  at  predetermined  intervals  in  the  tube 
and  upstream  of  the  model  position.  From  the  velocity  one  calculates  the 
overpressures  by  means  of  one-dimensional  shock  tube  theory. 

An  8-megacycle  recorder  registers  time  Intervals.  The  light  screen 
photocell  output  is  fed  to  an  amplifier  and  then  to  the  chronograph  where 
immediate  visual  indication  is  given  by  a  bank  of  lights. 

The  measurement  of  response  pressures  on  the  model  itself,  and  of 
the  shock  pressure  generated,  is  based  on  the  use  of  special  barium  titanate 
gages.  The  material  possesses  piezoelectric  sensitivity  and  produces  sun  out¬ 
put  signal  proportional  to  the  stress  applied.  This  type  of  gage  satisfies 
all  the  requirements  for  gages  usable  on  small  models,  being  diminutive  in 
terms  of  the  model  dimensions  and  effecting  no  change  in  model  geometry  when 
installed.  In  addition,  the  gages  respond  in  less  than  20  microseconds,  havs 
a  sensitivity  suited  to  conventional  amplifiers,  and  are  relatively  insensi¬ 
tive  to  acceleration.  The  gage  is  a  ceramic  cylinder  l/8-inch  long  by  l/8-inch 
outside  diameter,  with  a  wall  thickness  of  19  mils. 

8evcral  such  gages  are  Installed  in  the  faces  of  the  model  and  are 
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connected  tc  cathode  ray  oscilloscopes.  The  6-foot  tube  Is  equipped  for  eight 
channels,  bat  this  number  of  channels  can  be  Increased  as  required. 

Eash  oscilloscope  may  be  triggered  to  give  single  sweep  records  at 
faster  sweep. speeds  than  the  drum  camera  could  record.  An  S.  I.  M.  No. 731 
-amera,  with  a  writing  speed  of  1  millisecond  per  inch,  makes  eight  simultan¬ 
eous  records,  the  oscilloscopes  being  staggered  in  the  vertical  phase  and 
arranged  on  the  circumference  of  a  circle.  An  L.  F.  E.  oscilloscope  with  sig¬ 
nal  Input  from  a  gage  records  the  overpressure  for  each  test. 

The  gages  which  have  been  used  are  the  barium  titanate  gage,  the 
Klssler  gage,  Rutico  gage,  and  SR4  strain  gages.  In  the  future  other  more 
powerful  piezo  electric  materials,  likjj  zirconate  will  be  Investigated. 

Before  loading  models  one  must  be  reasonably  sure  to  obtain  the  desired 
Incidence  freestream  overs re ssures  and  durations. 

Figure  11  gives  the  shock  strength  as  a  function  of  scaled  distance  for 
any  atmospheric  pressure  using  Primacord  as  the  shock  generator.  The  solid  line 
has  also  been  obtained  from  an  analytic  solution  independent  of  the  test  cali¬ 
bration.  By  using  reduced  atmospheres  one  obtains  high  shock  strengths  for 
relatively  low  overpressures,  low  forces,  and  low  acceleration  signals,  etc. 

For  example,  one  can  obtain  shocks  of  {  =  7  corresponding  to  over  100  pel  in 
the  test  section  with  adequate  durations  for  overpressures  of  15  psi  corres¬ 
ponding  to  shocks  of  £  =  2  at  atmospheric  pressure. 

Figure  12  gives  the  scaled  duration  as  a  function  of  the  scaled  dis¬ 
tance. 

Figure  13  to  15  shows  typical  weapons  effects  models.  Figure  13  ex¬ 
hibits  a  scaled  beam;  Figure  14  a  multistory  building.  Figure  15,  the  same 
building  mounted  on  the  tilting  plate  for  study  in  the  regular  reflection 
region.  Figure  l6  and  17  exhibit  Sperry  equipment  installed  in  the  compres¬ 
sion  chamber  for  explos^  ,  decompression  testing.  Figure  16  shows  the  front 
of  the  unit;  Figure  17  lo.  king  towards. the  panels  exhibits  this  unit  inside 
the  tube  and  incidently  shows  the  bars  of  the  puncturing  device  which  hold 
the  diaphragm  in  place  insuring  a  more  uniform  rupture. 

Some  of  the  weapons  effects  loading  programs  include  a  study  of  long 
versus  short  durations,  shielding,  surface  roughness,  regular  reflection, 
hollow  models  with  and  without  interior  obstacles,  multistory  buildings, 
drag-type  buildings,  tanks,  domes,  a  study  of  shock  tube  temperatures,  etc. 

A  few  limited  response  studies  were  successfully  conducted  in  which 
the  effect  of  gravitational  forces  could  be  neglected,  at  in  case  of  wall 
panel  failure. 

As  an  example  of  equipment  checkout,  one  might  cite  a  few  instrumen¬ 
tation  checkout  of  side-on  gageB  like  the  Rutico  gage,  "q"  gages,  and  more 
complex  telemetering  and  oscillographic  equipment  to  be  used  on  large  scale 
field  tests. 
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A  test  to  be  conducted  soon  will  be  the  study  of  grids  to  eliminate 
the  reflections  from  the  back  end  of  the  tube  or  from  the  tilting  plate  thus 
permitting  longer  observations. 

In  conclusion,  the  Air  Force  shock  tube,  which  is  the  largest  shock 
tube  in  this  country  and  probably  in  the  world,  provides  an  excellent  research 
tool  for  (a)  weapon  effects  work  which  includes  basic  flow,  loading  on  simple 
and  complex  targets,  limited  response  tests,  and  checkout  and  calibration  of 
equipment  for  full  scale  tests,  (b)  a  6hort  duration  wind  tunnel  to  conduct 
muzzle  blast  tests  simulating  firing  from  planes  in  flight,  and  (c)  explosive 
decompression  studies  which  include  structural  tests  and  equipment  checkout  in 
the  compression  chambers.  These  are  some  of  the  present  applications.  Un¬ 
doubtedly  in  the  future  many  other  applications  will  be  found. 

Recently  the  6-foot  shock  tube  facility  has  been  modified  by  strength¬ 
ening  the  compression  section  and  nearly  doubling  the  tube  length.  This  will 
permit  overpressures  up  to  60  psi  in  the  test  section  ul  th  associated  durations 
up  to  120  milliseconds. 


SHOCK  MACH  NUM 


Various  Ratios  of  Sound  Velocity  in  High  and  Low  Pressure  Cases 
(Diatomic  Driver  Gas;  Monotomic  Driven  Gas) 


-  12  - 


» 


Pressure 


760  mm 
23*C 

Time  *  0  ms 

1 5  mm  23aC 

Distance  ► 


Figure  3.  Typical  Shock  Tube  Flow  Pressure-Time  Variation 
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Figure  4.  Princeton  Double  Block 
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SHOCK  TUBE  LABORATORY 
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Figure  6.  View  of  2  ft.  Diameter  (foreground)  and  6  ft  Diameter  Shock  Tube 
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Figure  7.  6  Foot  Diameter  Shock  Tube 
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Figure  7.  6  Foot  Diameter  Shock  Tube 
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Figure  10.  Instrument  Room 


figure  12.  Scaled  Duration  vs  Scaled  Distance 
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Figure  13.  Scaled  Beam 
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Figure  15.  Multistory  Building  in  Regular  Reflection  Region 


,;,.%v\vy 
*»\*  i.\*  <  AY 


P®iW?p"*eiSli 


V 


SWisppS: 


m  ■  :  n. 

i  .  f 

,\'»\>V\V"  •:’  •-■•■•.'.■  >••  WVVWrt'WV  '  .  .".W'*»V*\*iW*'  v».  >  ■.' V- 


m 

pm 


}.V-yVvV 


*4  Ifilii 

Kpil 

'■  . v\  I »\\ •  .  ‘A  t«\-  .UHMMVir  u«\ni^i 


w%mb. 

«  «w%u  *» 


i4i^f'vIw 

8*  .v%AAwS\wS.  AVtVVXK 


Z? 


Figure  17.  Sperry  Navigation  Unit  behind  Puncturing  Device 


'j’HS  DESIGN  AND  PK  RFORMANCE 

OF  TK5  GENERAL  ELECTRIC  SIX -INCH  SHOCK  TUNNEL  FACILITY 

Walter  R.  Warren 
General  Electric 


The  design  of  the  General  Electric  Missile  and  Ordnance  Systems  De¬ 
partment  six-inch  hypersonic  shock  tunnel  is  discussed  and  the  advantages  of 
several  unique  design  features  -  various  test  regions,  multiple  nozzles,  dia¬ 
phragm  clamping  section,  nozzle  opening  section,  sting  -  are  indicated.  The 
operation  of  the  tunnel  is  described  in  terms  of  its  maxiirum  performance.  In¬ 
strumentation  techniques  used  for  measuring  model  heat  transfer  rates  and 
pressure  distributions  and  fo:  observing  flow  fields  are  described.  The  per¬ 
formance  of  the  facility  is  presented  in  terms  of  model  simulation  of  free 
flight  conditions  and  the  problem  of  non-equilibrium  flow  in  the  different 
nczzies  is  analyzed.  Important  experimental  programs  that  can  be  run  in  the 
shock  tunnel  are  outlined. 

The  aerodynamic  phenomena  associated  with  high  velocity  flight  at 
various  altitudes  are  highly  complicated  because  of  the  many  processes  that 
occur.  For  example,  the  air  coving  near  the  surface  of  a  body  travelling 
through  the  atmosphere  at  a  hypersonic  Mach  number  (over  10)  will  be  highly 
dissociated  and  partially  ionized,  it  will  radiate  appreciable  amounts  of 
energy  to  cooler  surroundings,  and  it  may  be  in  a  chemical  and  thermodynamic 
state  appreciable  far  from  the  equilibrium  state  associated  with  its  transla¬ 
tional  temperature.  This  flow  situation  is  difficult  to  analyze  and  the  need 
for  experimental  facilities  that  es-ablish  a  suit-able  environment  for  model 
testing  is  apparent.  One  facility  that  is  well  suited  for  adaptation  to  the 
study  of  high  velocity-high  temperature  problems  is  the  shock  tube  or  its 
modification,  the  shock  tunnel. 

The  General  Electric  6“inch  shock  tunnel  haB  been  built  specifically 
for  the  study  of  realistic  hyper  sonic  flight  problems.  It  is  essentially  a 
high  performance  shock  tube  in  which  relatively  large  size  models  may  be 
tested  over  wide  ranges  of  Mach  number,  Reynolds  number,  and  stagnation 
enthalpy.  The  facility  is  subjected  to  the  usual  shock  tube  limitations  and 
difficulties  (short  testing  times,  Incident  shock  wave  attenuation,  non¬ 
equilibrium  flows,  etc.);  however,  it  is  believed  that  its  design  and  opera¬ 
tion  characteristics  tend  to  minimize  these  problems. 

The  purpose  of  this  paper  is  to  describe  the  shock  tunnel  and  Its 
performance,  to  discuss  t.he  programs  that  have  been  conducted  during  its 
development,  and  to  outline  Its  applicability  t,o  several  important  hypersonic 
studies. 


Figure  l  shows  the  original  conception  of  the  shook  tunnel.  Many  of 
the  details  have  been  changed  during  the  development  of  the  facility  out  the 
overall  size  and  general  layout  of  equipment  are  approximately  the  sane. 

Three  flow  regions  of  high  stagnation  enthalpy  are  used  fcr  model  ex¬ 
perimentation  (see  Figure  2):  the  flow  behind  the  incident  shock  wave  - 
region  2  -  where  the  pressure  and  Reynolds  number  are  high  but  the  flew  Mach 
number  is  lew  (below  3);  the  flow  in  a  non-ref lected  or  straight-through 
nozzle,  where  the  pressure  on  a  model  is  lower  than  it  is  in  the  region  2  flow 
but  the  flow  Mach  number  is  higher  (between  U  and  6);  and  the  flew  in  a 
ref  lected  nozzle  where  the  pressures  are  quiet  low  but  the  Mach  nutters  are 
high  (to  16} .  The  size  of  the  models  located  in  the  nozzle  flews  can- be  five 
to  six  times  that  of  the  models  in  the  region  2  flow. 

An  interesting  facility  modification  for  reflected  nozzle  meeting  is 
the  addition  of  several  nozzles  to  the  downstream  end  of  the  driver,  tube. 

This  is  shown  in  Figure  1  and  in  more  detail  in  Figure  3*  Such  &  scheme  would 
increase  the  productivity  of  an  experimental  program  and  would  reduce  the 
repeatability  problem  existing  in  high  performance  shock  tubes,  sir  e  it  allows 
the  taking  of  many  lata  during  one  test.  The  installation  now  has  cniy  one 
nozzle  but  it  may  be  easily  adapted  to  several. 

The  dimensions  of  the  shock  tunnel  and  its  supporting  shock  tube  are 
gjiven  in  Figure  h.  The  short  shock  tube  has  the  same  diametric  dimensions  as 
tiie  tunnel.  It  vas  first  used  in  the  high  pressure  combustion  study  that  is 
described  in  a  later  section.  The  primary  purpose  cf  the  tube  is  tc  act  as  a 
development  facility  for  the  tunnel,  although  it  will  also  be  used  fcr  a  few- 
basic  studies.  Frcgrams  being  run  with  it  now  are  intended  to  develop  and 
evaluate  'neat  transfer  and  pressure  instrumentation  and  to  determine  the 
durability  and  response  cf  various  gauge-model  combinations. 

The  driver  tube  lias  an  8- inch  inner  diameter  and  the  walls  are  6 
Inches  thick.  It  is  made  of  SAF  steel  and  is  designed  to  contain  a  work- 

ing  pressure  10,000  psi.  It  has  been  hydrostatically  tested  to  cl,~Z-2  psi;  tne 
ac:tnal  rupture  strength  should  be  well  over  100,000  psi.  The  ime.  surface  is 
chrome-plated  to  a  depth  of  0.GC5  inches.  The  access  ports  are  '.crated  spi¬ 
rally  every  H5  degrees  at  an  axial  spacing  of  6  inches.  With  this  system  of 
sp-ark  plugs,  gas  inlets,  vacuum  lines,  and  gauges  may  be  placed  In  almost  any 
desired  array.  Two  diametrically  opposed  access  ports  are  located  every  10 
fc?et  in  the  upstream  sections  of  the  driven  tube  and  every  foot  fcr  the  final 
2(5  feet  of  the  tube.  The  ports  allow  a  1-inch  diameter  entrance  tc  the  tube 
aijid  are  equipped  with  blanks  that  were  bored  and  honed  in  place.  The  inner 
surface  of  the  tubes  have  a  12  to  16  micro-inch  honed  finish.  Region  2  models 
ai;e  supported  in  a  special  6-inch  section  that  fits  between  any  two  large  sec¬ 
tions  of  the  driven  tube. 

The  reflected  nozzle  is  shown  in  Figure  k.  The  throat  section  is  re¬ 
movable  and  rv  j  be  replaced  by  sections  with  various  sized  throats.  The 
nozzle  is  a  «  iegree  cone  which  causes  some  flow  divergence  at  the  test 
section.  However,  it  was  thought  better  to  accept  a  slight  flow  divergence 
titan  to  attempt  the  design  of  a  parallel  flow  nozzle  for  the  complicated  flow 
conditions  of  interest.  A  maximum  model  diameter  of  1  foot  i3  believed 
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possible  in  the  30-inch  test  section.  The  model  is  supported  on  a  sting 
extending  through  the  dump  tank.  Two  10-inch  diameter  Schlieren  systems 
(horizontal  and  vertical)  are  available  for  flow  observation  in  the  test 
section. 


Figure  5  is  a  photograph  of  the  installation  showing  the  driven  sec¬ 
tions  of  the  tube  and  of  the  tunnel.  The  control  room  is  seen  in  the  center 
of  the  picture  beside  the  nozzle  end  of  the  tunnel.  A  view  of  the  tunnel  . 
looking  from  the  driver  house  towards  the  dump  tank  is  shown  in  Figure  6. 

The  piping  to  the  left  of  the  driver  is  the  gas  inlet  system.  The  location  of 
the  driver  tube  in  a  specially  designed  room  external  to  the  main  building  Is 
part  of  the  overall  safety  system  which  also  provides  for  the  complete  remote 
control  operation  of  an  experiment  from  an  Interlocked  panel  in  the  control 
room. 


The  high  driver  pressure  needed  to  operate  the  tube  is  obtained  by 
one  of  three  methods;  helium  pressurized  to  2000  psi,  constant  volume  com¬ 
bustion  of  helium  and  stoichiometric  hydrogen  and  oxygen  mixtures  to  10,000 
psi,  and  constant  volume  combustion  of  nitrogen  and  stoichrometric  hydrogen 
and  oxygen  mixtures  to  300C  psi.  These  three  driver  techniques  were  chosen 
because  of  their  ability  to  cover  the  complete  range  of  desirable  test  condi¬ 
tions  and  because  of  their  practicality,  reliability,  and  safety  compared  to 
other  techniques.  Maximum  performance  is  obtained  with  a  rich  helium,  hydro¬ 
gen,  and  oxygen  mixture  although  very  little  advantage  is  gained  by  going  to 
mixtures  of  less  than  7*#  helium  (mole  fraction). 

Various  design  performance  characteristics  are  illustrated  In  the 
following  discussions.  A  10,OCO-psi  driver  produced  by  an  initial  1100-psi, 
helium  combustion  mixture  will  drive  a  Mach  20  shock  into  air  originally 
at  a  pressure  of  1  mm  of  Hg.  Considering  the  reflected  nozzle  flow  with  a 
3 A -inch  throat,  a  shock  Mach  number  of  6,  and  air  at  atmospheric  pressure 
initially  in  the  driven  tube,  equilibrium  temperature  and  pressure  in  region 
5,  which  is  the  stagnation  chamber  for  the  nozzle,  are  5000  psi  and  3^00*K, 
respectively.  The  equilibrium  flew  Mach  number  at  the  model  is  10.6.  This 
condition  is  produced  with  an  initial  combustion  charge  in  the  driver  of 
approximately  450  psi.  At  higher  shock  Mach  numbers  the  stagnation  tempera¬ 
ture  goes  up,  of  course,  but  the  flow  Mach  number  at  the  model  drops  because 
the  expanded  gas  is  still  hot. 

Several  aspects  of  the  design  of  the  tunnel  are  interesting.  One  is 
the  large  inner  diameter  of  the  driven  tube.  This  allows  a  relatively  long 
tube,  and  therefore,  long  testing  times  with  a  reasonable  length  to  diameter 
ratio.  The  shock  wave  attenuation  losses  are  thus  kept  low.  Another  impor¬ 
tant  feature  is  the  separation  of  the  different  tube  sections  through  sleeve- 
shaft  seal  designs.  This  has  been  done  to  reduce  the  transmission  of  vibra¬ 
tions  through  the  structure.  It  also  allows  simple  operation  of  the  tunnel. 
Figure  8  shows  the  three  sections  of  interest.  Consider  the  diaphragm 
clamping  system.  The  driver  is  attached  to  a  sleeve  that  slides  over  the 
driven  tube  through  four  hydraulic  cylinders.  The  thrust  of  each  cylinder 
at  full  pressure,  7000  psi,  is  250,000  pounds.  In  this  arrangement  only  the 
sleeve  and  the  hydraulic  pistons  move  when  a  diaphragm  is  changed  after  a 
run  and  this  is  done  simply  by  throwing  a  switch.  The  seals  are  achieved  by 
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shaft  mounted  "0"  rings  between  the  sleeve  and  the  driven  tube  and  by  face 
mounted  "0"  rings  between  the  diaphragm  and  both  the  sleeve  and  the  driver 
tube.  When  the  diaphragm  breaks,  it  hits  only  the  sleeve  and  the Vansmlssion 
of  vibrations  to  the  driven  tube  is  not  large.  The  sleeve  has  a  removable 
insert  in  the  flow  passage  before  the  area  contraction  which  permits  the  dia¬ 
phragms  to  break  out  into  a  round  or  square  support  section.  A  photograph  of 
this  system  with  only  two  hydraulic  cylinders  is  shown  in  Figure  8. 

The  nozzle  end  of  the  driven  tube  is  also  a  sleeve  unit  activated  by 
hydraulic  cylinders.  The  second  diaphragm  is  located  between  the  end  of  the 
driven  tube  and  the  sliding  block  which  slips  over  the  removable  throat  sec¬ 
tion  on  another  "0"  ring  shaft  seal.  Figure  9  is  a.  photograph  of  the  test 
section  end  of  the  shock  tunnel. 

The  bellows  mounted  seai  between  the  dump  tank  and  the  model  support¬ 
ing  sting  also  helps  to  isolate  the  model  from  structural  vibrations.  A 
sketch  of  the  complete  sting  and  drive  mechanism  is  given  in  Figure  10. 

The  severe  shock  tube  instrumentation  problems  are  well  known,  as  are 
several  methods  that  are  now  being  employed  to  make  desirable  measurements. 

For  model  instrumentation  in  the  M.O.S.D.  program,  three  basic  types  of  in¬ 
strumentation  are  being  employed:  thin  resistance  thermometers  to  measure 
surface  temperatures  from  which  heating  rates  may  be  deduced,  piezoelectric 
pressure  transducers  to  measure  body  pressure  distributions,  and  high  speed 
photcgrapty  ;o  study  flow  establishment  and  the  character  of  the  developed  flow. 
More  or  less  standard  instrumentation  techniques  are  used  to  operate  the  fa¬ 
cility  and  to  determine  its  performance.  A  few  of  these  are  promoter-type 
spark  plugs  with  variable  ignition  energy  systems  to  initiate  combustion, 
pressure  gauges  and  ionization  probes  used  in  conjunction  with  cathod  ray 
oscilloscopes  and  a  multichannel  recording  oscillograph  to  determine  the 
shock  velocity  and  the  tunnel  pressures  as  a  function  of  time,- and  thermo¬ 
couple  vacuum  gauges  and  absolute  manometers  to  measure  the  initial  pressures 
in  the  driven  tube  and  in  the  test  section. 

The  thin  resistance-type  thermal  gauges  have  been  developed  to  a  high 
productivity  at  Princeton  University,  the  AVCO  Advanced  Development  Research 
Laboratory,  the  Cornell  Aeronautical  Laboratory,  and  other  laboratories.  It 
is  planned  to  use  the  very  thin  (the  order  of  0.1  micron),  or  surface  temper¬ 
ature  following  gauge  as  often  as  possible  because  of  the  relative  ease  of 
depositing  such  a  gauge  on  a  doubly  curved  surface  without  causing  flow  dis¬ 
turbance.  For  high  Reynolds  number  testing,  when  the  durability  of  the  thin 
film  becomes  critical,  the  calorimeter,  or  thick  film  gauge  imbedded  in  the 
model  surface  will  be  used. 

A  number  of  experiments  have  been  made  with  a  miniature  Kistler  quartz 
crystal  pressure  gauge  located  at  the  stagnation  point  of  a  2-inch  hemisphere 
in  the  region  2  flow.  From  the  results  of  these  tests  it  is  believed  that 
this  gauge  will  be  quite  useful  for  detailed  model  studies.  For  example,  the 
response  time  of  a  non-shock  mounted  gauge  in  the  described  flow  situation  is 
4  microseconds  and  the  ringing  envelope  is  about  +  5^  of  the  mean  value 
(p  *  200  psi).  The  mean  value  remains  constant  within  3^  until  the  reflected 
wave  appears  at  the  gauge  (200  microseconds).  Figure  11  is  a  photograph  of 


some  of  the  pressure  Instrumentation  being  used.  To  the  left  ia  the  regular 
Kistler  gauge  and  in  the  center  is  a  gauge  manufactured  by  the  Atlantic 
Research  Corporation.  This  instrument  has  a  barium  titanate  crystal  imbedded 
in  an  epoxy  resin.  The  miniature  Kistler  gauge  is  shown  to  the  right  next  to 
the  2-inch  hemisphere-cylinder  model  in  which  it  has  been  mounted. 

The  phctographic  and  optical  equipment  Includes  a  double  (horizontal 
and  vertical)  10-inch  Schlieren  system  in  conjunction  with  various  phctographic 
techniques.  In  addition  to  still  sub"-mlcrosecond  spark  and  Fastax  equipment, 
there  is  available  a  unique  rotating  mirror,  rotating  drum,  framing  camera  that 
allows  the  taking  of  sub-microsecond  exposures  every  few  microseconds  for  a 
continuous  period  of  a  few  milliseconds.  Long  duration  (l  to  3  milliseconds) 
exploding  wires  are  used  as  a  light  source  for  this  camera.  An  electronic 
shutter  (Kerr  cell)  is  also  available  and  can  be  used  for  still  photographs  and 
with  a  simple  drum  camera. 

It  has  been  mentioned  that  the  aerodynamic  properties  of  the  flew  In 
the  three  test  regions  are  complimentary  to  each  other  in  terms  of  the  flow 
Mach  and  Reynolds  numbers.  The  maximum  performance  that  can  be  achieved  in 
each  region  will  new  be  presented  in  terms  of  the  simulation  of  state  proper¬ 
ties  in  the  stagnation  region  of  a  blunt  body.  These  calculations  are  based 
upon  the  properties  of  air  in  thermodynamic  equilibrium  and  upon  a  10,000-psi 
driver  produced  by  the  constant  volume  combustion  of  a  7^  helium,  3 stoi¬ 
chiometric  hydrogen  and  oxygen  mixture. 

In  the  region  2  flow,  stagnation  region  properties  at  the  nose  of  a 
blunt  model  can  be  produced  to  match  those  on  a  free  flight  vehicle  traveling 
at  a  Mach  number  of  15  at  sea  level.  The  maximum  simulated  flight  Mach  number 
increases  with  altitude  until  at  130,000  feet  an  of  30  is  simulated.  Cor¬ 
responding  calculations  for  the  non-ref lected  nozzle  (A  /A2  =  25)  show 

simulation  at  **5,000  feet  for  en  M.  of  16  and  at  90*000  feet  for  an  Mf 

of  30.  For  the  reflected  nozzle  (A  /A  =  l6C0),  Mach  17  is  simulated 
at  120,000  feet  and  Mach  30  at  m  160,000  feet.  It  is  again  pointed 

out  that  the  test  section  flow  Mach  number  in  each  of  these  cases  is,  in  gen¬ 
eral,  different  ard  usually  lower  than  the  simulated  Mach  number.  The  useful¬ 
ness  of  the  nozzle  flows  is  seen,  therefore,  in  providing  a  considerable  Mach 
huE&er  range  for  an  experimental  program. 

The  thermodynamic  and  chemical  state  cf  the  gas  in  the  test  section 
is  an  important  characteristic  cf  the  nozzle  flows.  The  expanding  gas  will 
be  somewhere  between  the  frozen  (at  upstream  concentrations)  and  the  equili¬ 
brium  state  depending  upon  the  flow  times,  or  size  and  design  of  the  facility, 
and  upon  the  recombination  kinetics.  Calculations  have  been  made  fer  several 
different  flow  models  considering  air  as  the  expanding  gas  and  for  the  re¬ 
flected  and  non-ref  lected  nozzles.  Initial  conditions  correspond  to  a  shock 
Mach  number  of  12  and  an  initial  driven  tube  pressure  of  100  mm  of  Hg.  The 
flew  models  are  described  and  the  results  of  the  calculations  are  given  in 
Table  I.  It  is  seen  that  the  stagnation  pressure  at  the  model  nose  does  not 
differ  largely  between  the  different  flow  models.  Therefore,  the  Reynolds 
number  of  the  flow  are  not  greatly  affected  by  departures  from  equilibrium. 

The  flow  Mach  numbers  do  not  vary  from  the  equilibrium  value  by  more  than  lOJt 
except  for  the  last  frozen  model  in  the  reflected  nozzle  calculations.  At- 


-  32  - 


tempts  will  be  made  to  determine  the  departure  from  equilibrium  in  the  nozzle 
flows  by  measuring  the  stagnation  pressures  and  upstream  pressures  during  a 
run. 


The  test  program  being  followed  can  be  divided  into  three  somewhat 
distinct  phases,  the  first  two  of  which  are  preliminary.  The  initial  test 
series  was  a  study  of  the  combustion  driver  under  conditions  in  which  it  was 
safe  to  investigate  high  pressure  phenomena.  Experiments  were  run  with  the 
short  shock  tube  in  a  rocket  test  pit  at  the  General  Electric  Malta  Test 
Station.  The  test  installation  is  shown  in  Figure  12. 

There  were  f.hree  general  objectives  of  this  program.  The  first  was  to 
determine  the  test  conditions  required  to  produce  repeatable  combustion  over  a 
pressure  range  of  approximately  500  to  10,000  psi,  and  to  match  the  combustion 
pressure  rise  with  a  cleanly  breaking  diaphragm  near  peak  pressure.  The  second 
was  to  determine  the  hazard  limits  of  the  facility  such  as  the  detonation  limits 
and  the  ability  of  the  metal  diaphragms  to  stay  in  one  piece  after  rupturing. 

The  third  objective  was  to  become  familiar  with  the  operation  of  the  facility. 

The  first  runs  were  at  constant  volume  conditions  and  were  done  tc 
determine  the  pressure  history  in  the  driver  throughout  the  combustion- cycle. 
Promoter-type  spark  plugs  were  located  along  the  inner  tube  surface  in  a 
spiral  array  at  distances  varied  for  different  tests  from  1  to  1-1/2  feet. 

The  pressures  were  sensed  at  two  stations  in  the  tube  by  standard  Kistler 
gauges  and  were  recorded  on  cathode  ray  oscilloscopes.  Some  of  the  results  are 
shown  in  Figure  13.  The  maximum  pressure  rise  is  generally  about  90ft  of  the 
theoretical  value  which  is  calculated  on  the  basis  of  r.o  dissociation  and  no 
heat  loss  to  the  surroundings.  The  heat  loss  to  the  cold  tube  walls  is  respon¬ 
sible  for  the  drop  in  pressure  after  the  peak  value  is  reached. 

The  results  shown  in  Figure  13  are  considered  to  be  good.  To  obtain 
these  pressure  rises  it  was  necessary  to  provide  good  gas  mixing  (multiple 
injection  ports  through  which  the  gases  were  swirled  into  the  driver  tube), 
variable  ignition  energy  for  certain  pressure  ranges  to  retard  detonation 
(not  critical  below  initial  pressures  of  300  psi,  low  ignition  energy  at  and 
above  300  psi),  variable  ignitor  plug  locations  (l-foot  spacing  at  and  below 
300  psi,  l-l/2-foot  spacing  above  300  psi). 

Figure  lh  presents  pressure  data  obtained  when  detonation  or  auto¬ 
ignition  occurs  during  -onstant  volume  combustion.  This  is  undesirable  from 
the  points  of  view  of  both  safety  and  tube  performance.  The  100-psi,  50)1 
helium  data  is  not  strictly  detonation,  but  6trong  waves  are  produced  in  the 
tube.  The  delayed  detonation  at  80]t  helium  is  Interesting  and  was  repeatable. 
Because  of  its  general  occurrence  and  because  of  the  occurrence  of  strong 
detonation  at  lower  helium  percentages,  the  preferred  operating  percentages 
are  limited  to  704  to  754  helium  (except  at  =  100  and  200  psi  when  it  is 
possible  to  combust  cleanly  with  6056  to  65^  helium  mixtures). 

Pressure  records  obtained  from  the  driver  when  combustion  is  matched 
with  breaking  diaphragms  are  shown  in  Figure  15.  The  first  record  shows  a 
typical  response  for  a  detonable  mixture.  The  trace  for  *  900  psi  is 
distorted  somewhat  by  extraneous  electronic  noise.  For  approximately 
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75  of  the  runs  It  was  possible  to  obtain  good  matching  (diaphragm  breaks 
within  5$  of  the  peak  combustion  pressure). 

Four  different  diaphragm  materials  were  tested  -  annealed  copper, 

2hS0  aluminum,  SAE  1020  steel,  and  type  3CU  stainless  steel.  Figure  16  is  a 
photograph  of  copper  diaphragms  ruptured  in  the  circular  support  insert.  The 
thicknesses  shown  are  1/8,  1/4,  and  3/8  inch  and  the  diameter  is  l6  inches. 
The  diaphragms  have  a  milled  cross  on  one  face  so  that  a  controlled  rupture 
is  obtained.  The  depth  of  the  cut  is  kept  between  lOjt  and  20)t  of  the  thick¬ 
ness  so  that  the  required  diaphragms  thicknesses  are  not  excessive.  The  most 
desirable  diaphragm  materials  have  been  found  to  be  copper  to  3000  psi  and 
stainless  steel  from  3000  to  10,000  psi. 

The  objectives  of  the  second  phase  of  the  experimental  program,  which 
is  being  run  now,  are  to  calibrate  the  complete  facility  and  to  evaluate  the 
model  instrumentation  techniques.  Certain  aspects  of  this  work  will,  of 
course,  continue  for  some  time. 

The  third  phase,  which  involves  detailed  experimental  programs,  is 
about  to  begin.  The  following  list  presents  some  general  aerodynamic  studies 
for  which  the  facility  is  suited. 

1.  Laminar  boundary  layer  heat  transfer 

2.  Turbulent  boundary  layer  heat  transfer 

3.  Spearated  flow  heat  transfer 

4.  Body  pressure  distributions 

5-  Flow  field  characteristics 


ACKNOWLEDGEMENT 

The  construction  and  operation  of  this  facility  within  the  Missiles 
and  Ordnance  Systems  Department  was  sponsored  by  the  Air  Force  under 
Contract  AF  -04 (6U5) -24 .  The  author  gratefully  acknowledges  this  sponsorship. 

A  major  portion  of  the  design  of  the  facility  and  of  the  instrumenta¬ 
tion  was  originated  by  Dr.  Y.  A.  Yoler  of  the  Missiles  and  Ordnance  Systems 
Department. 


Table 


COMPARISON  OF  TSCT  SECTION  AIR  PROPERTIES  WITH 
FROZEN  AND  EQUILIBRIUM  EXPANDING  FLOW  MODELS 
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REGIONS  AVAILABLE  FOR  AERODYNAMIC  EXPERIMENTATION 
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Figure  2.  Regions  Available  for  Aerodynamic  Experimentation 
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SCHEMATIC  DRAWING  OF  SIX  INCH  SHOCK  TUNNEL  8  SHOCK  TUBE  FACILITIES 

M.O.  S.D.  AERODYNAMIC  LABORATORY 

Figure  4.  Schematic  Drawing  of  Six  Inch  Shock  Tunnel  and  Shock  Tube 
Facilities  -  MOSD  Aerodynamic  Laboratory 
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Figure  6.  MOSD  Shock  Tunnel  -  Driver  end 


Diaphragm  Clamping  Mechanism 
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Figure  9.  Second  Diaphragm  Clamping  Mechanism,  Nossle,  and  Teat  Section 


FIGUPE  5C 

STINC  SUPPORT  AND  RETRACTION  MECHANISM 

igu'-c-  10.  Sting  Support  and  Retraction  Mechanism 
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Combustion  Studies 


TYPICAL  CONSTANT  VOLUME  COMBUSTION  PRESSURE 


RECORDS  OVER  OPERATION  RANGE  OP  DRIVER 


CONSTANT  VOLUME  COMBUSTION  PRESSURE  RECOROS 
WHEN  DETONATION  OCCURS 
MALTA  TEST  PROGRAM,  OCT.  19  56 
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Figure  14.  Constant  Volume  Combustion  Pressure  Records  when 
Detonation  Occnrs 


COMBUSTION  PRESSURE  RECOROS  IN  DRIVER  TUBE 


WITH  GHhAKING  DIAPHRAGMS 


MALTA  TEST  PROGRAM,  OCT  1 9  5  6 
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Figure  16.  16  Inch  Diameter  Copper  Diaphragm  1/8,  1/4,  3/8  Inch  Thick 


DEVELOPMENT  CF  A  SHOCK  TUBE  TO 
GENERATE  VARIABLE  STEPPED  SHOCK  FRONTS 

Richard  I.  Condit  and  Robert  V.  Steele* 

Broadview  Research  and  Development 

As  reported  in  preliminary  form  at  the  December  1953  meeting  of  the 
American  Physical  Society  at  Stanford  University,  California,  we  have  had  some 
success  in  replacing  the  bursting  diaphragm  of  conventional  shock  tubes  with 
electromechanical  devices.  These  so-called  "shock  wave  valves"  allow  the  gas 
in  the  compression  chamber  to  be  dumped  repeatedly  with  merely  the  push  of  a 
button  and  with  no  intervening  requirement  for  disassembly  or,  replacement  of 
parts.  A  sketch  illustrating  the  operating  features  of  one  of  these  devices 
is  shown  in  Figure  1.  Initially  the  same  pressure  is  placed  in  both  the  com¬ 
pression  chamber  and  the  actuating  chamber.  The  resulting  pressure  differen¬ 
tial  across  the  central  part  of  the  large  diaphragm  seals  off  the  compression 
chamber  from  the  expansion  chamber.  Operation  of  the  valve  occurs  when  the 
quick-release  solenoid  is  suddenly  tripped,  thereby  releasing  the  pressure  in 
the  actuating  chamber  and  causing  the  large  flexible  diaphragm  to  be  thrown 
open  to  release  the  gas  in  the  compression  chamber  into  the  expansion  chamber. 
If  the  rate  of  appearan'e  of  gas' in  the  expansion  chamber  is  suff  icientl^\ 
fast,  the  compression  wave  that  results  will  develop  into  a  shock  wave  as  it 
moves  down  the  tube.  Thus  shock  Waves  can  be  .conveniently  produced  without 
the  inconvenience  of  disassembly  after  each  shot,  or  flying  fragments  from 
broken  diaphragms.  (Incidentally  large  arrays  of  such  shock  wave  valves 
would  seem  to  be  especially  attractive  for  giant  shock  tubes,  where  both  dis¬ 
assembly  and  diaphragm  puncture  tend  to  become  either  difficult  or  tedious.) 

In  searching  for  a  useful  application  for  such  devices,  we  were  in¬ 
formed  that  there  was  a  need  for  a  laboratory  arrangement  which  could  produce 
a  variable  stepped  shock  front  whose  pressure  characteristics  were  as  follows: 


final  pressure!  5  “  20  psi;  total  duration:  20  milliseconds 

bt  of  1st  step:  25  -  75$  of  final  pressure  duration  of  1st  step:  ^illSseeonds 

*  Bow  with  Pratt  and  Whitney  Co.,  Inc.)  Livermore,  California 
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It  Is  to  be  noted  that  the  intended  use  was  in  calibrating  pressure  gauges,  so 
that  the  pressure  profile  was  all  important,  and  the  conditions  of  gas  flow  of 
little  importance.  Consequently  in  our  development  efforts  we  have  considered 
only  pressure-time  characteristics  and  have  ignored  the  associated  gas  flov 
phenomena . 


It  seemed  to  us  that  here  was  an  interesting  opportunity  for  the  use 
of  Pho-'k  wave  valves,  and  we  therefore  proposed  to  develop  an  arrangement  to 
produce  the  desired  pressure  profile  using  two  separate  shock  wave  valves. 

Their  times  of  actuation  were  to  be  so  coordinated  as  to  dump  first  one  and  then 
the  other  into  the  extension  chamber.  The  variable  duration  of  the  first  step 
was  to  be  provided  by  the  spacing  between  the  times  of  actuation  of  the  two 
units.  The  variable  heights  of  the  two  steps  were  to  be  provided  by  the  dif¬ 
ferent  pressures  placed  in  the  separate  compression  chambers.  On  this  basis 
it  was  decided  to  proceed  with  the  development  of  a  generator  of  variable 
stepped  shock  fronts.  The  size  of  the  installation  was  to  be  such  as  to 
supply  the  desired  stepped  wave  to  a  shock  tube  having  an  inside  diameter  of 
about  4-1/2  inches.  The  scheme  is  illustrated  below. 


shock  wave 
source  #1 


shock  wave  . 

source  #2  — 

Time  of  actuation  of  1  with  respect  to  §2  is  to  be  variable  and  such  as  to 
produce  the  desired  stepped  shock  front. 

Early  efforts  so'n  demonstrated  that  a  single  shock  wave  valve  could 
be  readily  constructed  to  produce  satisfactory  shock  waves  for  a  4-1/2  inch 
tube  in  the  pressure  range  of  interest.  However  the  possibility  of  readily 
synchronizing  the  tine  of  operation  of  two  such  shock  wave  valves  to  the 
desired  accuracy  turned  out  to  te  very  remote  indeed.  Since  the  shortest 
separation  in  time  of  the  two  steps  was  to  be  2  milliseconds,  this  meant  that 
the  time  of  operation  of  each  shock  wave  valve  had  to  be  reproducible  to  a 
small  fraction  of  2  milliseconds.  Thus  one  requires  the  two  electromechanical- 
pneumatic  assemblies  to  maintain  a  reproducible  relative  time  of  operation  of 
something  like  100  microseconds.  Such  a  high  degree  of  accuracy  in  the  time 
of  operation  of  the  simple  devices  constituting  our  shock  wave  valves  was  soon 
shown  to  be  several  orders  of  magnitude  away  from  that  which  was  readily 
attainable.  Rather  than  attempt  to  carry  out  the  development  along  these 
lines  with  such  extreme  requirements  for  improvement,  it  seemed  better  to  re¬ 
consider,  and  to  search  for  some  other  principle  for  generating  stepped  shock 
fronts,  some  principle  that  would  seem  to  be  more  naturally  compatible  with 
the  specified  requirements.  Thus  part  way  through  the  development  it  seemed 
necessary  to  abandon  this  procedure  and  start  over  again. 


In  our  uminations  about  what  to  try  next  we  gradually  became  con¬ 
scious  of  the  significant  apparent  difference  in  the  timing  requirements  if 
the  base  could  be  changed  from  time  to  distance.  When  1  millisecond  or  so 
is  viewed  as  the  time  interval  between  the  operation  of  electromechanical 
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devices  it  appears  as  a  7ery  severe  requirement.  However,  since  in  1  milli¬ 
second  a  shock  wave  travels  roughly  1  foot,  the  accurate  positioning  which 
is  readily  possible  with  respect  to  1  foot  suggests  that  a  fine  equivalent 
breakdown  in  time  should  be  easy  if  it  could  be  based  on  the  distance  travelled 
by  the  shock  wave. 

Pushing  on  this  favorable  thought,  it  seemed  that  if  the  two  compres¬ 
sion  chambers  for  the  stepped  shock  front  were  placed  at  different  distances, 
then  they  could  be  initiated  simultaneously  and  still  get  the  desired  output 
pressure  profile.  This  scheme  is  illustrated  below. 


Time  of  actuation  of  #1  and  §2  to  be  simultaneous;  their 
variable  separation  in  distance  is  to  provide  the  desired 
stepped  shock  front. 


Although  this  scheme  was  more  attractive  than  the  original  one  It  3till  has  a 
significant  requirement  for  synchronized  timing,  since  two  operations  must  be 
done  simultaneously.  Clearly  an  improvement  would  seem  to  follow  if  all  re¬ 
quirements  for  synchronize!  operations  were  eliminated.  This  would  seem  to 
necessitate  the  use  of  but  one  source  of  shock  waves.  Thus  the  need  was  for 
an  arrangement  employing  but  one  source  of  shock  waves  (and  therefore  no 
requirements  for  synchronization)  and  a  variable  inferential  path  length 
(to  allow  the  necessary  separation  in  time  of  the  stepped  shock  wave). 


Such  an  arrangement  is  possible  with  the  scheme  illustrated  below. 


Thus  all  requirements  would  seem  to  he  met  in  principle  by  using  but  one  source 
of  shock  waves  and  two  branch  ducting  between  the  source  and  the  measurement 
section.  The  difference  in  pa*h  length  could  be  made  variable  by  extending  the 
longer  branch  in  a  manner  similar  to  that  of  a  trombone.  (Because  of  this 
similarity  this  shock  wave  generator  ras  been  referred  to  as  a  "tromrone". ) 
Actually,  th£  variable  length  of  hue  first  s*ep  in  the  pressure  profile  Is 
provided  by  an  array  cf  flanged  pipes  of  different  lengths  which  can  be  in¬ 
serted  into  the  longer  branch.  This  procedure  of  a  series  of  fixed  steps  in 
the  longer  branch  was  employed  ratoer  tnan  a  con* nuously  variable  sliding 
arrangement  for  simplicity,  to  provide  a  smooth  inferior  wall,  and  to  more 
readily  withstand  the  20-psi  sb.ock  forces. 

For  the  single  source  cf  shock  waves,  it  was  decided  to  utilize  a 
conventional  shock  tube  compression  chapter  and  frangible  diaphragm  assembly 
since  there  were  no  longer  any  requirements  for  synchronized  times  of  operation 
and  since  such  a  source  is  probably  superior  for  calibration  purposes.  The 
general  arrangement  of  the  resulting  shock  wave  generator  is  shown  in  Figure  2. 

Fig.  2  shows  the  assembly  for  a  difference  of  path  length  of  8  feet 
(two  4-foot  pipes).  Fa*h  lengtn  differences  of  4  feet  (two  2-foot  pipes)  and 
12  feet  (two  2-foot  plus  k-fo c*-  pipe  combinations)  were  also  provided.  Tne 
compression  chamber  pipe  was  made  sufficiently  long  to  provide  the  specified 
20  milliseconds  total  duration.  The  diaphragm  and  diaplu*agm-punctur  ing 
mechanism  were  conventional  and  located  just  ahead  of  the  "Y"  branching. 

To  control  the  distribution  of  pressure  between  the  two  ducts,  a 
hallow  tube  called  a  "scoop"  was  placed  wi*hin  each  passageway  of  the  inpu* 

"Y".  Each  of  these  was  moveable  as  required  from  1,he  outside.  To  tuild  up 
the  pressure  of  the  first  step  the  scoor>  in  the  leading  duct  was  pushed  for¬ 
ward.  Similarly  to  build  up  the  press^e  in  the  second  step  (and  reduce  the 
pressure  in  the  firs*,  step)  *h@  scoop  In  the  leading  duct  would  be  withdrawn 
and  the  scoop  in  the  trailing  duct  advanced.  Only  one  scoop  was  effective  at 
a  time;  the  scoop  in  *he  trailing  duct  served  for  first  step  pressures  from 
C fjo  to  50/G  of  the  total  pressure,  while  *hat  in  The  leading  due*  served  for 
first  step  pressures  from  to  100i  of  the  total. 

Samples  of  oscillograms  of  pressure  versus  time  obtained  with  this 
apparatus  are  shown  in  Figure  3-  Tee  upper  series  shows  how  tee  duration  of 
the  first  s+ep  gets  grea'er  as  the  dif ferer.tlal  path  length  is  increased.  The 
lower  series  shews  how  *he  relative  height  of  *te  first  s  +  ep  decreases  as  the 
scoop  in  the  leading  duct  is  withdrawn  and  that  in  the  trailing  duct  is  .ad¬ 
vanced  . 


From  an  extended  series  of  such  data  we  can  obtain  the  pressure  height 
of  the  first  step  as  a  function  cf  the  total  pressure,  for  various  positions 
of  the  scoops.  These  relations  are  shown  in  Figure  4.  in  Figure  4,  the  path 
length  difference  has  been  held  constant  at  8  feet.  -The  fact  that  each  line 
is  straight  reveals  that  the  distribution  of  gas  between  the  two  ducts  is 
fairly  independent  of  the  gas  pressure  (within  the  limits  of  the  available 
data).  The  spread  among  the  lines  shows  how  the  total  pressure  can  be  dis¬ 
tributed  at  will  among  the  two  ducts  by  adjustments  in  scoop  positions. 
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While  the  main  effect  of  the  differential  path  length  is  to  determine 
the  duration  of  the  first  step,  and  the  main  effect  of  scoop  position  is  to 
determine  the  relative  height  of  the  fir6t  to  the  second  step,  there  is  ac¬ 
tually  some  interaction  between  these  two  effects.  This  interaction  of  scoop 
position  on  the  duration  of  the  first  step  stems  from  the  variations  in  shock 
wave  velocity  with  pressure  level.  If  all  shock  waves  progressed  through  the 
trombone  passages  at  a  constant  velocity,  then  the  height  and  duration  of  the 
first  step  would  be  strictly  independent.  Actually,  as  the  pressure  rises, 
the  shock  wave  velocity  increases  somewhat,  so  that  the  same  lifferential  path 
length  will  lead  to  first  steps  having  slightly  different  durations. 

These  interactions  between  length  and  height  of  the  first  step  are 
documented  in  Figure  5*  Figure  5  gives  the  duration  of  the  first  step  as  a 
function  of  path  length  difference  for  various  scoop  positions  and  several 
different  total  pressures  (compression  chamber  pressures).  It  can  be  seen 
that  the  general  average  duration  decreases  as  the  pressure  is  raised  (arrow 
positions  drop  with  increasing  pressure).  Additionally,  the  spread  in  dura¬ 
tion  for  different  scoop  settings  tends  to  increase  as  pressures  get  higher. 
Ths-se  interactions  complicate  somewhat  the  adjustment  of  the  apparatus  but  in 
no  way  prevent  it  from  satisfying  the  given  specifications. 

As  described,  the  trombone  was  found  to  yield  satisfactory  stepped 
shock  fronts  over  the  specified  range  of  interest. 


However  with  this  goal  achieved,  it  was  determined  that  a  "slight" 
deviation  from  the  original  specifications  was  highly  desirable.  To  wit,  the 
first  step  should  have  a  slowly-rising  compression  wave  for  the  leading  edge 
rather  than  a  discontinuous  shock  wave.  However  the  leading  edge  of  the 
second  step  should  remain  a  shock  wave.  Thus  the  desired  pressure  profile  was 


changed  from  the  original 


to  the  new 


was  not  yet  finished. 


So  the  development  effort 


Since  all  other  specifications  seemed  to  be  met  by  the  trombone  @8 
described,  it  seemed  desirable  to  attempt  to  supply  the  additional  feature  of 
a  slowly  rising  front  by  adding  some  new  component  in  such  a  way  that  the 
favorable  results  already  achieved  were  not  sacrificed.  Therefore  the  initial 
attempt  (and  the  approach  actually  adopted)  was  to  leave  the  trombone  as  it 
was- -arranged  to  produce  two  shock  vaves--and  then  to  add  something  which 
would  degrade  the  leading  shock  wave  into  a  compression  wave.  It  seemed  con¬ 
ceivable  that  the  leading  shock  wave  could  be  degraded  into  a  compression 
wave  by  causing  it  to  open  some  kind  of  a  door  which  initially  blocked  the 
shock  tube.  With  such  arrangement,  the  shock  might  leak  past  Ju.-.t  a  little  at 
a  time  as  the  door  Just  opened,  and  then  later  it  might  be  nearly  unimpeded 
after  the  door  was  wide  open.  It  was  hoped  that  the  impedance  to  gas  flow 
presented  by  such  an  opening  door  would  yield  a  compression  wave  when  acted 
on  by  the  incident  shock  wave.  However  the  second  shock  front  would  pre¬ 
sumably  pass  through  the  doorway  unmodified,  since  the  door  would  already  be 
opened  when  it  arrived. 


Schematically,  the  leading  shock  wave  1b  forced  to  open  up  a  swinging 
door  Bomevhat  as  follows: 

Before  1st  shock  strikes 
swinging  door 

2nd  1st  closed 

shock  shock  swinging 

front  front  door  % 


2nd  1st  front 

shock  leaking  through 

front  partially  open  ioor 


As  1st  shock  strikes 
swinging  door 


2nd  1st 

shock  compression 
front  front 


Door  held  open  by  gas 
flow  as  2nd  shock  arrives 


swinging 

door 


2nd  shock  front  passes 
through  open  door 
unchanged 
shock  compression 
front  front 


Actually,  with  a  tube  of  circular  cross  section  It  Is  not  convenient 
to  hinge  the  swinging  door  at  the  vail,  so  a  central  post  was  added  as  a  sup¬ 
port  and  pivot  point.  The  first  swinging  doors  for  the  trombone  were  made  of 
metal  door  hlngles. 


Unfortunately  these  metal  assemblies  could  not  be  made  to  withstand  more  than 
a  few  shocks.  After  that  they  were  generally  damaged  beyond  repair.  However, 
rubber  diaphragms  clamped  between  two  central  posts  would  survive  many  shocks 
and  when  properly  arranged  would  produce  the  desired  compression  wavo. 
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The  proper  arrangement  of  these  "flappers"  consisted  of  a  suitable  ratio  of 
mass  to  restoring  force  for  a  given  region  of  operating  pressures.  If  the  mass 
was  too  small  for  the  given  shock  pressure,  the  shc~k  front  would  pass  through 
the  flapper  assembly  with  scarcely  any  degradation.  If  the  restoring  force  of 
the  rubber  were  too  great  the  flapper  would  partially  close  against  the  flowing 
gas  and  the  resulting  constriction  would  show  up  as  a  dip  in  the  pressure  pro¬ 
file  in  place  of  the  desired  flat  top. 

Figure  6  shows  two  oscillograms  of  the  stepped  fronts  from  the  trombone 
when  the  first  step  was  degraded  to  a  compression  wave.  Correct  and  incorrect 
operation  of  the  flapper  is  6hown. 

With  a  suitable  combination  of  metal  weights  added  to  various  rubber 
diaphragms,  it  was  possible  by  empirical  trials  to  evolve  a  family  of  flappers 
which  would  degrade  the  leading  edge  of  the  first  shock  to  a  compression  wave 
over  the  entire  pressure  range  of  interest. 

The  only  apparent  drawback  to  this  procedure  was  that  the  majority  of 
flappers  were  not  self-cocking.  After  each  employment  they  had  to  be  reset  by 
hand,  a  process  requiring  disassembly  of  that  part  of  the  trombone. 

In  a  final  development  effort  some  success  was  achieved  with  rubber 
diaphragms  cut  in  the  form  of  a  star  and  mounted  between  flanges.  The  arrange¬ 
ment  of  these  rubber  "stars"  is  shown  in  Figure  7»  These  devices  tended  to 
reset  themselves  automatically,  so  they  could  be  used  repeatedly  without  dis¬ 
assembly.  Here  again  it  was  necessary  by  trial  to  determine  the  proper  ratio 
of  mass-to-restoring  force  for  the  petals  of  these  star-cut  diaphragms.  Mass 
was  increased  by  adding  metal -weights  to  the  tips  of  the  petals.  Restoring 
force  was  reduced  by  undercutting  the  rvibber  In  the  region  at  the  base  of  the 
petals  where  the  extreme  flexing  was  to  be  expected.  In  the  limited  time 
available  for  these  finishing  touches  it  was  not  possible  to  achieve  re¬ 
cocking  star  assemblies  for  the  entire  pressure  range  of  interest.  However 
a  significant  replacement  of  flappers  by  stars  was  accomplished. 

It  Is  to  he  noted  that  the  compression  wave  in  the  leading  edge  of  the 
first  step  following  the  degradation  by  the  star  or  flapper  assembly  is  only 
transitory,  since  it  tends  to  build  back  up  into  a  steep-fronted  shock  wave  as 
it  travels  on  down  the  expansion  chamber.  Because  of  these  inevitable  changes, 
as  long  as  a  compression  wave  is  required,  it  is  on~y  possible  to  pro  uce  >4 
desired  pressure' prof  lie  at  one  distance  from  where  degradation  occurs.  In  our 
set-up,  the  point  of  observation,  i.e.  the  point  for  pressure  gauge  calibration, 
was  to  be  about  10  feet  beyond  the  position  of  the  star  or  flapper  assembly. 

At  distances  closer  than  1C  feet,  the  compression  wave  would  have  a  slower  rise 
time  and  at  greater  distances  a  faster  rise  time  than  at  the  designated  point 


of  observation.  j 

The  final  array  of  stars  and  flappers  which  were  supplied  to  provide  the 
degraded  leading  edge  for  the  stepped  shock  front  throughout  the  pressure  range 
of  interest  is  shown  in  Figure  8.  It  can  be  seen  that  the  range  of  effective¬ 
ness  of  any  ’one  degrader  Is  not  extensive  and  consequently  a  sizeable  number  of 
different  units  are  required. 

Thus  the  specifications  created  for  the  variable  stepped  shock  front 
were  finally  satisfied  by  this  development.  It  is  interesting  to  note  that 
the  resulting  equipment  was  a  far  cry  from  the  pair  of  shock  wave  valves  ori¬ 
ginally  envisioned. 
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Figure  2.  General  Arrangement  of  Special  Shock  Tube  to 
Generate  Variable  Stepped  Shock  Fronts 


Figure  3*  Variation  of  1st  Step  Duration  by  Changing  Differential  Spool  Lengths 

Main  Pressure  Constant  at  v.o  pal  Leading  Duet  Scoop  Fixed  at  2-5/8 lnche* 
(Dual  beam  oscillograms  of  pressure  vs  time  are  shown;  two  gauges  per  shot) 
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Figure  4.  Variation  of  l3t  Step  Pressure  as  a  Function  of  Scoop  Setting 
at  Different  b'aln  Pressures 


note  undesirable 
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Prenature  Closing 


shock  wave  front 


compression  wave  front 


Correct  Opening 


Figure  6.  Examples  of  Prenature  Closing  and  Correct  Opening 

of  Slope  Generator 


2  4  6  8  10  15  pal 

1st  Step  Pressure 

Figure  8.  Bar  Graph  of  Slope  Generator  Assemblies  for  Different  Pressures 


SjSjglgS  OF  TRAHSIENT  AIR  FORCES 
tfj  IV  0  -  L  IMn  NS  I  OKA  I.  AIRFOILS* 

Frunett  A.  Witmer 

Massachusetts  Institute  of  Technology 


Th«  shock  tube  facility  at  M.I.T.  was  designed  to  serve  a  three-fold 
purposes  (a)  to  study  blast  airloads  on  aircraft,  (b)  to  Investigate  struc¬ 
tural  response  to  such  airloads,  and  (c)  to  conduct  basic  fluid  mechanics 
studies.  Ac-ordlngly,  the  design  of  this  facility  was  preceded  by  a  study  of 
the  aerodynamics  areas  to  be  investigated  and  appropriate  measurement  tech¬ 
niques  for  each  experiment  study.  For  studies  of  two-dimensional  airflow 
development  (about  airfoils,  simple  shapes,  the  wall  boundary  layer,  and  so 
forth),  the  interferometer  ar peered  to  offer  many  advantageous  features.  For 
bcHjiuary-]ayer  work,  a  narrow  tube  was  Indicated.  Sting  mounting  of  two- 
dimensional  models  was  rejected  because  of  the  associated  optical  and  flow 
interference  effects.  Clamping  these  models  between  the  test  section  windows  . 
themselves  was  the  selected  technique;  this  technique.  In  turn,  Indicated  the 
need  for  a  narrow  tube.  In  addition,  a  very  deep  tube  was  Indicated  In  order 
to  reduce  the  interference  effects  from  the  top  and  bottom  walls.  In  the 
final  analysis,  space  limitations  for  the  facility  limited  the  tube  depth. 
Studies  of  air  forces  on  finite-span  wings  and  of  response  of  structural 
models  to  blast-induced  air  forces  requires  a  tube  of  relative]/  large  cross- 
section  in  o’der  to  minimise  respectively  aerodynamic  Interference  effects 
and  the  problems  of  bui'  Slr.g  suitable  structural -response  models.  The  result 
of  these  and  other  similar  considerations  is  the  8  by  2h-lnch  croS3  section 
of  this  tube. 

The  requirement  that  transient  air  forces  be  measured  from  the  Instant 
of  shock  impingement  until  steady  state  air  forces  were  achieved  on  the  model 
in  the  (theoretically)  uniform  hot  gas  flow  determined  the  length  of  the 
shock  tube.  Theory^  predicts  that  for  low  Mach  number  flows,  the  30ught-for 
steady  state  air  force  condition  is  essentially  achieved  In  roughly  twelve- 
chordlc-ngths  of  flow.  Taking  a  h -inch  chord  airfoil  and,  rather  arbitrarily, 
a  90  fcot-per-se cond  material  velocity  condition  as  a  basis,  the  shock  tube 
length(s)  resulting  are  as  given  in  Figure  1.  The  selection  of  a  smaller 


*  Since  the  group  attending  the  symposium  visited  the  M.I.T.-W .A.D.C.  shock 
tube  facility  at  which  time  a  number  of  shock  tube  features  and  experiments 
were  discussed,  a  number  of  items  were  deleted  from  this  paper  as  presented. 
The  chitted  uaterlal  is  included  here,  briefly,  primarily  for  those  readers 
who  were  unable  to  attend  the  symposium. 
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material  velocity  (to  meet  this  flow-length  criterion)  would  have  greatly 
increased  the  required  length  of  the  chock  tube.  It  should  be  pointed  out  that 
the  Bteady  state  air  force  condition  is  achieved  in  even  fewer  flow  chord 
lengths  for  the  higher  subsonic  ana  Tor  supersonic  flow  Mach  numbers,  excepting 
the  region  near  Mach  1.  Figure  2  shows  the  theoretically  available  flow  dura¬ 
tion  in  terms  of  flow  chord  lengths  for  a  **-inch  chord  model.  The  flow  duration 
at  the  test  section  is  limited  at  low  Mach  numbers  by  the  arrival  of  the  shock 
reflected  from  the  closed  end  of  the  expansion  chamber.  At  higher  Mach  numbers, 
the  arrival  of  the  contact  surface  (or  region)  terminates  the  hot  gas  flow  at 
the  test  section.  In  practice,  it  is  found  that  the  contract  surface  arrives 
earlier  than  predicted  by  ideal  shock  tube  theory-3'  **»  5.  For  the  present  tube 
conf iguration,  this  earlier  contact  surface  arrival  is  indicated  by  the  experi¬ 
mental  data  in  Figure  2p  air  is  the  working  medium  in  all  cases  shown. 

For  aerodynamic  studies,  it  is  desirable  that  measurements  be  made  at 
the  largest  feasible  Reynolds  number  in  order  that  extrapolation  to  free  flight 
conditions  of  interest  can  be  carried  out  with  better  reliability.  With  a 
model  of  fixed  chord  and  U6ing  (room)  air  as  a  working  medium,  these  higher 
Reynolds  numbers  can  be  achieved  most  conveniently  by  an  increase  of  the  flow 
density.  This,  in  turn,  means  pressurizing  the  expansion  chamber  of  the  shock 
tube.  At  the  250-psi  a  maximum  design  pressure  selected  for  this  shock  tube, 
the  available  flow  Reynolds  numbers  in  air  per  foot  of  model  chord  length  is 
given  in  Figure  3  as  a  function  of  flow  Mach  number.  Note  that  Reynolds  num¬ 
bers  in  excess  of  3  x  10®  are  available  for  flow  Mach  numbers  ranging  from  0.05 
to  1.2.  If  helium,  for  example,  were  used  in  the  compression  chamber,  even 
larger  Reynolds  mi«ber3  could  be  realized  with  air  as  the  working  medium  in  the 
hot  gas  region;  these  are  also  shown  in  Figure  3« 

Figures  U,  5,  and  6  show  respectively  an  overall  view  of  the  shock  tube, 
the  diaphragm  section,  and  the  test  section  together  with  the  interferometer 
and  the  shock  detecting  stations.  Incidentally,  the  nominal  flow  Mach  number 
produced  in  the  shock  tube  is  obtained  by  a  conventional  measurement  of  the 
shock  velocity  (Just  ahead  of  the  test  section)  and  measurement  of  temperature 
of  the  expansion-chamber  air. 

In  the  following,  a  brief  outline  of  some  of  the  problems,  in  each  of 
the  originally  mentioned  three  areas,  which  have  been  or  are  planned  to  be 
investigated  are  given.  Results  of  interferometric  measurements  of  transient 
airflow  development  about  a  double -wedge  airfoil  in  M  ■  0.4*  hot  gas  flow 
will  be  reviewed  in  somewhat  mere  detail.  This  will  serve  to  point  up  some 
of  the  problem  areas  in  the  measurement  techniques  realm  as  well  as  problems 
posed  by  the  basic  shock  tube  flow  itself. 

In  general,  the  response  of  aircraft  structures  to  transient  air  loads 
can  be  predicted  satisfactorily  without  the  inclusion  of  gravitational  forces. 
Only  elastic  and  inertial  simulation  of  an  aircraft  structure  Is  required. 
Generally  speaking,  scaling  of  a  prototype  aircraft  structure  cannot  be 
achieved'-  reaii- tlcally  when  the  scaling  reduction  factor  exceeds  roughly  *». 


*  In  this  paper,  Mach  number  M  refers  to  the  (particle)  Mach  number  of  the 
hot  gas  flow  region.  Shock  Mach  number  is  never  used. 


Hovever,  a  mathematical  model  describing  the  dynamic  behavior  of  prototype  air¬ 
craft  structures  can  generally  be  devised.  The  dynamic  characteristics  of  this 
mathematical  model  are  often  simulated  by  much  smaller  than  one  -quarter  scale 
models)  at  the  same  time,  the  Internal  structure  may  be  quite  dissimilar  to 
that  of  the  prototype.  If,  on  the  other  hand,  one  wishes  to  Investigate  failure 
modes  of  a  prototype  structure  under  dynamic  loading  conditions,  the  model  si¬ 
mulation  must  be  carried  out  to  a  much  more  complete  degree. 

Figure  7  Illustrates  a  simple  cantilever-mounted  rectangular  airfoil 
model  consisting  of  an  aluminum  sheet-spar  covered  vith  balsa  wood  which  haa 
been  contoured  to  a  typical  airfoil  shape.  This  model  was  set  at  a  small  angle 
of  attack  and  was  then  subjected  to  a  shock  wave  and  its  afterflow.  The  re¬ 
sponse  of  this  elastic  model  was  measured  by  means  of  Btraln  gages  mounted  on 
opposite  sides  of  the  spar  near  the  root)  a  typical  strain-time  record  is  shown 
In  Figure  8.  Upon  being  struck  by  the  shock,  the  model  begins  to  vibrate  and 
tends  to  seek  an  equilibrium  strained  position  corresponding  to  tbs  steady 
state  air  loads  produced  on  the  model  In  the  "uniform"  hot  gas  flow  flsld  be¬ 
hind  the  shock.  For  this  case,  the  shock  reflected  from  the  closed  end  of  the 
expansion  chamber  arrives  at  the  test  section  at  about  0.043  second  after  the 
arrival  of  the  incident  shock  and  stops  the  flow.  The  model  then  tends  to 
seek  Its  initial  "unstrained"  position.  Subsequently,  the  rodel  is  observed 
to  respond  to  a  whole  succession  of  waves  which,  for  some  time,  continue  to 
reflect  from  the  closed  ends  of  the  shock  tuba. 

Similar  studies  have  been  conducted  wherein  a  smaller  cantilever  model, 
similarly  instrumented,  has  been  placed  at  a  number  of  positions  downstream  as 
well  as  above  and  below  the  streamwlse  plane  of  the  model  shown  In  Figure  7* 
Tests  were  run  to  study  the  effect  of  the  "starting  vortex”  (and  the  subsequent 
dovnwash  field)  from  the  upstream  airfoil  on  the  response  of  the  downstream 
airfoil  (tall).  For  each  position  of  the  tall  model,  tests  were  run  at  the 
same  initial  flow  conditions  both  in  the  presence  and  In  the  absence  of  the  up¬ 
stream  cantilever-wing  model.  The  response  of  the  tail  vas  observed  to  vary 
quite  markedly,  depending  upon  its  location  both  downstream  and  normal  to  the 
wing  stream  plane.  In  Reference  6  the  characteristics  of  such  downwash  for  tha 
two-dimensional  incompressible  flow  case  are  treated,  and  the  response  of  sim¬ 
ple  systems  to  the  associated  airloading  is  discussed. 

Another  example  of  aeroelastlc  and  structural  response  studies  con¬ 
ducted  In  the  shock  tube  follows.  The  question  posed  was  that  of  the  time 
required  to  produce  flutter  of  a  wing  should  the  wing  find  Itself  abruptly 
In  a  field  of  flow  in  which  the  material  velocity  were  In  excess  of  that  ra- 
qulred  to  produce  flutter  in  steady  state  flow'.  A  cantilever  flutter  model 
consisting  of  an  aluminum  B-spar,  covered  with  balsa  wood,  and  airfoil  con¬ 
toured,  was  designed  and  built  (6ee  Figure  9)  •  1°  order  to  control  the  tor¬ 

sional  stiffness  of  the  model,  the  balsa  wood  contouring  was  slit  chordwlse 
at  s.  number  of  spanwlse  stations.  The  contour  was  preserved  and  the  balsa 
covering  protected  (from  diaphragm  particles  and  other  foreign  matter)  by  a 
thin  rubber  sheet  applied  to  the  model.  For  this  model,  the  ratio  of  the 
first  bending  to  the  first  torsion-mode  frequency  was  approximately  unity. 

The  model  was  strain-gage  instrumented  both  (a)  to  reveal  In  vibration  teata, 
the  bending  and  torsional  frequencies  of  the  model  and  (b)  to  permit  mea¬ 
suring  model  response  In  the  actual  shock  tube  experiments.  The  model  was 
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cantilever  mounted  and  was  subjected  to  shock-initiated  flow  starting  at  small 
material  velocities.  The  material  velocity  of  the  shock-initiated  flow  was 
increased  by  small  increments  in  succeeding  runs.  During  each  run,  the  (strain) 
response  of  the  model  was  oscillographically  recorded.  At  small  material  velo¬ 
cities,  the  vibratory  response  decayed  essentially  monotonically  until  the  ini¬ 
tial  hot  gas  flow  wa3  terminated  by  the  reflected  shock.  On  the  other  hand, 
when  the  shock-initiated  material  velocity  exceeded  the  steady  state  flutter 
speed  c  ?  the  model,1  the  vibratory  response  was  observed  to  decay  for  a  short 
period  of  time  and  then  diverged.  The  early  arrival  of  the  shock  reflected 
from  the  closed  end  of  the  expansion  chamber  stopped  the  flow  and  prevented 
catastrophic  failure  of  the  flutter  model.  Typical  response  records  are  given 
in  Figure  10  for  material  velocities  both  well  below  and  well  above  the  flutter 
speed.  The  time  at  which  the  response  began  to  diverge  noticeably  In  the 
Figure  10  record  corresponds  to  approximately  22  flow  chord  lengths;  steady 
state  air  forces  are  observed  tc  be  developed,  for  such  flow  conditions,  in 
roughly  10  flow  chord  lengths.  The  aerodynamic  time  lag  observed  may  be  en¬ 
tirely  different  for  other  combinations  of  steady  state  flutter  flow  Mach  num¬ 
ber,  bending  frequency,  torsional  frequency,  mass  ratio,  and  so  forth. 

An  additional  aeroelastic  response  study  conducted  consisted  of  sub¬ 
jecting  thin-walled  cylindrical  shell  models  to  the  shock  wave  and  measuring 
the  model's  strain-time  response  at  a  number  of  points  along  the  periphery”. 

The  shock  wave  overpressures  employed  were  such  that  the  observed  response  re¬ 
mained  in  the  (linear)  elastic  region.  For  these  models  the  response  modes 
are  of  the  "dilatory"  variety,  whereas,  the  modes  of.  most  importance  in  the 
case  of  the  wing-tail  cantilever  models  and  the  aforementioned  flutter  models 
were  those  of  bending  and  torsion;  the  responses  again  in  these  latter  cases 
ware  confined  to  the  elastic  region.  Sections  of  typical  aircraft  lifting- 
surface  structures  were  suitably  mounted  between  the  shock  tube  walls  and 
placed  at  various  orientations  with  respect  to  the  incident  shock.  For  the 
case  in  which  the  shock  struck  the  nose  of  the  airfoil  and  enveloped  the 
structure  symmetrically,  skin-panel  vibrations  and  vibrations  of  the  rib 
structure  were  produced;  these  modes  of  vibration  may  be  thought  of  as 
"dilatory"  modes.  Such  modes  were  excited  to  the  extent  that  permanent  damage 
resulted  in  some  cases  and  only  elastic  response  in  others.  Figure  11a  shows 
a  (nose-on)  symmetrical  exposure  condition  wherein  dilatory  nodes  were  domi¬ 
nant;  rib-crushing  was  a  typical  type  of  damage  produced,  and  skin-dishing  was 
another.  Such  damage  is  termed  "overpressure  damage".  In  Figure  lib,  a  typical 
model  arrangement  is  illustrated  for  which  "beam-bending"  modes  were  dominant. 
For  the  latter  case,  first,  elastic-response  tests  were  conducted.  These  were 
followed  by  tests  in  which  the  response  produced  extended  far  beyond  the  elas¬ 
tic  region.  The  incipient  "failure"  of  the  structure  in  the  beam-bending  mode 
of  response  occurred  through  local  buckling  of  the  structure.  After  such  a 
buckle  occurs,  the  linear  elastic-restoring  force  characteristics  of  the 
model  no  longer  exist;  as  the  model  continues  to  deflect,  the  rertoring 
force  characteristics  continually  change  in  an  irreversible  non-linear  manner. 
Damage  resulting  from  such  beam-bending  modes  of  vibration  is  termed  "gust 
damage" . 


An  analysis  to  predict  the  elastic  and  tne  post-elastic  responses  of 
such  structural  models  must  properly  include  the  pertinent  restoring  force 
characteristic  of  the  structure  a3  a  function  of  the  controlling  deflection 


paraseter(s) j  the  air  loads  causing  the  response  must  also  be  properly  taken 
Into  account.  The  latter  are  unquestionably  the  more  difficult  to  treat* 

A  primary  purpose  in  developing  this  shock  tube  facility  was  to  inves¬ 
tigate  air  loads  produced  on  airfoils  and  other  aerodynamic  bodies  by  a  shock 
wave  and  the  associated  flow  field.  Two-dimensional  airfoils  (an  KACA  6J-010 
and  a  10%  thick  symmetrical  double-wedge)  vere  mounted  between  the  test  section 
windows  as  illustrated  in  Figure  12.  The  airfoil  was  set  at  some  pre-selected 
angle  of  attack  relative  to  the  floor  of  the  shock  tube,  and  vas  then  subjected 
to  shock-initiated  flow.  Shock- initiated  flow  Mach  numbers  of  0.40,  1.28,  and 
1.5  were  employed  in  tests  on  the  double-wedge  airfoil  which  was  set  at  angles 
of  attack  of  0,  4,  and  8.3  degrees.  The  KACA  65-OIO  airfoil  has  been  tested 
at  M  ■  0.4  and  angles  of  attack  of  0,  5,  10,  15,  and  30  degrees*  At  least 
two  Reynolds  numbers  vere  employed  in  testing  each  airfoil.  During  each  firing 
of  the  shock  tube,  in  these  tests,  a  single  lnterferogran  of  the  flow  about  the 
airfoil  vas  made  at  a  pre-selected  time  after  the  flow-initiating  shock  im¬ 
pinged  upon  the  leading  edg  of  the  model,  using  a  spark  light  source  of 
1-micrcseccnd  duration  and  filtered  to  a  30*A  band  centered  at  448l*A.  By 
repeating  such  runs  and  varying  the  time  at  which  the  lnterferogran  vas  re¬ 
corded,  a  reasonably  complete  measure  of  the  flow  pattern  development  fron 
shock  incidence  to  steady  state  flow  conditions  vas  obtained* 

To  illustrate  some  of  the  measurements  obtained  and  some  of  the  inter¬ 
pretation  problems  involved  In  such  measurements,  a  set  of  runs  for  the  double¬ 
wedge  airfoil  at  M  -  0.4  will  be  discussed.  Figure  13  illustrates  a  series  of 
infinite  fringe  flow  interf erograms  taken  at  vaious  chord  lengths  of  particle 
flow,  3/c,  after  the  shock  struck  the  leading  edge  of  the  model.  For  this  set 
of  runs,  the  angle  of  attack  is  4  degrees  and  the  flew  Reynolds  nuaber  based 
on  the  model  chord  is  1,210,000.  This  and  other  sets  of  runs  for  tbe  double- 
wedge  airfoil  are  discussed  in  some  detail  in  Reference  9* 

Having  information  available  on  the  effective  wave  length  of  the  light 
used,  the  test  section  width,  and  the  pertinent  value  of  the  0 lads tone -Dale 
constant,  the  difference  in  mean  density  from  one  fringe  to  the  next  in  the 
infinite  fringe  flow  interferogram  is  a  known  value.  Thus  if  the  fid.d  den¬ 
sity  at  any  point  in  the  flow  is  known,  the  density  at  all  fringe  positions 
throughout  the  flow  field  are  readily  determined.  Having  the  density  distri¬ 
bution,  one  proceeds  to  determine  the  surface  pressure  distribution  on  the 
airfoil.  Except  where  dir, continuities  are  present  in  the  flow  field  It  la 
reasonable,  under  these  test  conditions,  to  regard  the  flow  as  lsentroplc  In 
translating  density  to  pressure.  In  assessing  the  pressure  distribution  on 
the  airfoil,  one  observes  that  the  boundary  layer  in  some  regions  has  n  re¬ 
spectable  thickness.  It  16  assumed  that  the  pressure  gradient  across  the 
boundary  layer  is  negligible;  hence,  the  surface  pressure  is  the  same  as  that 
at  the  outer  edge  of  the  boundary  layer.  Kote  that  on  the  upper  surface  be¬ 
hind  the  sharp  leading  edge  of  the  airfoil,  there  is  a  region  of  separated  flour 
in  the  near  steady  state  flow  condition.  This  is  a  region  in  which  a  reliable 
determination  of  surface  pressure  distribution  is  most  difficult.  Incidentally, 
in  order  to  provide  independent  information  on  the  complete  pressure  distribu¬ 
tion  and  to  guide  interferogram  interpretation  in  these  separated  flow  regions, 
direct  pressure  distribution  measurements  are  planned  using  a  miniature  pressure 
sensor  which  is  currently  being  developed  in  the  laboratory.  A  discussion  of 
interferogram  interpretation  for  those  cases  Involving  discontinuities  In  the 
flow  field  may  be  found  in  Reference  9* 
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Once  the  surface  pressure  distribution  Is  known,  one  can  evaluate  the  lift,  the 
moment,  and  the  pressure  drag  (and  other  Interesting  quantities)  experienced  by 
the  airfoil.  This  information  can  converted  to  the  more  convenient  coeffi¬ 
cient  form  if  the  freestream  flow  properties  are  known.  Ideally,  ^he  hot  gas 
flow  is  a  step  function.  Actually,  non-ideal  flow  behavior  is  observed;  non¬ 
ideal  diaphragm  removal  and* boundary  layer  growth  are  some  of  the  factors 
responsible.  Measuring  white  light  finite  (width)  fringe  displacement  time 
histories  in  the  flow  by  means  of  a  strip  film  camera  and  by  means  of  a  series 
of  single-spark  white  light  interferograms,  time  histories  of  the  free  stream 
flow  density  were  determined.  If  one  postulates  that  at  some  (late)  time 
after  shock  passage,  for  example,  the  density  increase  over  that  observed  im¬ 
mediately  behind  the  shock  were  due  to  one  shock,  two  shocks,  or  a  large  num¬ 
ber  of  weak  compression  waves,  the  re*  llting  predicted  material  velocity  is 
found  to  vary  but  little  from  one  assumed  model  to  any  other.  Hence,  if  such 
a  flow  model  is  postulated,  one  may,  throughout  the  flow,  calculate  the  mater¬ 
ial  velocity  and  thus  the  (tentative)  time  history  of  the  dynamic  pressure  of 
the  freestream  flow.  Figure  14  illustrates  the  resulting  time  history  of  tha 
free  stream  dynamic  pressure  for  the  conditions  of  this  experiment;  the  shock- 
velocity- indicated  material  Mach  number  immediately  behind  the  shock  at'  the 
measuring  station  was  0.4.  Figure  15  shows  the  variation  of  lift  coefficient 
with  elapsed  flow  chord  lengths  (s/c)  when  the  lift  coefficient  is  based  upon 
the  dynamic  pressure  immediately  behind  the  shock.  Figure  16  illustrates  the 
aerodynamic  time  variation  of  the  lift  coefficient  based  upon  local  (time) 
values  of  the  density-indicated  dynamic  pressure.  In  both  figures,  curve*  of 
lift  response  to  a  step-function  imcompresBible  flow®  are  given;  these  provide 
only  rough  theoretical-experimental  comparisons  since  the  theoretical  flow 
model  is  Bomevhat  different  from  the  conditions  applying  in  the  6hock  tube  ex¬ 
periments.  Figure  17  6howa  similar  results  for  the  (early)  60-called  "shock 
diffraction  phase"  of  the  loading. 

Figure  18  illustrates  a  series  of  infinite  fringe  Interferograms  for 
the  double-wedge  airfoil  at  a  nominal  shock- initiated  flow  Mach  number  of 
1.28  and  an  angle  of. attack  of  4  degrees.  The  shock  should  be  detached  at  this 
Mach  number.  The  last  picture  in  this  series  was  taken  in  the  (much  higher 
Mach  number)  cold  gas  flow  following  the  contact  surface.  Note  that  this  cold 
gas  flow,  in  contrast  to  the  hot  gas  flow,  is  quite  non-uniform.  It  is  inter¬ 
esting  also  to  observe  the  oscillatory  wake  which  appears  in  the  latter  stages 
of  the  hot  gas  flow  pictures. 

In  addition  to  the  tvo-dimensional  airfoil  experiments  mentioned,  plans 
have  been  made  for  conducting  similar  experiments  on  finite-span  wings  using 
the  aforementioned  miniature  pressure  sensorB  for  pressure  distribution  mea¬ 
surements.  These  measurements  will  be  supplemented  by  high  speed  photographic 
studies  of  the  behavior  of  tufts  distributed  over  the  surface  of  the  wing. 
Incidentally,  such  tuft  studies  have  been  made  on  the  two-dimensional 
HACA  65-OIO  airfoil  using  a  Fastax  camera  operated  at  16,000  frames  per  second. 
The  tuft  behavior  has  been  correlated  with  measurements  made  at  corresponding 
times  from  interferograms. 

Mention  should  be  made  of  the  fact  that  in  such  shock  tube  studies, 
the  question  of  (classical,  non-viscous)  wall  interference  in  addition  to 
wall-boundary-layer  interference  must  be  taken  into  consideration.  Such 
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interference  is  not  serious  in  the  very  early  stagee  of  the  flow;  however, 
multiple  wave  reflections  and  wall-boundary- layer  interference  sire  clearly 
evident  later  on  in  the  flow.  Although  steady  state  wall  interference  effects 
have  been  treated  extensively  in  the  literature,  the  transient  condition  of 
Interest  here  has  net  been  studied. 

Other  simple  models  have  been  subjected  to  shock-initiated  flow  and 
the  resulting  flow  patterns  observed.  These  have  included  bodies  of  square, 
rectangular,  and  circular  cross  section.  Figure  19  shows  a  sequence  of  in¬ 
finite  fringe  interferograms  of  the  flow  about  a  l/4  by  2  inch  rectangular 
model  for  a  Mach  number  of  0.17.  Tn  Figure  20,  Schlieren  picture  and  an 
infinite-fringe  flow  lnterferogram  taken  in  separate  runs  but  at  nearly  the  same 
aerodynamic  time  in  the  flow  are  shown;  here  the  flow  Mach  number  is  0.45 . 

A  5/l6- inch  diameter  model  was  subjected  to  M  ■*  0.17  flow  and  interferograms 
obtained;  a  sequence  of  these  pictures  is  shown  in  Figure  21.  Note  that 
symmetrical  vortices  are  first  shed;  later  in  the  flow,  the  Karman-type  vortex 
street  develops. 

Several  tests  were  also  conducted  on  an  axially  symmetric  body  con¬ 
sisting  of  a  simple  cone-type  nose  followed  by  a  cylindrical  afterbody^;  a 
half  model  was  attached  .0  the  side  wall  of  the  shock  tube.  Surface  pressures 
were  sensed  at  several  streamwise  positions  along  the  body  by  means  of  barium 
titanite  pressure  transducers.  The  pressure-time  histories  were  recorded  by 
photographing  the  sweeping  trace  of  a  cathode  ray  oscilloscope.  Similar  mea¬ 
surements  have  been  made  of  pressure -distribution  time  histories  on  two  low- 
aspect-ratio  finite-span  rectangular  wings  cantilever-mounted  from  the  shock 
tube  wall^j  however,  the  pressure  leads  from  the  model  surface  to  the  barium 
titanite  pressure  transducers  were  of  such  a  length  and  diameter  that' the 
frequency  response  of  the  sensing  system  was  reduced  to  the  point  where  the 
diffraction  loading  could  not  he  observed.  Only  the  latter  stages  of  the  flow 
development  and  the  steady  state  air  forces  were  successfully  measured.  This 
measurement  limitation  will  he  removed  by  the  miniature  pressure  sensors  now 
being  developed  in  the  laboratory. 

In  support  of  theoretical  work  by  Hobbs  in  Reference  6,  experiments 
were  conducted  in  which  a  two-dimensional  airfoil  was  placed  Just  upstream  of 
the  test  section  and  set  at  a  small  angle  of  attack.  A  schematic  of  this 
arrangement  is  shown  in  Figure  22,  part  a.  The  path  of  the  "starting  vortex" 
was  observed  both  in  space  and  In  time  by  photographing  the  vortex  at  succes¬ 
sive  time  delays  in  a  series  of  repetitive  tests.  Excellent  repetition  was 
observed.  One  prime  end  result  was  that  of  measuring  the  downwash  induced 
by  the  vortex  at  various  points  in  its  vicinity.  Following  a  suggestion  by 
Dr.  Ruetenlk,  a  small  diameter  wire  was  strung  between  the  test-section  win¬ 
dows.  Photographs  were  taken  defining  the  wake  behind  this  fine  wire  when  the 
subject  vortex  was  at  various  positions  with  respect  to  the  wire;  the  inclina¬ 
tion  of  this  wake  served  to  indicate  the  "downwash  velocity"  component  pro¬ 
duced  at  the  sensing  wire,  since  the  freestream  flow  velocity  was  reasonably 
well  known.  While  this  technique  has  yielded  useful  results,  the  Inherent 
uncertainties  indicate  the  need  for  a  more  precise  measurement  method.  A 
satisfactory  prediction  of  this  "lndicial"  downwash  field  permits,  in  turn, 
with  the  aid  of  available  theory,  the  prediction  of  air  forces  on  airfoils 
which  come  under  the  influence  of  such  an  "indicial  vortlcity  field". 


Incidentally,  in  the  earller-deBcribed  wingtail  els. a tic  response  tests,  the 
influence  on  tail  response  of  the  "starting  vortex"  frc-m  the  ving  vas  observed 
while  that  vortex  passed  near  the  tail.  Part  b  of  Figure  22  shows  the  space- 
time  history  of  the  vortex  path.  The  vortex  does  not  imsediately  upon  its 
initial  appearance  acquire  a  downstream  velocity  equal  to  the  material  velocity 
itself.  It  starts  its  Journey  slowly  anl  soon  attains  essentially  a  velocity 
equal  to  the  material  velocity.  For  the  case  illustrated  in  Figure  22,  the 
vortex  arrived  at  the  test  section  center  line  late  by  approximately  one-eigth 
of  a  chord  length  based  upon  a  constant  value  of  shock-velocity-indicated 
material,  velocity. 

Both  on  its  own  merits  and  for  a  meaningful  reduction  of  the  previously 
described  transient  air  force  measurements,  it  is  Important  that  the  basic  free- 
Btream  flow  in  the  shock  tube  be  known  accurately.  Over  the  complete  range  of 
available  and/or  practical  operating  conditions,  one  is  Interested  in  a  com¬ 
plete  time  history  oft 
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Measurements  have  been  made  of  items  (a)  anu  (b)  in  many  shock  tubes^'  . 

Time  histories  of  these  quantities  have  been  very  strongly  linked  to  the  di¬ 
mensions  of  the  tube  cross  sections  and  lengths  employed  through  the  associated 
boundary  layer  development,  non-ideal  diaphragm  bursting,  and  other  effects. 

All  of  the  items  listed  above  are  interdependent.  Numerous  papers5#  13*16 
have  been  written  on  non-ideal  shock  tube  flow,  shock  attenuation,  and  boundary 
layer  development;  review  of  these  papers  will  not  be  attempted  here. 

8ome  exploratory  measurements  have  been  made  at  M.I.T.  by  Dr.  Ruetenik 
seeking  to  measure  stream  temperature,  material  velocity,  and  mas6  flow  with  a 
hot  wire  anemometer;  he  vlll  discuss  this  topic  in  hia  paper. 

The  quet tion  of  the  turbulence  level  of  the  shock-initiated  hot  gas 
flow  is  one  that  appears  as  of  now  to  be  unanswered.  The  cold  gas  flow,  cm 
the  other  hand,  has  been  observed  to  be  quite  irregular  — -  irregular  to  the 
point  where  the  value  of  using  this  flow  for  air  force  studies  on  models  is 
open  to  question.  As  to  the  hot  gas  flow,  this  laboratory  has  proposed  to 
conduct  a  series  of  tests  designed  to  examine  the  turbulence  level  present. 
Also  to  be  explored  is  the  degree  of  repeatability  of  such  turbulence  levels 
and  time  histories  in  "repetitive"  runs. 

Flow  "cleanliness",  in  the  sense  used  here,  pertains  to  the  number 
and  strengtha  of  trailing  shocks  and  rarefaction  waves  appearing  in  the  hot 
gas  flow  which  is  considered  to  be  the  flow  region  of  primary  Interest  for 
transient  aerodynamic  studies  on  models.  The  trailing  wave  family  appears  to 
depend  upon  the  diaphragm  support  configuration  employed,  the  pre-burst  cur¬ 
vature  of  these  diaphragms,  diaphragm  rupture  rate,  amount  of  diaphragm  re¬ 
maining,  and  distance  of  the  diaphragm  from  the  test  station,  as  well  as  any 
wall  protuberances  or  irregularities.  Such  effects  have  been  noted  to  a 
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greater  or  lesser  degree  In  all  shock  tube  facilities.  Everyone  agrees  that 
every  feasible  step  should  be  taken  to  minimize  these  effects. 

Wall-boundary  layer  studies  have  been  undertaken  here*®*  *®*  *^  and  at 
Princeton^  .  Loth  wall  temperature  and  Interferometric  measurements*’  have 
been  reported  at  Princeton.  Only  interferometric  measurements  have  thus  far 
been  made  at  M.I.T.  Some  of  these  results  will  be  discussed  in  Dr.  Ruetenlk's 
paper20. 


Aside  from  the  ideal  step-function  flow  previously  considered,  there 
are  certain  prob  lems  where  in  flow  time  histories  other  than  a  step-function  are 
of  interest.  Blast-type  profiles  are  of  particular  interest.  Such  profiles 
can  be  produced  in  the  Bhock  tube  by  explosions2!  (gaseous  mixtures,  solids, 
etc.),  by  arranging  rarefaction  overtaking  of  the  primary  shock22,  by  the  use 
of  an  auxiliary  diaphragm  placed  upstream  of  the  test  statlon23,  and  by  the 
shock  and  flow  expansion  through  a  slit2’,  to  mention  a  few.  Such  profiles  do 
not  exactly  duplicate  those  observed  in  free-air  three-dimensional  explosions; 
the  greater  part  of  the  positive  phase  Is  similar  but  differences  arise  in  the 
latter  stages  of  these  shock  tube  blast-like  profiles. 

In  order  to  better  simulate  the  condition  in  which  a  wing  in  steady 
(high-speed)  flight  is  struck  by  a  blast  vave,  work  is  In  progress  to  devise 
appropriate  flows  in  the  shock  tube.  A  number  of  techniques  have  been  consi¬ 
dered.  One  technique  which  shews  promise  has  been  explored  here.  It  is  shown 
schematically  In  Figure  23}  the  geometry  must  be  tailored  to  the  primary  oper¬ 
ating  conditions  selected.  This  technique  consists  of  producing  shock-initia¬ 
ted  .flow  In  the  usual  manner;  near  the  test  station,  a  part  of  the  initial 
shock  is  split-off  and  entrapped.  This  entrapped  9hock  is  then  reflected  at 
a  time  such  that  when  tbe  resulting  wave  expands  into  the  primary  flow,  this 
auxiliary  blast-like  wave  arrives  at  the  model  being  tested  at  a  time  when 
essentially  steady  state  air  forces  have  already  developed  on  the  model. 

An  attempt  has  been  made  herein  to  review  some  of  the  shock  tube  ac¬ 
tivities  carried  out  at  M.I.T.  A  number  of  the  measurement  problems  encoun¬ 
tered  at  M.I.T.  are  common  to  all  shock  tube  facilities.  By  reviewing  the 
deficiencies  in  (and  of)  measurement  techniques  particularly,  as  well  as 
Ideas  on  unresolved  shock  tube  f~ow  characteristics,  vich  those  who  face  com- 
non  and  allied  problems,  it  is  hoped  that  a  mutual  interchange  of  Ideas  and 
experience  will  expedite  studies  which  face  each  of  us. 

If  we  seek  ultimately  to  predict  the  responses  of  various  structures 
to  transient  air  forces,  we  must  have  an  accurate  knowledge  of  the  applied 
air  loads.  It  is  clear  that  such  air  force  information  cannot  be  determined 
to  the  precision  required  or  in  a  form  reeded  for  practical  use  by  attempting 
to  infer  such  airloads  information  from  structural  response  data.  Finally, 
aside  from  questions  of  basic  shock  tube  flow  characteristics  themselves,  the 
usefulness  of  the  shock  tube  as  a  tool  for  providing  air  force  information 
and/or  for  developing  appropriate  methoda  for  predicting  structural  response 
is  critically  dependent  upon  tbe  availability  of  rapidly  responding  precision 
instrumentation . 
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COMPRESSION 

CHAMBER 


EXPANSION  CHAMBER 


INTERNAL  CROSS-SECTION  8“  X 24" 

INTERNAL  LENGTH  981-!  j 

APPROXIMATE  TOTAL  WEIGHT  OF  SH0CK-TU8E 

PROPER  "'V'IS  Ton* 

MAXIMUM  DESIGN  PRESSURE  250  ptlo 

MAXIMUM  REYNOLDS  NO.  WITH  A  4- IN.  CHORD 

MODEL  WITH  AIR-AIR  OPERATION  6,750,000 

MAXIMUM  MACH  NO.  CURRENTLY  OBTAINABLE 

WITH  AIR  IN  THE  EXPANSION  CHAMBER  1.5 


Figure  II  Summary  of  Physical  Data  on  the 
M.I.T.-W.A.D.C.  Shock  Tube 
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MATERIAL  MACH  NUMBER 


Figure  k.  Overall  Flew  of  the  Shock  Tube 
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Figure  7*  Rectangular  Elastic  Ving  Cantilevered 
from  the  Shock  Tube  Vail 


Figure  8.  Response  of  Elastic  Cantilevered  Wing  to 
Shock-Initiated  Flow 
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(o)  CONSTRUCTION  FEATURES 


u  =  184  MPH 

(o)  AT  LESS  THAN  FLUTTER  SPEED 


I 


u  =  298  MPH 

(b)  BEYOND  FLUTTER  SPEED 


Figure  10.  typical  Strain-Time  Responses  of  the 

Flutter  Model 
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Figure  12.  Double-Wedge  Airfoil  Mounted  between  Shock  Tube  Windows 


U  =  SHOCK  VELOCITY 

u  =  MATERIAL  VELOCITY 

t  =  MEASURED  FROM  SHOCK  IMPING EMENT  AT  L.E. 

»/c  »  (uO/c  =  FLOW  CHORD  LENGTHS  ELAPSED  FOR 

ANY  GIVEN  PICTURE 


NOMENCLATURE 


s/c  3  0.03 


Figure  13.  Sequence  of  Infinite  Fringe  Flow 
Interferograas  for  the  Double- 
Wedge  Airfoilj  M  -  0.4,  am  k9 


Figure  13*  Continued 
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HOBBS 

EXPERIMENTAL 


FLOW  CHORD  LENGTHS  .  s/C 


s/c  *  0.08 


s/c  »  0.83 


Figure  lb.  Flow  Interferogram  Sequence;  Double-Wedge 
Airfoil;  M  «  1.28,  ac  «  4* 
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Figure  19.  Flow  Interferograra  Sequence;  Flat 
Plate  Model;  M  *  0.17 
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(a)  SCHEMATIC  LAYOUT 
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(b)  VORTEX  MOTION  DATA 


Figure  22.  Schematic  for  Vertex  Motion  and  Vortex-Induced 

Dovnvask  Study 
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EARLY  PATTERN 


STUDY  OF  GRIDS  IN  SHOCK  TUBES 
A.  Wiederaann 

Ano'Jr  Research  Foundation 


The  work  reported  In  thir  talk  was  sponsored  by  the  Air  Force  Special 
Weapons  Center  in  connection  with  the  operation  of  the  Air  Force  Shock  Tube 
Laboratory.  The  objective  was  to  determine  quantitatively  the  effectiveness 
of  grids  In  shock  tubes  to  eliminate,  or  reduce,  the  reflected  signals  from  the 
end  of  the  shock  tube.  We  considered  two  specific  problems!  (l)  placing  the 
grid  at  the  end  of  a  constant  area  channel,  and  (2)  placing  the  grid  inside  a 
constant  area  channel.  I  will  present  the  results  of  both, but  because  of  the 
time  limitation  I  will  omit  the  details  of  the  second  case  which  are  quit® 
similar  to  those  of  the  first. 

It  is  well  known  that  when  a  flat-topped  blast  wave  (that  Is,  &  shock 
wave  followed  by  a  uniform  flow  field)  is  incident  upon  the  closed  end  of  a 
shock  tube,  the  shock  wave  Is  reflected  and  moves  back  upstream  leaving  a  uni¬ 
form  flow  field  downstream.  Also,  if  the  end  of  the  tube  is  open,  the  ensuing 
interaction  results  in  the  propagation  of  a  rarefaction  wave  back  upstream, 
such  that  either  the  exit  pressure  is  equal  to  the  original  ambient  pressure, 
or  the  flow  at  the  exit  is  choked  (that  Is,  the  exit  Mach  number  Is  unity). 

The  pressure  variations  chat  would  be  observed  a  short  distance  upstream  from, 
the  end  of  the  tube  are  illustrated  in  Figure  1  for  these  cases. 

For  any  grid  with  an  area  ratio,  say,  sigma  (•*),  which  Is  defined  as 
the  ratio  of  open  area  to  total  area,  between  the  limits  zero  (the  closed  case) 
and  unity  (the  open  case),  a  weaker  shock  wave  or  rarefaction  wave  will  be  pro¬ 
pagated  back  upstream.  It  appears  reasonable  to  assume  that  there  exists  a 
critical  area  ratio,  say,  sigma  bar  (?)  fer  which  no  signal  will  be  propagated 
back  upstream. 

Observations  of  the  interaction  of  shock  waves  and  grids  in  constant 
area  channels  indicate  that  the  resulting  flow  can  be  approximated  by  steady 
flow  tlvrough  the  grid  (or  nezeies)  ana  quasi -steady  flow  throughout  the  re¬ 
mainder  of  the  flow  field.  In  other  words,  the  diffraction  signals  coalesce 
after  they  have  propagated  several  nozzle  lengths  (some  characteristic  dimen¬ 
sion  of  the  nozzle).  The  radial  gradients  are  s'nall  compared  to  the  axial 
gradients  (except  in  the  nozzle),  hence,  a  one-dimensional  analysis  is  appli¬ 
cable.  Figure  2  illustrates  the  flow  fields  which  result.  Regions  1  through 
4  are  relatively  free  of  disturbances  and  are  assumed  to  be  uniform  in  the 
analysis.  For  the  case  where  the  grid  is  placed  at  the  end  of  a  constant  area 
channel,  the  region  to  the  right  of  the  "grid  position"  does  not  exist,  and 
either  a  shock  wave  is  reflected  back  upstream--as  shown  here--or  a  rarefac¬ 
tion  is  reflected  back  upstream. 
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The  nozzles  placed  at  the  end  of  the  tube  can  be  either  converging- 
diverging  or  converging  with  a  sudden  expansion.  The  boundary  condition*  la- 
posed  by  the  grid  aret  (l)  The  exit  pressure  is  equal  to  the  original  ambient 
pressure,  or  (2)  the  flow  In  the  throat  of  the  nozzle  Is  choked.  The  exit 
pressure  state  is  In  the  throat  for  the  converging  nozzle  and  at  the  end  of  the 
nozzle  for  the  ccnverglng*-dlverging  nozzle.  The  equations  used  aret  (l)  the 
Rankine -Hugonlot  equations  for  the  states  separated  by  shock  waves,  (2)  the 
lsentroplc  non- steady  flow  equations  for  states  separated  by  a  centered  rare¬ 
faction  wave,  and  (3)  the  steady  one -dimensional  lsentroplc  flow  equations 
betveen  the  entrance,  throat,  and  exit  states  in  the  nozzle.  I  will  not  pre¬ 
sent  these  equations  here  since  moBt  of  you  are  quite  familiar  with  them.  If 
ve  go  back  to  Figure  1,  ve  can  define  a  quantity  which  will  measure  the  degree 
of  the  disturbance.  This  quantity  will  be  called  £  and  is  defined  as 


If  t,  is  the  shock  strength  or  the  Incident  shock  wave  (that  is,  the  ratio  of 
the  absolute  pressures  across  the  shock  wave),  4'  18  the  strength  of  the 
reflected  shock  wave,  and  Tj  Is  the  strength  of  the  rarefaction  wave,  then 


for  the  reflected  shock  wave 


and 


«?- 1) 


for  the  reflected  rarefaction  wave. 


A  numerical  solution  of  the  governing  system  of  equations  results  in 
the  plots  of  Figure  3,  which  was  made  for  a  gas  with  a  ratio  of  specific  heats 
of  l.H.  All  the  numerical  results  presented  here  are  based  upon  this  value. 

The  results  for  both  types  of  nozzles  are  presented.  The  dashed  line 
Is  the  locus  of  solutions  where  the  flow  Just  becomes  choked.  Cnee  the  throat 
becor.es  choked,  the  solutions  for  both  nozzles  are  Identical.  The  region  to 
the  left  of  the  dotted  line  Is  the  region  where  the  flow  Is  not  choked.  The 
numerical  results  Indicate  that  there  Is  a  «*  for  which  5  vanishes,  and  an 
explicit  relation  la  obtainable.  The  solution  Is  presented  in  Figure  4.  The 
solution  for  T-  is  broken  down  Into  three  ranges  of  %  >  for  1  <  %  $  , 

vhich  Is  the  region  where  the  flow  is  not  choked|  for?  $  $  £,  which  Is  the 
region  where  the  flow  In  the  throat  1b  choked  hut  the  flow  behind  the  Inci¬ 
dent  shock  wave  la  still  subsonic;  end  for  ?  £  £  ,  which  Is  for  the  case 
where  the  flow  behind  the  Incident  shock  wave  Is  supersonic.  These  equations 
are  presented  In  Figure  5,  and  you  will  note  the  small  difference  between  the 
two  types  of  nozzles.  The  solution  for  the  latter  range  (|  (  (  }  i«  the  tri¬ 
vial  one,  namely  «  1. 

The  strength  of  the  reflected  wave  (  for  4  open  end  as  no  signal 

ecme*  back,  either  shock  ways  or  rarefaction  wave)  were  also  obtained  and  are 
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presented  in  Figures  6  and  7.  The  difference  between  the  two  types  of  nozzles 
is  small,  and  therefore  only  these  two  plots  are  presented  here. 

Since  we  were  primarily  interested  in  the  interaction  of  a  peaked  blast 
wave,  (which,  Incidentally,  is  generated  by  detonating  a  "plane"  charge  of  high 
explosive)  a  solution  was  obtained  fora  typical  blast  wave  as  encountered  in 
the  Shock  Tube  Laboratory.  For  these  calculations  a  shock  strength  of  two  was 
selected  so  that  the  entire  flow  field  could  be  considered  isentropic.  The 
solutions  were  obtained  by  the  method  of  characteristics;  a  graphical  treatment 
was  employed  in  this  particular  case.  As  a  basis  for  comparison,  solutions 
were  obtained  for  the  given  blast  wave  incident  upon  (l)  a  closed  end,  (2)  an 
open  end,  and  (3)  a  partially  open  end.  The  area  ratio  used  in  the  latter  case 
was  one  which  completely  eliminates  the  reflected  signal  from  the  incident 
shock  wave  (  «*  =  0.6).  The  flow  in  the  throat  of  the  nozzle  was  initially 
choked.  The  wave  diagrams  for  these  three  cases  are  presented  in  Figure  8. 

Time  is  measured  in  units  of  the  positive  phase  duration  (t0)  of  the  blast 
wave  (at  the  initial  position),  and  distance  is  measured  in  units  of  aQt0  where 
a0  is  the  ambient  sound  velocity.  The  wave  attenuates  as  It  propagates  down¬ 
stream  toward  the  end  of  the  shock  tube.  In  Figure  8  you  see  the  reflected 
shock  wave  from  the  closed  end,  the  reflected  rarefaction  wave  followed  by  a 
compression  wave  for  the  open  end,  and  finally  the  much  less  disturbed  flow 
field  for  the  critical  area  ratio  grid.  The  flow  is  initially  choked  in  the 
nozzle--the  dashed  line  represents  rarefaction  waves.  Then  the  boundary  con¬ 
dition  changes  as  the  pressure  drops  so  that  the  pressure  in  the  threat  is 
equal  to  ambient  pressure  (the  pressure  Just  outBide  the  tube);  the  solid  lines 
represent  compression  waves. 

The  effecti  .eness  of  the  grid  can  be  seen  in  Figure  9,  whore  the  pres¬ 
sure  variation  is  presented  for  the  three  cases  at  three  positions  upstream 
from  the  grid.  First  we  have  the  position  x/a0to  ■  0,  which  is  the  position 
where  we  imposed  the  given  blast  wave,  then  at  x/a0t0  =  0.5,  and  finally  at 
x/a0t0  ■  0.75*  The  solid  line  represents  the  pressure  variation  if  the  tube 
were  infinitely  long  (  no  end  effect).  The  dot-dash  line  is  for  the  case  of 
a  closed  end.  The  dotted  line  is  for  the  case  of  an  open  end,  and  the  last 
line  is  for  the  case  where  the  grid  area  ratio  was  approximately  0.6.  The 
grid  lias  the  effect  of  decreasing  the  magnitude  of  the  disturbance  and  in  de¬ 
laying  its  arrival  at  a  given  station  back  upstreoa. 

I  will  now  discuss  briefly  the  results  of  the  case  where  a  grid  la 
placed  in  a  constant  area  channel  and  1b  struck  by  a  flat-topped  blast  wave. 

We  will  refer  back  to  Figure  2.  The  nozzle  used  in  treating  this  case  was 
the  converging-diverging  nozzle.  If  we  neglect  friction,  as  was  done  pre¬ 
viously,  we  would  expect  no  effect  of  the  grid  until  the  strength  of  the 
incident  shock  wave  was  great  enough  to  choke  the  flow  in  the  throat  of  the 
nozzle.  Thor,  both  a  reflected  and  transmitted  shock  wave  exist.  As  the 
strength  of  the  incident  shock  wave  is  increased  a  third  shock  wave,  a 
standing  shock  wave,  exists  in  the  nozzle.  As  the  shock  strength  of  the  in¬ 
cident  wave  4  further  increased,  this  auxiliary  shock  wave  is  swept  down¬ 
stream  behind  he  transmitted  shock  wave.  The  same  simple  one -dimensional 
equations  are  .:sed  as  before,  however,  due  to  the  shock  pattern  there  must 
exist  an  entropy  discontinuity.  It  is  shown  in  Figure  2  aa  the  "interface*. 
The  conditions  at  the  Interface  require  that  the  pressure  and  particle  velo¬ 
city  of  ste  tea  3  be  Identical. 
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Again  a  numerical  procedure  must  he  used  to  obtain  the- results.  The 
strength  of  the  reflected  shock  wave  (t,  '}  Is  the  same  as  shovn  before  in 
Figure  6.  The  strength  of  the  transmitted  shock  wave  (<£")  *8  presented  In 
Figure  10.  The  region  to  the  left  of  the  dotted  line  Is  the  region  where  a 
standing  shock  wave  exists  In  the  nozzle.  The  strength  of  the  auxiliary 
shock  wave  was  determined  but  is  not  presented  here.  The  strength  of  this 
shock  wave  was  quite  strong  and  had  Its  maximum  value  at  the  breakaway  point, 
that  Is,  along  the  dotted  line. 

Some  crude  estimates  were  made  as  to  the  effect  of  friction  and  these 
estimates  Indicated  that  friction  could  be  neglected.  The  area  ratio  la  the 
principal  variable  governing  the  interaction  of  the  blast  wave  vith  the  grid. 
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Figure  U.  Solution  for  Complete  Elimination  of  Reflected  Signal 
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Figure  5.  Equation  of  Solution  for  Complete  Elimination  of  Reflected  Signal 
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Figure  6.  Reflected  Shock  Strength 
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action  Wave 


Figure  8.  Wave  Diagrams  for  Peaked  Wave  Incident  upon  End  of  Shock  Tube 


Figure  9*  Pressure  Variations  near  the  End  of  the  Shock  Tube 

for  Peaked  Wave  Blast 


Figure  10.  Transaitted  Shock  Strength 


DEVELOPMENT  0?  THE  SHOCK  TUBE  FACILITY  FOR  AIRFOIL  STUDTKB 

J.  R.  Ruetenik 

Massachusetts  Institute  of  Technology 


The  M.I.T.-V.A.D.C.  8  by  24-inch  shock  tube  vas  built  principally  to 
study  the  transient  forces  on  airfoils  and  other  aerodynamic  objects  upon 
immersion  in  a  blast  wave,  although  structural  response  experiments  also  occu¬ 
pied  a  prominent  role  in  the  planning  of  this  facility. 


A  delay  until  complete  measurements  of  the  flow  field  were  made  vas 
felt  inadvisable,  so  shortly  after  completion  of  the  shock  tube,  a  4-inch 
chord  double-wedge  airfoil  vas  installed  in  the  tube  and  the  measurements  pre¬ 
viously  described  by  Dr.  Witmerl  vere  made.  In  order  to  test  the  reliability 
of  the  measurements,  the  data  after  the  flow  field  had  achieved  a  "steady 
state",  which  was  conservatively  chosen  as  20  chord  lengths,  vere  compared  to 
wind  tunnel  measurements  made  at  Cornell  Aero.  Lsb.^  and  Langley  Field’, 

Figure  1.  At  an  angle  of  attack  of  4  and  8.3  degrees,  the  lift  coefficient  in 
the  shock  tube,  which  Is  represented  by  circles  in  Figure  1,  was  respectively 
58  and  48)1  greater  than  in  the  wind  tunnels. 


A  number  of  explanations  came  to  mindt  the  possibility  of  reduced 
separation  at  the  comers  of  the  airfoil  by  the  heat  transfer  from  the  hot 
gas  to  the  cool  model,  or  by  a  larger  turbulence  level  due  to,  among  other 
causes,  diaphragm  particles  and  the  rupturing  process..  But  these  explanations 
were  considered  insufficient  since  the  values  of  the  lift  data  were  greater 
than  the  flat-plate  theory.  The  uncertainty  in  the  interferogram  fringe  mea¬ 
surement  and  the  Glad3tcne-Dale  constant  were  estimated  to  he  only  5  to  10)1. 
This  was  2-1/2  years  ago,  before  much  was  known  about  the  flow  following  the 
shock  in  a  shock  tube,  so  a  series  of  measurements  were  made  of  the  density 
history  without  an  airfoil  following  shock  passage,  similar  to  Mack's1*  mea¬ 
surements  at  Lehigh. 

The  interferometer  vas  located  at  the  test  section  and  adjusted  for 
white  light  fringes  oriented  perpendicular  to  the  shock  tube  axis.  The  cam¬ 
era  view  was  blocked,  except  for  a  narrow  slit  perpendicular  to  the  fringes. 
Numerous  strip  films  were  made  over  a  range  of  initial  pressure  ratios,  as 
well  as  slngle-3p \rk  exposures  at  various  delay  times.  Figure  2,  vhlch  vas 
from  a  run  with  M2  -  0.5,  shows  the  first  2-1/2  milliseconds  following  shock 
passage,  and  is  typical  of  the  runs.  In  Figure  2,  the  slit  image  is  vertical 
and  sweeps  frcm  left  to  right.  The  shock  passes  from  top  to  bottom  so  that 
the  upper  frii  .»s  which  are  behind  the  shock  appear  before  the  lower  fringes 
ahead  of  the  shock  terminate. 

To  test  for  vibrations  transmitted  through  the  interferometer  struc- 
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ture,  the  interferometer  was  raised  above  the  shock  tube,  the  tube  fired,  and 
a  strip-film  was  taken.  The  fringes  were  as  straight  during  che  45  millisec¬ 
onds  following  shock  passage,  which  is  the  longest  test  period,  as  they  are  in 
this  figure  preceding  shock  arrival. 

Consider  first  the  details  of  the  fringes  in  Figure  2  following  shock 
passage.  Note  the  steps  in  the  fringes  which  correspond  to  .weak  shocks  which 
have  not  yet  overtaken  the  initial  shock.  These  weak  shocks  are  often  seen 
in  the  single-exposure  interferograms.  In  addition,  there  is  a  considerable 
waviness  in  the  fringes,  which  indicates  that  the  flow  is  not  as  smooth  as 
shock  tube  flow  was  once  believed  to  be.  Of  course  the  interferometer  mea¬ 
sures  fluctuations  in  wall-boundary  layer  as  well  as  the  so-called  "free- 
streara",  but  calculations  indicate  that  the  fringe  shift  due  to  the  window¬ 
boundary  layer  is  quite  smalls  so  it  is  concluded,  for  the  present,  that  the 
fringes  Indicate  the  density  variation  in  the  freestream. 

These  fluctuations  originate  in  two  principal  sources.  First  is  the 
diaphragm  rupture  with  its  ncn  -planar  release  cf  the  high  pressure  gas,  and 
curved  Initial  shock  with  Mach  stems  and  reflected  waves,  etc.  But  correla¬ 
tion  of  attenuation  data  with  the  attenuation  theories  of  both  Hollyer*  and 
Mirels®  indicate  that  the  wall-boundary  layei  induces  pressure  waves  which 
are  more  important  to  attenuation,  and  could  quite  reasonably  be  causing  the 
density  fluctuations  in  the  freestream.  And,  of  course,  theso  waves  do  over¬ 
take  the  shock  since  the  shock  is  subsonic  with  reference  to  the  gas  behind. 
Further  evidence  of  the  character  of  the  freestream  flow  is  the  much  lower 
values  of  the  transition  Reynolds  number  of  the  wall-boundary  layer,  as  ob¬ 
served  at  this  laboratory  and  at  others,  than  would  be  expected  by  the 
increased. stability  of  the  laminar-boundary  layer  by  the  wall  cooling.  M.I.P. 
measurements?  at  M2  -  0.7  indicate  that  the  transition  Reynolds  number,  based 
on  the  particle  velocity  and  dis*ance  traveled  and  properties  behind  the 
shock,  is  about  2  million.  And  yet  Lee's  calculations”  for  a  flat-plate 
boundary  layer  at  this  Mach  number  and  wall  to  freestream  temperature  ratio 
appropriate  for  the  test  condition  cited  give  a  minimum  critical  Reynolds 
number,  based  on  distance  from  the  leading  edge,  cf  60  million.  Of  course, 
the  velocity  and  temperature  profile  are  come-.nat  different,  but  this  effect 
is  not  believed  to  be  significant.  Incidentally,  In  other  shock  tubes,  the 
transition  Reynolds  number  was  somewhat  lower,  which  may  be  due  to  the 
smaller  tube  diameter  and  the  consequent  increase  in  boundary  effects.  The 
M.I.T.  wall-boundary  layer  measurements  were  reported  at  the  annual  meeting 
of  the  American  Physical  Society?,  and  are  *0  be  submitted  as  a  W.A.D.C. 
report.  The  present  conclusion  would  be  that  the  freestream  needs  further 
investigation  in  order  to  compare  shock  tube  flow  to  wind-tunnel  flow.  For 
blast  wave  studies,  however,  the  flow  may  be  quite  similar. 

Return  now  to  the  strip-film  tests.  It  was  desired  to  obtain  the  time 
development  of  the  dynamic  pressure  during  +he  45  ’millisecond  test  period  from 
the  strip-films  in  order  to  explain  the  large  steady  state  lift  values.  The 
dynamic  pressure  was  computed  at  the  end  of  the  time  period  from  the  density 
using  three  models.  It  was  postulated  that  the  density  Increase  was  due  to 
(l)  a  stronger  shock,  (2)  a  second  shock  following  the  initial  shock,  and  (3) 
an  infinity  of  shocks,  an  isentropic  compression,  following  the  Initial  shock. 
The  values  of  the  dynamic  pressure  agree  within  one  percent  for  all  three  flow 
models. 
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The  density  history  compiled  from  ihm-/  runs  Is  shown  in  Figure  3. 

Here  s/c  is  the  chord  lengths  of  flow  at  the  ‘est  station,  and  c  =  1*  inches 
corresponding  to  the  two  airfoils  tested.  An  s/c  of  60  corresponds  to  40 
milliseconds  of  flew  for  this  case.  The  quantity  q  is  *he  dynamic  pressure 
following  Ehock  passage,  and  q2  is  the  value  immediately  behind  the  shock 
at  the  test  section  station.  The  solid  cur*e  is  an  average  of  the  data;  the 
shaded  region  indicates  the  data  scatter.  It  is  interestirg  to  note  that  the 
subsequent  boundary  layer  measurements  showed  that  the  boundary  layer  was 
laminar  until  about  s/c  =■  6.  The  greater  scatter  of  data  thereafter  apparent¬ 
ly  represents  the  effect,  of  the  irregular  growth  of  the  turbulent  boundary 
layer  during  each  run,  and  from  run  to  run. 

Braun  and  Mirels  have  computed  the  change  in  properties  at  a  fixed 
station  following  shock,  passage,  an  extension  of  Mirels  NACA  TN  3278  attenua¬ 
tion  work*.  They  worked  out  two  cases:  one  a  completely  laminar-boundary 
layer  from  the  shock  to  the  rarefaction  wave,  and  the  other  completely  tur¬ 
bulent.  We  computed  the  dynamic  pressure  rise  from  their  equations  for 
Mg  -  0.4  at  the  test  station  and  with  their  kind  permission  have  compared  the 
result  in  Figure  3*  The  agreement  with  the  theory  based  upon  a  turbulent 
boundary  layer  is  very  encouraging,  and  indicates  the  importance  of  the  wall¬ 
boundary  layer  on  the  flow.  The  value  with  a  laminar  boundary  layer  (at 
M2  ■  0.35)  io  much  nearer  unity. 

The  solid  curve  has  been  employed  to  correct  the  airfoil  measurements 
in  the  shock  tube,  as  represented  by  the  squares  in  Figure  1.  The  steady 
state  Cl  values  at  4  ani  8.3  degree  were  brought  to  within  30  and  24^ 
respectively  of  the  wind  tunnel  data.  It  is  consoling  however  to  note  that 
the  agreement  with  flat -plate  theory  is  excellent^  More  recently  similar 
measurements  have  been  conducted  in  the  M.I.T.  facility  on  a  more  conventional 
subsonic  airfoil,  the  NACA  6?. 010  profile.  At  an  angle  cf  attack  of  5  degrees, 
the  corrected  steady  state  data  are  only  pjl  above  wind  tunnel  data,  which  is 
quite  good. 

But  near  stall,  the  interferometer  that  has  been  employed  for  these 
measurements  is  at  a  disadvantage  for  several  reasons.  First,  it  measures 
all  along  the  span  giving  a  spanwise  average.  The  boundary  layer  on  the  air¬ 
foil  tends  to  separate  first  at  the  wlxdcw  ends  where  it  interacts  with  the 
wall-boundary  layer.  The  wind  t.urnel  measurements,  on  the  other  hand,  are 
generally  made  near  the  center  of  the  airfoil.  Secondly,  near  stall  a  vortex 
forma  on  the  top  surface,  and  there  is  some  uncertainty  as  to  how  to  determine 
the  surface  pressure  from  the  vortex  fringes  in  this  region  cf  strong  viscous 
effects  and  large  density  gradients.  Therefore  a  progiam  to  develop  a  pres¬ 
sure  Bensor  to  measure  the  mid-span  pressure  directly  is  proceeding  currently 
in  this  laboratory.  Incidentally,  the  difference  between  the  double-wedge 
pressure  profiles  and  the  wind  tunnel  profiles  was  almost  entirely  limited  to 
the  front  upper  surface.  The  direct  pressure  measurement  is  expected  to  re- 


*  Described  in  a  private  communication,  and  to  be  forthcoming  in  an  NACA 

technical  note  entitled  "Nonuniformities  in  Shock  Tube  Flow  Due  to  Unsteady 
Boundary-Layer  Action",  by  W.  H.  Braun  and  H.  Mirels. 


move  this  uncertainty.  The  possibility  of  temporarily  obtaining  very  high 
lift  values  before  the  onset  of  stall,  which  has  been  indicated  by  interfero- 
grams,  has  intensified  the  interest  in  this  measurement. 

The  pressure  sensor  being  developed  here  consists  of  two  barium 
titanate  trystals,  one  6n  the  upper  surface  cf  the  airfoil  and  cue  on  the 
lower.  The  signals  are  fed  into  a  high  Impedance  differential  amplifier  so 
that  only  a  difference  signal  1b  recorded.  It  is  important  to  obtain  the 
difference  as  early  as  possible  in  the  network  since  at  M2  =»  O.t  the  dyna¬ 
mic  pressure  is  only  25jt  of  the  shock  overpressure,  so  to  measure  the  lift  to 
10JG  accuracy,  and  taking  Cl  ■  0(1),  would  require  measuring  the  difference 
to  2j  of  the  pressure  rise.  At  higher  Mach  numbers  this  requirement  is  less 
severe.  Crystals  were  selected  primarily  to  miniaturize  the  pickups.  The 
electronic  portion  of  the  pressure  measuring  is  now  satisfactory.  The  input 
Impedance  of  the  pickup,  lead,  and  amplifier  input  is  about  2  by  10  inch  ohms, 
and  the  capacitance  is  about  300  micro-micro  frrads,  with  a  time  constant  of 
one  minute— a  5  -second  time  constant  gives  only  a  one  percent  fall  off  in  a 
50- millisecond  test  period.  These  values  may  be  Juggled  later  to  reduce  the 
noise  from  the  large  grid  resistance  to  less  than  one  percent.  Efforts  are 
now  being  made  to  remove  the  heat  transfer  effects  on  the  gage  from  the  hot 
gas,  and  the  program  will  then  be  aimed  at  determining  the  accuracy  of  the 
sensor. 


Another  program  will  be  conducted  to  measure  the  mass  flow  and  temper¬ 
ature  behind  the  shock  with  the  hot  wire  anemometer.  The  equipment  is  ready 
for  testing  in  the  shock  tube,  and  will  follow  *he  program  of  the  pressure 
sensor  development. 

A  total-pre6sure  probe  with  a  fairly  short  inlet  passage  to  a  strain 
gage  pickup  has  also  been  constructed.  It  is  planned  to  measure  the  desired 
dynamic  pressure  more  directly. 

Finally,  plans  have  been  made  to  test  a  "moving"  airfoil  subjected  to 
a  blast  loading.  The  plan  is  to  use  the  initial  shock  to  establish  a  steady 
state  flow,  and  entrap  a  portion  of  the  shock  in  an  auxiliary  channel.  The 
channel  would  be  constructed  in  such  a  manner  V>at  the  reflecting  waves  would 
then  form  into  a  "blast"  wave  and  pass  over  the  airfoil.  Tests  involving 
varying  "blast"  wave  strengths  and  orientations  with  respect  to  the  airfoil 
chord  lino  will  be  included. 
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LIFT  COEFFICIENT 
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Figure  1.  Comparison  of  Steady-State  Life  Coefficient 
Versus  Angle  of  Attack,  Mg  3  0.4 
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Figure  2. 


White  Light  Strip  Interferogram  of  Flow  Field 
Arrival  with  M  ■  0.5 •  Interval  between  Marks 

One  Millisecond  - 


at  Shock 
at  Top  is 
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- Average  of  M.l.T.-W.A.D.C.  8X24  inch 

shock*tube  data.  *  0.4 

Broun  ond  Mirels  Theory  (unpublished) 

A  Boundary  layer  completely  turbulent  Mf*0.4 
•  Boundary  layer  completely  lominar  *0.35 


Figure  3.  Freestream  Dynamic  Pressure  Following  Shock 
Passage.  Computed  from  Interferometric 
Measurements 


TBE  SURFACE  FIIM  THERMOMETER:  A  VERSATILE  SHOCK  TUBE  TECHHlflUg 


Daniel  Bershader 
Lockheed  Aircraft  Corporation 


Thermal  measurements  In  the  shock  tube  must,  of  necessity,  he  performed 
with  rapid  response  instrumentation,  preferably  In  the  microsecond  range.  In 
recent  years  groups  at  Princeton  and  Lehigh,  as  well  as  several  other  labora¬ 
tories,  have  employed  the  thin  film  resistance  thermometer  In  order  to  measure 
surface  terperature  and  heat  transfer.  Such  films  are  typically  several  hun¬ 
dred  atoms  thick  and  display  a  resistance  of  the  order  of  100  ohms.  The  tech¬ 
nology  of  their  preparation  is  now  well  developed. 


.rfuen  current  xs  supplied  from  a  source  uhose  impedance  Is  high  compared 
to  that  of  the  film,  a  change  in  film  temperature  produces  a  proportional  change 
in  potential  difference  across  its  terminals.  Further,  if  Joule  heating  of  the 
film  can  be  neglected  by  comparison  with  the  heat  transfer  taking  place,  then 
the  latter  may  be  determined  in  a  straightforward  manner  from  the  time  variation 
of  film  temperature. 


The  present  discussion  centers  on  some  physical  aspects  of  the  trans¬ 
ient  thermal  behavior  of  film  thermometers  and  Includes  an  application  to  the 
study  of  shock-induced  laminar-boundary  layer  heat  transfer. 


Application  of  a  temperature  step  to  a  system  is  a  standard  way  to 
study  Its  transient  response.  We  begin  with  a  semi-infinite  slab  whose  sur¬ 
face  lies  In  the  plane  y  =  const.  If  the  temperature  i3  abruptly  changed  from 
T0  to  Ta  at  time  t  *  0,  then  the  surface  heat  transfer  varies  with 

-1/2  1 
t  '  according  to  the  relation 

3T\  ,  (Ta  *  V 

*y'8urf  ’  /r?r- 

One  can  as  well  formulate  the  reverse  problen^  which 
-1/2 

at  '  type  heat  transfer  as  the  forcing  function, 
the  surface  temperature  adjusts  stepwise. 

Consider  next  a  finite  slab  of  thickness  6  lying  on  a  perfect  in¬ 
sulator.  The  exact  solution  to  a  temperature  s'tep  applied  at  the  free  sur¬ 
face  of  ths  slab  is  given  on  page  84  of  Reference  1,  but  a  good  approximation 
is  obtained  by  equating  the  average  rate  of  heat  transfer  into  the  film  with 
its  thermal  energy  increase  rate.  To  do  this  we  denote  by  T(t)  the  average 


(1) 

consists  of  specifying 
It  then  turns  out  that 


temperature  of  the  slab  at  time,  t.  Then 


k 


Ta  -  T(t) 
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For  the  Initial  condition  T(o)  =  T0,  the  solution  is 


(2) 


T  -  T(t)  »  (T  -  T 
&  a  o 

where  T  Is  a  time  constant  given  by 


t  a 


X 

2 


(3) 


If  the  thermal  resistance,  R,  be  so  defined.  For  a  gold  film  2  micro-inches 
thick,  t  is  about  10"^  seconds.  In  other  words,  in  10“^  seconds,  the 
slab's  (film  would  be  a  more  appropriate  word  here)  average  temperature  has 
l/e  of  the  total  applied  step  yet  to  accomplish.  Through  an  obvious  rearrange¬ 
ment  of  the  exponent  in  Eq.  3,  the  state  of  affairs  Just  described  can  be  said 
to  apply  when  the  diffusion  dept  hv^a  equals  the  film  thickness  &  . 


The  time  constant  discussed  above  does  not  correctly  represent  the 
step-response  for  the  more  general  case  of  a  film  (or  finite  slab)  on  a  semi- 
infinite  backing  material  of  arbitrary  thermal  properties.  One  finds  then^- 
that  the  surface  heat  transfer  is  obtainable  from 
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A  few  special  cases  of  Eq.  4  are  pertinent.  Firstly,  if  A  -  0,  the  resul¬ 
tant  expression  applies  to  a  film  on  a  perfect  insulator  as  discussed  above. 
Secondly,  we  note  that  the  above  equation  corresponds  to  Eq.  1  for  times 
t  «  T  ,  This  is  reasonable  since  the  film  would  appear  as  a  semi-infinite 
slab  until  the  thermal  diffusion  penetrated  to  the  interface. 


Thirdly,  it  is  evident  that  for  very  long  times  the  exponentials  in 
Eq.  4  approach  unity.  The  series  may  then  be  summed  to  give 
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For  sufficiently  long  times  it  is  evident  (see  next  paragraph)  that  the  heat 
transfer  to  the  film  is  being  transmitted  to  the  backing  with  very  little 
further  increase  in  film  temperature.  Thus, 
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and  Eq.  4a  may  be  rewritten  as  follows 

iT\  a  £T\  B.  ^Ta  '  V  (kb) 
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Equation  4b  is  an  exact  duplicate  of  Eq.  1  for  the  step-response  of  a  semi¬ 
infinite  slab  made  of  the  backing  material.  This  is  the  desired  situation  in 
which  the  film  plays  a  negligible  role  in  the  heat  transfer,  but  behaves  as  a 
surface  thermometer  by  recording  the  true  surface  temperature  of  its  backing 
material.  It  is  important  to  rote,  however,  that  the  backing  material  is 
normally  some  insulator  such  as  glass,  Lucite  or  quartz,  with  thermal  proper¬ 
ties  considerably  different  than,  say,  the  m^tal  wall  of  a  shock  tube.  For 
this  reason  the  surface  heat  transfer  for  a  given  applied  temperature,  or  con¬ 
versely  and  more  important,  the  observed  temperature  change  for  a  given  applied 
heat  transfer  must  be  corrected  accordingly.  If  the  latter  is  of  the  t-1/2 
type  mentioned  earlier,  the  observed  modified  temperature  "step"  of  the  surface 
(see  subsequent  discussion  in  this  section)  is  inversely  proportional  to  the 
quantity  \/k/9cp.  A  film  deposited  on  a  cellulose  tape  backing  will,  for  ex¬ 
ample,  show  a  temperature  rise  approximately  50  times  that  taking  place  at  an 
adjacent  steel  surface.  Indeed,  such  "amplification"  proves  to  be  a  welcome 
factor  for  the  experimentalist. 

Residual  ambiguity  in  the  above  discussion  centers  in  the  terms  "very 
long  times"  and  "sufficiently  long  times".  To  resolve  these  matters,  we  in¬ 
quire  first  for  the  film-backing  interface  temperature  as  a  function  of  time 
following  application  of  a  temperature  step  to  the  film's  exposed  surface. 

From  Ipage  263,  we  have 
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It  turns  out  that  a  gold  film  with  J  =*  50QA  (  X  «  1.2  x  10  sec)  placed  on 
a  cellulose  tape  backing  displays  an  interface  temperature  rise  equal  to  8 9^ 
of  the  applied  step  T  -  T  in  a  time  2  t  ( =RC ) .  Thus,  the  average  tempera¬ 
ture  of  the  film  rises* to  considerably  over  ?Of>  of  the  applied  temperature 
step  in  this  negligible  period  of  time.  This  i6  generally  consistent  with 
the  result  of  the  more  simplified  calculation  resulting  in  Eq.  3»  Further 
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calculations,  however,  would  show  a  longer  asymptotic  tail  on  the  temperature 
rise  curve  for  the  case  of  the  arbitrary  tacking  material,  reflecting  the 
heat  transfer  across  the  interface. 

Again,  for  purposes  of  the  application  discussed  in  the  next  section, 
there  is  more  immediate  interest  in  the  surface  heat  transfer  as  a  forcing 
function  rather  than  the  temperature.  Neglecting  thermal  resistance  of  the 
film,  Bromberg^  has  obtained  a  solution  for  a  type  heat  transfer  into 

a  small  thermal  capacitance  which,  in  turn,  rests  on  a  backing  material  of 
arbitrary  thermal  properties.  The  energy  transfer  is  expressed  by 
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His  result  reveals  a  characteristic  time  parameter  for  the  asymptotic  rise 
in  question  given  by 
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The  time  response,  T(t),  of  the  film  is  of  the  "inertial"  type  with  a  steep 
rise  in  the  beginning,  followed  by  a  highly  tapered  curve.  It  is  found  that 
57^  of  the  asymptotic  rise  takes  place  in  a  time  t  =»t*,  02$  for  t  *  9“**# 
and  9I49E.  for  t  »  100  x*.  For  the  gold  -  "Scotch"  tape  combination  of  our 
previous  example,  x*  =•  0. 27^  sec  onus .  The  total  rise,  incidentally.  Is 
identical  with  that  represented  by  Eq.  1  in  which  the  backing  material  plays 
the  role  of  the  slab. 


By  comparing  the  above  result  with  the  earlier  discussion  of  the 
semi-infinite  slab,  we  see  the  effect  of  the  film  on  the  transient  response. 
The  reversibility  of  forcing  functions  and  effects  that  obtains  between  a 
surface  temperature  s*ep  and  a  t"l/2  surface  heat  transfer  for  the  semi- 
infinite  slab  no  longer  strictly  holds  for  the  film  and  backing  combination, 
as  both  Eq.  4  and  the  Brouiberg  calculation  show.  However,  outside  the  time 
resolution  represented  by,  say,  lOO*':*,  one  does  have  an  effective  etep- 
temperature  reaction  to  a  fl/2  surface  heat  input. 

3  4 

Recent  experiments  *  have  shown  that  the  temperature  rise  observed 
when  a  travelling  wave  passes  over  a  film  thermometer  located  on  the  side 
wall  of  a  shock  tube  is  essentially  stepwise  in  nature*.  This  responee,  as 
we  have  seen,  implies  a  surface  heat  input  varying  with  t"l/2.  Examination 
of  the  laminar-boundary  layer  formed  by  the  flow  following  the  wave  confirms 
this  implication.  The  layer  in  question  is  a  quasi-steady  type  since  its 
structure  is  time- independent  when  viewed  from  a  frame  moving  with  the 


Figure  1  is  an  oscillogram  demonstrating  the  rise 


shock  wave.  There  are  several  similarities  between  this  layer  and  the  steady 
state  boundary  layer  formed,  say,  on  a  flat  plate  in  a  wind  tunnel  flew.  In 
particular,  the  shock-generated  layer  displays  a  heat  transfer  rate  varying 
with  the  inverse  square  root  of  distance  from  the  shock  wave.  Since  the  shock 
itself  is  assumed  to  move  at  uniform  velocity,  distance  from  shock  is  propor¬ 
tional  to  time  to  an  observer  in  the  laboratory  reference  frame.  Thus 
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The  thin  film  thermometer  provides  two  bits  of  information  when  applied 
to  the  present  problem.  Firstly,  it  gives  a  direct  check  on  the  wall  tempera¬ 
ture  rise,  providing  the  wall  temperature  boundary  condition  required  for  solu¬ 
tion  of  the  energy  equation  and  subsequent  comparison  between  theory  and  ex¬ 
periment.  Secondly,  one  can  deduce  from  the  measured  temperature  step  the  mag¬ 
nitude  of  the  heat  transfer,  leading  to  a  numerical  determination  of  the  para¬ 
metric  constant  for  the  process. 

In  the  experiments  described  in  Reference  3>  the  shock  Mach  number  was 
about  2.5  and  the  observed  temperature  rise  was  2*C  for  the  gold -cellulose  tape 
film  assembly,  corresponding  to  around  O.OU^C  for  the  rise  at  the  steel  wall. 
Clearly,  then,  one  could  use  a  constant  temperature  wall  approximation  for 
solution  of  the  boundary  layer  equations.  Combination  of  the  boundary  layer 
heat  transfer  expression  with  Eq.  1  yields  the  following  expression^  for  the 
desired  parametric  heat  transfer  -skin  friction  constant  (please  refer  to 
list  of  symbols). 
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Satisfactory  agreement  was  found  among  values  of  determined  in  this 
fashion,  by  interferometric  evaluation  of  the  boundary  layer,  and  by  theory. 

By  neglecting  heat  conduction  in  the  solid  in  the  direction  of  travel 
of  the  shock  wave,  we  have  been  able  to  treat  the  present  problem  as  one¬ 
dimensional.  Justification  for  this  approach  has  been  given  by  Rott  and 
Hartunians^  with  respect  to  the  shock  wave  the  .two-dimensional  heat  equation 
reads 
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In  analogy  with  Prandtl'a  classical  treatment  of  the  incompressible,  isothermal- 
boundary  layer,  one  non -dimens  ionalizes  the  above  equation  such  that  deriva¬ 
tives  are  of  equivalent  order.  This  is  done  by  introducing  reference  values  x 
and  S  respectively  and  cet  ° 
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Equation  9  then  becomes 
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In  line  with  the  assumption  that  6<<  x  ,  certainly  Justified  for  values  of 
x  larger  than,  say,  10  shock  wave  thicknesses,  the  second  term  on  the  right 
may  be  discarded.  To  show  this  more  explicitly,  ve  note  that  by  adjusting 
the  derivatives  to  have  the  same  orders  of  magnitude,  one  implies  from  Eq.  10 
that  the  non-dimensional  diffusivity,^*,  is  of  unity  order.  This  leads  di¬ 
rectly  to  the  following  estimate  for  4  /x 


For  a  wave  travelling  at  Mach  2.9  over  a  steel  wall,  one  computes  a  value  of 
6/x  ■  0.0016  for  xQ  ■  1  inch.  The  only  difference  from  the  Prandtl  treat¬ 

ment  is  that  here  we°deal  with  the  thermal  diffusivity  instead  of  the  kine¬ 
matic  viscosity  so  that  the  quantity  u_xo  may  be  interpreted  as  a  type  of 
thermal  Reynolds  number.  Because  "~2~~  thermal  diffurivities  of  typical 

solids  are  smaller  than  kinematic  ^  viscosities  of  typical  fluids, 

the  thermal-boundary  layer  approximation  in  the  solid  is,  as  Rott  points  out, 
better  than  the  equivalent  viscous  boundary  layer  approximation  in  the  fluid 
for  this  type  of  problea. 
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external  stream  adjacent  to  boundary  layer 


indicates  non-dimensionalization,  according  to  discussion  in  text 

distance  parallel  to  wall  of  shock  tube,  measured  backward  from 
shock 

coordinate  perpendicular  to  vail,  measured  from  surface 
time 

velocity  (with  respect  to  shock  wave  unless  otherwise  noted) 
temperature 

temperature  rise  observed  on  film 

thermal  conductivity 

average 

(a)  thickness  of  film  resistance  thermometer 

(b)  thickness  of  thermal  layer  in  solid 
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iC  P 

time  constant  of  film  «  ~ 
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specific  heat  at  constant  pressure 

WA  T  _9 

Angstrom  unit  =  1*  x  10  7  inches 
heat  flux 

-l/2 

time  parameter  for  t  '  type  heat  transfer  U>  a  film  on  the 
surface  of  an  arbitrary  backing  material  ■  cf/  ). 

i #  p  0 

parametric  constant  for  boundary  layer  viscous  shear  and  heat 
transfer  .  ,/9 
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viscosity  coefficient 
wall  shearing  stress^ 


Prandtl  number 


enthalpy 
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be  assumed  constant 


C  T _ 

P  P 

non-dimensional  stagnation  value  of  enthalpy  with  respect  to 
laboratory  frame 
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SOME  OBSERVATIONS  OF  BIFURCATED  REFLECTED  SHOCK  WAVES 


Arthur  Cohen  and  Roger  Strehlov 
Ballistic  Research  Laboratories 


The  bifurcation,  or  splitting,  of  the  foot  of  a  normal  shock  has  been 
studied  in  wind  tunnels*-. 

This  paper  will  describe  the  effects  ,of  R^  and  Mg  on  the  bifurcation 
of  a  normal  reflected  shock  traveling  down  a  tube  of  uniform  cross  section. 

The  investigation  was  undertaken  to  obtain  information  that  would  help  inter¬ 
pret  strip-film  Schlieren  pictures  of  shock  waves  reflecting  from  the  end  wall 
of  a  shock  tube.  Velocities  measured  from  these  Schlieren  pictures  can  be 
used  to  furnish  kinetic  information  about  the  gas  behind  the  reflected  shock. 

For  a  given  incident  shock  velocity,  the  velocity  of  a  reflected  shock 
depends  on  the  enthalpy  of  the  gas  behind  it.  The  enthalpy  of  an  ideal, 
reacting  gas  is  a  function  only  of  temperature  and  a  reaction  variable.  The 
relation  between  the  incident  and  reflected  shock  velocities  and  the  reaction 
variable  can  be  determined  by  numerical  methods.  By  experimentally  deter¬ 
mining  incident  and  reflected  shock  velocities,  information  concerning  the 
reaction  can  be  obtained** 

Comparisons  of  theoretical  and  experimental  results  in  non-reactive 
systems  appear  in  Figure  1.  The  theoretical  velocities  were  calculated1*  using 
variable  heat  capacities  and  the  ideal  equation  of  state. 

Figure  2  contains  strip-film  Schlieren  and  multiple-spark  shadowgraph 
pictures  taken  in  various  gases,  under  similar  conditions.  All  incident  shock 
traces  appear  on  the  Schlieren  pictures  as  thin  straight  lines.  This  allows 
incident  velocities  to  be  measured  to  $  l£.  The  reflected  shock  trace  in 
argon  is  thin  and  straight,  but  its  appearance  changes  in  the  other  gases. 

In  nitrogen,  the  reflected  shock  line  thickens  and  small  density  gradients 
appear  around  the  shock.  In  carbon  dioxide,  density  gradients  appear  near 
the  reflected  shock.  The  bending  of  the  reflected  shock  trace  indicates 
accelerations.  In  propane,  large  density  gradients  appear  in  front  and  be¬ 
hind  the  rerlected  shock.  At  higher  velocities  this  behavior  is  accentuated, 
and  accurate  velocity  measurements  in  the  latter  three  gases  become  difficult. 
Below  each  strip-film  is  the  corresponding  multiple-spark  shadowgraph.  From 
these,  it  appears  that  the  pressure  gradients  are  caused  by  the  bifurcation 
of  the  normal  reflected  shock  and  that  the  extent  of  bifurcation  is  greater 
in  gases  with  larger  heat  capacities. 

Conditions  under  which  bifurcation  occurs  in  wind  tunnels  has  been 
investigated?'”.  Examination  of  the  boundary  layer  shock  wave  interaction 


revealed  that  bifurcated  shocks  occur  more  readily  with  laminar-boundary  layers 
than  with  turbulent-boundary  layers  and  that  separation  of  the  boundary  layer 
accompanied  the  blf'ircatlon  process.  At  present  a  complete  theoretical  treat¬ 
ment  of  the  boundary  layer-shock  wave  Interaction  In  shock  tubes  has  not  been 
made.  A  simplified  dimensional  analysis  of  this  complex  Interaction?  Indicates 
that  the  pressure  coefficient  necessary  for  boundary  layer  separation  and  bi¬ 
furcation  Is  proportional  to  the  local  skin  friction  coefficient,  l.e.  Q  oc  C-. 
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constant  for  laminar 


boundary  layers  and  C^,  x  R  '  *  constant  for  turbulent  boundary  layers. 


Therefore 
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At  above  room  temperatures,  w<  1  and  the  products  are 


Inversely  dependent  on  Mg.  At  constant  R^,  is  Inversely  dependent  on 
For  shock  tubes.  It  can  be  shown  that  R^,  based  on  the  distance  a  gas  particle 


travels  before  meeting  the  reflected  shock  may  be  written  ast 
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The  pressure  coefficient  can  be  expressed  ass 
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Graphs  of  - — j  and  Q  as  a  function  of  Mg  for  argon,  nitrogen  ai.e  arbon 
o 

dioxide  appear  In  Figure  3»  The  data  for  these  graphs  were  taken  from 


Reference  2. 


All  pictures  were  taken  with  a  6  inch  off-axis  Schlleren  system.  The 
focal  lengths  of  the  two  parabolic  mirrors  are  30  Inches.  The  data  was  re¬ 
corded  on  a  35  “a  drum  camera  which  could  be  rotated  at  a  speed  of  2  inches 
per  millisecond.  The  Schlleren  pictures  were  taken,  using  a  No.  30  flash  bulb 
as  the  light  source,  a  0.025-lnch  rod  perpendicular  to  the  flow  as  the 
Schlleren  stop,  and  a  0.025-inch  silt  across  the  center  of  the  window  parallel 
to  the  flow  to  define  the  shock.  The  multiple -spark  shadowgraphs  were  taken 
using  an  Edgerton,  Germeshausen,  Grier  Stroboscope  Unit.  The  spa  c  duration 
was  ~  1.5  microsecond.  The  interval  between  pictures  was  1C'  nlcroseconds. 
The  sketches  in  Figure  2  explain  the  Schlleren  and  shadowgraph  r  tures. 

Figure  4  is  a  sketch  of  the  shock  tube.  It  is  made  fr  .3-lnch 
ID  seamless  steel  pipe  curved  into  11-foot  diameter  sections.  compression 


section  is  1$  feet  long.  The  expansion  section  consists  of  36  feet  of  curved 
pipe  and  2  feet  of  transition  section  which  leads  from  the  pipe  to  a  6-foot 
section  of  3  *>y  l-l/2-inch  "S"  hand  wave  guide.  The  length  of  the  window 
section  is  1  foot  and  contains  two  6-inch  diameter  plate  glass  vlndows.  The 
location  of  the  hack  wall  could  he  varied  from  8  inches  inBide  the  window 
section  to  2  feet  past  the  window  section.  Ms  was  varied  hy  varying  the 
pressure  ratio  aorose  the  diaphragm.  Variations  in  R^  were  obtained  hy 
varying  the  location  of  the  hack  wall  and  hy  changing  P0.  Since  an  area 
reduction  of  10£  is  encountered  going  from  the  circular  to  the  rectangular 
sections,  pictures  were  taken  of  gas  originally  in  the  wave  guide  section. 
Pressure  probes  were  used  for  synchronization  purposes. 

The  effects  of  Ri  and  Hs  on  the  appearance  of  the  reflected  shock  in 
argon,  nitrogen  and  carbon  dioxide  are  shown  in  Figures  5>  6  and  7»  The 
dependence  of  a  separation  distance  on  Q,  for  shocks  generated  hy  supersonic 
flow  past  concave  corners,  has  been  reported”.  A  qualitative  comparison  of 
shock  tube  results  with  the  results  derived  from  the  dimensional  analysis? 
can  he  made  hy  assuming  that  the  separation  distance,  S  (Figure  2)  depends 
on  Q  -  Qc. 

At  constant  Ms,  Q  is  constant  and,  according  to  theory,  decreases 
with  increasing  R^.  Examination  of  the  shadowgraphs  in  Figures  5  and  6  do 
not  reveal  any  strong  dependence  of  S  on  R^.  Examination  of  Figure?  indicates 
that  In  carbon  dioxide,  contrary  to  theory,  S  decreases  with  increasing  R^. 

Figures  5*  8  and  7  indicate  that  at  constant  R^,  S  and  therefore 
Q  -  die  Increases  with  increasing  Ms.  This  implies  that  Qc  is  inversely  de¬ 
pendent  on  Ma,  since  ^  decreases  with  Increasing  Ms  (Figure  3)«  This 
implication  is  in  agreement  with  the  dimensional  analysis  theory. 

The  effect  of  Ms  and  the  lack  of  R^  effect  on  the  separation  process 
has  been  noted  In  wind  tunnel  experiments  ”*H.  The  Inverse  dependence  of 
S  on  R^  has  also  been  noted  5>8  and  attributed  to  a  transition  from  a  laminar 
to  a  turbulent  boundary  layer. 

Some  characteristics  of  the  boundary  layer  behind  the  incident  shock 
appear  in  Figures  8  and  9«  Multiple-flash  Schlieren  pictures  In  propane 
taken  with  the  stroboscope  unit,  using  a  knife  edge  parallel  to  the  flow  as 
the  Schlieren  6top,  appear  in  Figure  8.  Strong  horizontal  density  gradients 
appear  at  the  reflected  shock  due  to  the  boundary  layer-shock  wave  interaction 
The  presence  of  the  incident  shock  is  manifested  by  the  appearance  of  the 
boundary  layer  in  the  third  frame  of  the  picture  on  the  left.  Turbulent 
motion  occurs  at  the  floor  of  the  tube  behind  the  reflected  shock. 

Figure  9  contains  oscilloscope  traces  of  pulses  from  a  thin  gold  foil, 
mounted  flush  with  the  top  wall  in  the  window  section  of  the  tube.  The  foil 
acts  as  a  temperature  gage9> 10  and  signals  the  passage  of  Incident  and  re¬ 
flected  shocks.  The  scope  was  triggered  by  the  Incident  shock.  The  initial 
rise  In  voltage  is  due  to  the  boundary  layer  traveling  behind  the  incident 
shock  and  the  large  rise  in  voltage  is  due  to  the  passing  of  the  reflected 
shock.  The  rise  time  of  the  latter  pulse  iB  much  longer  in  propane  than  in 
argon,  which  suggests  that  the  thickness  of  the  reflected  shock  at  the 


boundary  layer  Is  greater  In  propane 


The  effects  of  and  Ms  on  the  bifurcation  of  a  normal  reflected 
shock  traveling  down  a  shock  tube  appeared  similar  to  those  found  in  wind 
tunnel  experiments.  The  extent  of  bifurcation  was  found  to  be  greater  in 
gases  with  larger  heat  capacities.  The  dependence  of  the  reflected  shock 
appearance  on  R.  and  M3  was  compared  to  results  predicted  by  a  simplified 
dimensional  analysis  of  the  shock  wave  -  boundary  layer  interaction.  The 
shock  tube  results  agreed  with  the  predicted  Ms  effect  but  disagreed  with 
the  predicted  R^  effect.  The  initial  results  of  a  gold  foil  technique 
Indicates  its  usefulness  for  investigating  boundary  layer  phenomena. 
However,  a  theoretical  analysis  of  heat  transfer  effects  in  boundary  layer 
flow  is  necessary  before  a  correlation  of  the  experimental  results  with 
boundary  layer  flow  can  be  made  with  this  technique. 
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SYMBOLS 

=  =»  Reynolds  number 

•«  gas  velocity,  cm/sec 

=  distance  gas  travels  before  meeting  reflected  shock,  cm 
■  viscosity  coefficient,  poise 
=•  density,  gm/cm^ 

=  temperature,  degrees  Kelvin 
=  local  skin-friction  coefficient 
P  -  Pg 

*  a  pressure  coefficient 

1/2 

=*  critical  pressure  coefficient 

»  Mach  number  of  shock  wave  relative  to  oi 

o 

=  sound  velocity,  cm/sec 

»  static  pressure,  cm  of  Hg 

»  distance  traveled  by  reflected  shock,  cm 

-  dimensionless  density  ratio 

-  dimensionless  pressure  ratio 


W  s  molecular  weight,  gm/mole 

•  =»  ratio  of  specific  heats  at  constant  pressure  and  volune 


SUBSCRIPTS 


o  initial  conditions 

s  3  shock  conditions 

r.s.  =  reflected  shock  conditions 
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Figure  Shock  Tube  Diagram 


Multiple  Spark  Shadowgraph  Pictures  in  Nitrogen 
Effects  of  Rx  at  constant  Mg 

Effects  of  M  at  constant  R. 

s  1 

(Arrows  indicate  the  foot  of  the  reflected  shock) 
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Figure  7*  Multiple -Spark  Shadowgraph  Pictures  in  Carbon  Dioxide 
(Arrows  indicate  the  foot  of  the  reflected  shock) 
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Figure  8*  Multiple-Spark  Schlleren  Pictures  of  the  Boundary  Layer  In  Propane 
Left  -  Boundary  Layer  behind  Incident  Shock 

Right  -  Boundary  Layer  behind  Reflected  Shock  Knife  Edge  Parallel 
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Figure  9*  Oscilloscope  Records  from  Surface  Film  Thermometer 
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NUMERICAL  SOLUTIONS  FOR  THE  KSFLBCTIOW 
OF  A  COMPRESSION  WAVE  FROM  A  RIO  ID  WALL 


Theodore  H.  Schiffnan 
Armour  Research  Foundation 


Recent  observations  have  shown  that  in  some  instances  a  relatively 
slow  pressure  rise  exists  where  a  shock  front  Is  expected,  i.e.,  the  shock 
wave  associated  with  blasts  does  not  form  until  Borne  distance  from  the  ex¬ 
plosion  center.* 

The  purpose  of  this  report.  Is  twofoldj  first,  the  determination  of 
the  reflected  pressure  for  such  a  compression  wave  hitting  a  rigii,  infinite 
plate,  and  secondly,-  the  length  of  time  it  takes  the  pressure  on  the  wall  to 
build  up  to  ltB  fully  reflected  value. 

For  an  11.6-psi  compression  wave  the  reflected  pressure  turns  out  to 
be  slightly  higher  than  that  for  an  equivalent  shock  wave  and  the  build-up 
time  is  roughly  10$  greater  than  the  10-milllsecond  rise  time. 


It  must  be  emphasized  that  the  plate  vas  assumed  Infinite  In  extent, 
and  If  relief  takes  place  on  the  wall  of  the  actual  structure,  the  average 
reflected  pressure  will  not  build  up  entirely  to  the  predicted  value. 


S.  Chandrasekhar  has  outlined  a  numerical  method  for  determining  the 
reflection  of  a  blast  wave  from  a  rigid  wall.1  In  this  approach  he  utilized 
the  method  of  characteristics  for  one -dimensional,  lsentropic,  non-steady 
state  gas  flow.  The  non-linear  differential  equations  governing  this  type 
of  flow  can  be  written  In  operational  form  as2 


*  (»  *  c) 

*  <■*  - c)  -£rj 


-  0 


(1) 


-  0 


*  The  work  reported  In  this  talk  was  tyonsored  by  the  Air  Materiel  Command.  • 
18.  Chandrasekhar,  "The  Normal  Reflection  of  a  Blast  Wave",  BRL  Report  No.  439, 
o  W 

See  Section  2  for  meaning  of  symbols  used  In  this  report. 


Hence,  the  so-called  Rlemann  Invariants 


2  dx 

y-y  c+u  =  5c+u  =  constant  along 


■  u  +  c 


and 


(2) 


2  dx 

1  •  ■  j  • £  c  -  u  =  5c”U  =  constant  along  — 


■  u  -  c 


If  the  Initial  conditions  R(x,t  )  and  S(x,t  )  are  given  at  t  ,  R(x, t)  and 
S(x,t)  can  he  determined  numerically.  Choosing  a  time  increment,  At,  and 
defining 


Axr  ■  (u  +  c)  At 
dXg  =  (u  -  c)  At 

it  follows  from  Eq.  2  that  approximately 
R(x,t)  a  R(x  +  Axr,  t  +  At) 
S(x,t)  =  S(x  +  AXg,  t  +  At) 


(3) 


(U) 


From  such  Iterations  the  Rlemann  Invariants,  R  and  S,  can  now  he  computed  for 
all  times,  t,  at  a  new  range  of  distance,  x. 


namely 


The  boundary  conditions  are  stated  at  the  vail,  x  =>  0,  for  values  of  t 
u(0,t)  =«  0 
R(0,t)  =i  S(0,t) 

From  Eq.  2  we  may  solve  for  u  and  c  as  a  function  of  R  and  S, 


(5) 


„  R  +  S  a  „ 

C  3  "■■■■'  a  C 

10  Pc 


,  irt 


u  =■ 


R  -  S 


5  <=, 


W  ■ '] 


(6) 


Utilizing  isentropic  pressure-density  relations,  P  =  K/I^,  the  side-on 
pressure  becomes 


t-1]  -'•[fc)T-1]-  -N7-;]; 


psi 
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Table  1 


i 


ITERATION  PROCEDURE  FOR  OBTAINING  THE  REFLECTION  OF  THE  GIVEN  COMPRESSION  WAVE 


«<~c) 

»(ft) 

*Vfl> 

K(  ft/»«e) 

S(ft/Mc) 

u(f t/**c) 

c(rt/**«) 

«P*0  », 

5 

r 

0 

... 

5650 

5650 

0 

U)0 

16.7 

0 

0 

7.) 

— 

... 

5650 

619) 

-278 

118) 

20.5 

5.8 

ir.2 

— 

... 

5650 

6628 

-688 

1228 

26.) 

11.6 

( 

0 

2.26 

•2.26 

5860 

5860 

0 

1170 

18.7 

6.0 

0.002 

2.26 

1.81 

-2.92 

5650 

6060  . 

-195 

1169 

18.7 

6.0 

U.8 

1.68 

-3.63 

5650 

6628 

.688 

1228 

26.) 

11.6 

0 

2.J6 

-2.J6 

6060 

6060 

0 

1210 

2).7 

9.0 

o.oo6 

6.2 

1.95 

-2.73 

5650 

C3C0 

-350 

1200 

22.7 

8.0 

10.% 

1.68 

-3-6J 

5650 

6620 

-668 

1228 

26.) 

11.6 

r 

0 

2.62 

•2.62 

6260 

626o 

0 

1250 

30.5 

15.8 

o.oo6  ‘ 

6.0 

1.6J 

-3.11 

5650 

6620 

.685 

1227 

26.) 

11.6 

7.0 

1.68 

-  3-*»  3 

5650 

6628 

-688 

1228 

26.) 

u.6 

0 

2.50 

-2.50 

6650 

66)0 

0 

1290 

)6.9 

22.2 

2.0 

2.00 

•2.80 

6)10 

6580 

-1)5 

1289 

36.7 

22.0 

o.oot 

1  55 

1.68 

-3.61 

6150 

6628 

1277 

36.7 

20.0 

6.0 

1.70 

-3.6) 

5080 

6628 

-233 

12  78 

)0.7 

16.0 

7,4 

1.58 

-3.68 

5650 

6628 

-683 

1228 

26.) 

11.6 

0 

2.58 

-2.58 

6610 

6610 

0 

1)22 

66.1 

29.6 

0.6 

2. 35 

-3.02 

6560 

6628 

-  35 

1)13 

66.0 

29.) 

2.0 

2.20 

•  ).10 

6690 

6628 

-  70 

1)12 

%1.6 

26.9 

0.010 

6.0 

1.90 

-3.6) 

63IO 

6628 

-160 

1296 

)8.) 

2). 6 

6.0 

1.70 

-3.6) 

6100 

6628 

-265 

12f) 

33.5 

18.8 

s.o 

1.68 

-3.6) 

5850 

6628 

-390 

1268 

)0.6 

15.7 

0.1 

1.68 

-3.6) 

5650 

6628 

•608 

1228 

26.) 

11.6 

0 

2.66 

-2.66 

6628 

6628 

0 

1)26 

65.0 

30.) 

2.66 

1.75 

-3.16 

6600 

66?8 

-  15 

1)26 

65.0 

30.) 

0.012 

w.o 

1.50 

-3-6) 

6500 

6620 

-  66 

1)20 

66.1 

29.6 

6.0 

1.68 

-3.6) 

6260 

6628 

-185 

1289 

36.8 

22.1 

B.O 

1.68 

-3.6) 

60)0 

6628 

-)00 

1266 

32.7 

18.0 

10.6 

1.68 

-3.6) 

5650 

6628 

-688 

1278 

26.) 

11.6 
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where  c  is  expressed  lit  ft/sec,  for  standard  atmospheric  conditions. 

Figure  1  shows  the  form  of  the  given  compression  wave,  both  In  the 
pressure-time  as  well  as  In  the  pressure-distance  plane.  This  wave  form.  Just 
touching  the  wall  (x  =  0)  at  t  =  0  forms  the  starting  point  for  the  Iteration 
process.  It  Is  believed  that  the  chosen  time  Interval 

At  .  0.002  -  (8) 

leads  to  pressure  and  velocity  profiles  of  sufficient  accuracy. 

From  Eqs.  2  to  6  an  Iteration  table  associated  with  Figure  2  for  the 
Riemann  invariant  R  and  Figure  3  for  the  Riemann  Invariant  S  has  been  pre¬ 
pared.  This  iteration  process  is  shown  in  Table  1.  The  resulting  pressure 
and  velocity  profiles  are  shown  in  Figures  4  and  5/  respectively.  The  pres¬ 
sure  on  the  wall  as  a  function  of  time  is  shown  In  Table  2  and  Figure  6. 


Table  2 

PRESSURE  AM)  SOUND  VELOCITY  ON  WALL  AS  A  FUNCTION  OF  TIME 


Time,  t  Sound  Velocity,  c 
(6ec)  (ft/sec) 

Absolute  Pressure,  P 
(psi) 

Overpressure,  p 
(psi) 

0 

1130 

14.7 

0 

0.002 

1170 

18.7 

4.0 

0.004 

1210 

23.7 

9.0 

0.006 

1250 

30.5 

15.8 

0.008 

1299  ..  - 

36.9 

22.2 

0.010 

1322 

44.1 

29.4 

0.012 

1326 

45.0 

30.3 

o,.oi2 

1326 

45.0 

30.3 

We  note  that  the  maximum  reflected  pressure,  p  ...  -  30.3/  is  almost 
equal  to  that  which  would  ariBe  from  a  normal  ref lectibn  "of  a  flat -top  shock 
wave  having  the  same  pressure  as  the  flat  portion  of  the  incoming  compression 
wave.  This  is  not  accidental.  It  can  be  shown  that  the  reflection  coeffi¬ 
cient  of  a  compression  wave 


1 


(9) 


i 


prefl 

p^(o) 


1/7 

-  1 

<  -  1 


-  1 


lies  only  slightly  above  that  for  a  shock  wave  of  the  same  value  of  ^  , 

/  which  could  be  obtained  from  the  compression  wave  by  letting  the  rise  time 

vanish 


profl  6  +  8  C 

pJ.(o)  "  6  +  if  ' 


shock  ratio  (10) 


as  given  in  Table  3 


and  Figure  7« 


Table  3 

REFLECTION  COEFFICIENTS  FOR  COMPRESSION  WAVES  AND  SHOCK  WAVES 


Shock  Ratio 
or 

Compression-Wave  Ratio 


Incidence 
Overpressure,  p^. 
(psi) 


Reflection  Coefficients, 

_ Prcfl/  p<r _ 

Compression  Shock 

Wave  Wave 


1.00 

0 

2.000 

2.000 

1.79 

11.6 

2.611 

2.608 

2.00 

14.7 

2.75^ 

2.750 

3.00 

29.4 

3.379 

3.333 

4.00 

44.1 

3.905 

3.800 

5.00 

CO 

• 

co 

irv 

4.351 

4.182 

Figures  4  and  5  summarize  the  pressure  and  velocity  profiles,  respec¬ 
tively,  near  the  wall;  while  Figure  6  summarizes  the  pressure-time  relation¬ 
ship  on  the  wall.  From  Figure  4  we  note  that  the  compression  wave  becomes 
much  steeper  upon  reflection  and  Irs  a  new  rise  distance  of  9  feet.  The  peak 


I 


» 


\ 
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moves  away  from  the  wall  at  a  constant  rate  of  2.64  ft/At,  while  the  tall 
moves  at  only  1.48  ft/ At;  hence  a  shock  would  form  roughly  at  nAt,  where 


n  ■  r:&  n  “.itSit  ®  8  <u) 

i.e.,  from  the  time  t  =  0,  where  the  wave  Just  touches  the  wall  It  becomes 
a  shock  In 


t 


^incident  +  *refl 


10  +  16  =  26  msec 


(12) 


Without  colliding  with  the  obstacle,  the  longest  time  for  the  compression 
wave  to  turn  into  a  shcck  wave  would  be 


t 


Xrlse  distance 


(u  +  c) 

-  c 

J  (1228  + 

Pr-  11.6 

P«-  “  0 

4^8)  -  (1130)  29  msec 

(13) 


Thus  a  compression  wave  bee  ones  a  shock  wave  at  roughly  the  same  time 
regardless  of  whether  or  not  It  collides  with  a  rigid  wall. 


LIST  OF  SYMBOLS' 

local  velocity  of  sound,  ft/sec 

c  in  undistrubed  region  «*  1130  ft/sec 

Integer 

side-on  pressure,  gauge  psi 
reflected  pressure  on  wall,  gauge  psi 
absolute  pressure  In  compression  wave,  psi 
absolute  pressure  In  flat-top  portion,  psi 
atmospheric  pressure  =*  14.7  psi 
Riemann  Invariant  -  5c  +  u,  ft/sec 
Riemann  invariant  ■  5c  -  u,  ft/sec 


t  3  time,  sec 

u  =  particle  velocity,  ft/sec,  positive  when  fluid  moves  in  positive 
x -direction 

x  =  distance  from  vail,  ft 

/  =  ratio  of  specific  heat  for  air  -  1.4 

At  3  iteration  increment,  sec 

AxR  =  (u  +  c)  At,  ft 

AXg  «  (u  -  c)  At,  ft 

P  =  density,  slugs/inJ 

for  compression  wave 

for  shock  wave  (shock  ratio) 
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Figure  1.  Fora  of  Coopresslon  Wave  Colliding  with  Rigid  Vail 
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Figure  2.  Rlemann  Invariant,  R,  versus  Distance  from  Wall,  x,  for  Various  Times 
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Figure  3*  Riemann  Invariant,  S,  versus  Distance  from  Wall,  x,  for  Various  Times 
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Figure  6.  Pressure-Time  Relationship  on  Wall 
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Figure  7«  Reflection  Coefficients  for  Cotapression  Wave  and  Shock  Wave 


THE  USE  OF  THE  SHOCK  rPE  IN  HYPERSONIC  RESEARCH 


F.  R.  Riddell  and  Mac  C.  Adams* 
AVCO  Research  Laboratory 


During  the  past  two  years  the  AVCO  Research  Laboratory  has  been  en¬ 
gaged  in  the  study  of  the  new  physical  phenomena  associated  with  flight  speeds 
^up  to  the  satellite  velocity.  The  shock  tube  has  been  used  to  simulate  the 
conditions  (e.g.  tempera ture  and  pressure)  of -these  very  high  speeds  and 
various  measurements  liave  been  made  to  understand  the  new  phenomena. 

The  areas  of  research  include: 

1.  Laminar  and  turbulent  heat  transfer  through  high  temperature 
boundary  layers  with  dissociation. 

2.  Heat  transfer  by  radiation  from  high  temperature  air. 

3.  The  determination  of  the  electrical  properties  of  ionized  air. 

4.  Measurements  of  the  recombination  rates  for  air  atoms. 

It  is  not  the  intent  of  this  paper  to  give  specific  results  of  any  of 
the  investigations.  Many  of  these  results  are  still  classified  and  could  not 
be  presented  in  any  event.  The  purpose  of  the  paper  is  therefore  to  illustrate 
some  of  the  shock  tube  techniques  which  have  been  employed  for  research  in  the 
above  areas  and  to  utilize  some  of  the  results  of  these  investigations  to 
define  flight  regimes  where  new  physical  effects  occur. 

Three  basic  experimental  te'hniques  have  been  used  to  study  the  high 
temperature  gas  dynamics  cf  air.  The  first  is  to  make  measurements  in  the  air 
which  is  set  in  motion  and  heated  to  high  temperature  by  the  incident  shock. 

An  example  of  an  experiment  which  utilizes  the  shock  tube  in  this  way  is  illus¬ 
trated  in  Figure  1.  In  this  experiment  the  electrical  conductivity  of  ionized 
air  is  measured. 

A  second  technique  is  to  use  the  shock  t-,be  like  a  wind  tunnel.  In 
this  case  a  model  is  placed  in  the  flow  field  generated  by  the  Incident  shock. 
By  adjusting  the  incident  shock  speed  and  the  initial  pressure  in  the  tube, 
the  stagnation  conditions  for  different  flight  velocities  and  altitudes  can  be 
simulated.  In  this  way  flight  velocities  as  high  as  the  satellite  velocity 
have  been  simulated.  Heat  transfer  measurements  can  be  made  at  any  position  on 
the  model.  A  schematic  of  this  technique  is  illustrated  in  Figure  2.  In  order 
to  obtain  larger  freestream  Mach  numbers  and  higher  simulated  altitudes  this 


*  Thi6  paper  compiles  partial  results  of  several  areas  of  research  by  the 
scientific  staff  of  the  AVCO  Research  Laboratory. 
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technique  can  be  modified  by  the  use  of  an  expansion  nozzle  attached  to  the 
shock  tube*. 

In  the  case  of  experiments  vnere  an  appreciable  quantity  of  high 
temperature  gas  at  rest  is  desired,  the  Incident  shack  is  reflected  from  the 
end  of  the  shock  tube*  3%e  region  of  the  gas  at  rest  grows  with  time  as  tbs 
shock  wave  moves  back  upstream  and  observations  can  be  made  in  this  region* 
This  technique  has.  been  used  An  the  study  of  radiation  from  air  and  is  illus¬ 
trated  in  Figure  3* 

As  a  result  of  research  in  areas  sentiomei  previously,  flight  regimes 
can  be  defined  where  new  physical  effects  bsccm*.  important*  Sane  of  these 
effects  are  illustrated  on  Fig-ore  4  which  presents,  on  an  a Ititude -versus - 
flight  Mach  neater  graph,  boundary  lines  fcr  various  flight  regimes*  This 
figure  is  intended  to  give  a  qualitative  picture  of  the  velocity  and  altitude 
ranges  where  various  known  physical  phenomena  occur*  Calculations  are  all 
based  on  the  properties  of  air  in  thermodynamic  equilibr  ium,  evaluated  at 
stagnation  conditions* 

The  ordinate  is  a  logarithmic  vcale  ot  the  aacieat  air  density  rela¬ 
tive  to  sea  level  derrsity.  This  ray  he  interpreted  in  tai-a-s  of  altitude  from 
the  following  table: 


Ratio  of  ante  lent  to  sea  level  density 


—  V 

10  ~ 


Altitude  in  feet 

60,000 
110,000 
160,000  • 
240,000  * 
000 

360,000 


It  should  be  rccted  that  the  altitude  scale  is  very  nearly  linear*  Ihe  ab¬ 
scissa  is  flight  Mach,  nvafrer  based  oa  room  temperat^ire  speed  of  sound  or, 
approximately,  flight  velocity  in  thousands  of  feet  per  second* 


The  limits  of  the  reglens  shown  were  determined  as  follows! 

Jjt  Dissociation:  yf>  c t  the  original  a.ir  molecules  are  dissociated 
under  stagiatlon  conditions*  A  similar  line  at  higher  velocity  shows 
where  4 Ojt  of  the  molecules  are  dissociated* 

Ionization:  The  degree  of  ionization,  i*e*  tit  nu nfijer  of  electrons 
per  air  molecule  reaches  a  value  of  1*5  x  10"®  under  stagnation  con¬ 
ditions*  The  electrical  0*C*  conductivity  of  air  is  then  equal  to 
that  of  sea  water*  This  is  also  about  100  times  the  degree  of  ioni¬ 
zation  in  the  8  layer* 

Thermal  "Farrier":  The  heat  flux  to  a  c me  foot  nose  radius  is  equal 
to  the  thermal  radiation  fro*  a  surface  at  2300*8* 
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Radiative  Heat  Transfer:  Heat  transfer  by  radiation  from  the  heated 
gases  is  an  appreciable  fractlcn  of  the  aerodynamic  heat  flux. 

Surface  Catalysis:  Heat  flux  can  be  reduced  It#  if  atom  recombination 
at  the  bcdy  surface  could  be  prevented. 

For  orientation,  the  present  and  immediate  future  flight  regimes  of 
airplanes  Is  shown.  In  addition.-  typical  ballistic  r^-entry  trajectories  are 
indicated.  These  trajectories  are  actually  those  followed  by  a  spherical  body 
of  5*foot  radius  with  a  density  equal  to  that  of  water.  Steady  flight  in  re¬ 
gions  below  the  thermal  barrier  is  complicated  by  severe  heating  problems  * 
Ballistic  missiles  are  able  to  penetrate  deeply  into  this  region  only  because 
the  heating  is  a  transient  problem. 

Figure  5  shows  for  comparison  the  flight  regimes  which  can  be  simu¬ 
lated  in  present  laboratory  equipment,  including  wind  tunnels,  ballistic 
ranges,  and  shock  tubes.  The  ordinate  and  abscissa  are  the  same  as  in  Figure  4. 
For  reference,  so..e  of  the  curves  in  Figure  1  have  been  repeated. 

It  should  be  mentioned  that  the  flight  velocities  of  ballistic  ranges 
are  being  increased,  and  it  is  believed  that  some  have  already  operated  at 
around  Mach  number  14.  . 
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Figure  2.  (a)  Schematic  Diagram  of  Ccrventional  Shock  Tube 
(b)  X-T  (Space-Time)  Diagram  Showing  the  Progress 
of  Shock  Wave,  Expansion  Wave  and  Interface 
(Contact  Disc-mtinuity)  Following  the 
Diaphragm  Bunt 


i 

1 


Figure  (a)  Sketch  of  Shock  Tube  C  s*rj  for  Radiation  Studies.  The  Ionization  Pick  Upa  Were  Used  to 
Determine  the  Incident  Shock  Velocx  -y  (Mach  Nunber) 

Figure  (b)  Spacetime  Diagram  Showing  Progress  of  the  Z  iock  and  Expansion  Waves  Following  the  Diaphragm 
Burst.  The  Mixed  Zone  Separates  the  Driver  Gas  from  the  Air.  The  Highly  luminous  Region  on  Which 
Observations  Were  Made  Is  Numbered  (5) 
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FLIGHT  REGIMES  ACCESSIBLE  TO 
PRESENT  LABORATORY  EQUIPMENT 


v.- 


MACH  NO 


PRESSURE  MUL  HUICATIOM  IN  RE  -  ENTRAP  T  CORNERS 

T.  Schiffman,  R.  J.  Heyman,  A.  Sherman,  D.  Weiner* 
Armour  Research  Foundation 


In  recent  years  a  great  deal  of  analytical  and  experimental  work  has 
been  done  on  the  reflection  and  interaction  of  shock  waves  and  on  the  diffrac¬ 
tion  of  shock  waves  around  structures.  Of  prime  interest  in  this  study  is 
the  pressure  multiplication  occurring  when  a  shock  front  reflects  oblique 
from  a  corner,  that  is  when  the  incident  shock  front  intersects  one  of  the 
sides  of  the  corner  at  an  angle  less  than  a  right  angle.** 

The  above  phenomena  can  occur  with  a  fair  amount  of  frequency  in  the 
diffraction  of  shock  waves  around  obstacles.  For  instance,  the  iront,  or  up¬ 
stream,  wall  of  a  block  in  the  regular  reflection  region  of  a  shock  wave  may 
feel  instantaneously  much  high'-  '"’essures  than  the  front  wall  of  an  obstacle 
struck  by  a  shock  wave  of  the  s>. strength  at  normal  incidence.  If  an  ob¬ 
stacle  is  shielded  by  another  that  is  struck  either  by  a  wave  reflecting  nor¬ 
mally  from  the  ground  or  by  a  Mach  wave,  the  corner  formed  by  the  front  wall 
of  the  second  obstacle  and  the  ground  will  undergo  this  type  of  pressure  mul¬ 
tiplication. 

The  problem  of  an  infinite  two-dimensional  comer  struck  obliquely  by 
a  plane  shock  is  amenable  to  analytical  treatment  which  is,  however,  very 
cumbersome.  Hence,  an  approx5.ua te  solution  of  a  much  simpler  form  has  been 
obtained.  If  the  two-dimensional  corner  is  finite,  no  exact  treatment  is 
possible;  however,  an  approximation  is  presented  which  agrees  with  experi¬ 
mental  results. 


The  peak  pressures  resulting  from  oblique  reflection  of  a  plane  shock 
on  an  infinite  plane  were  first  given  by  J.  V.  Neumann.  The  basic  equations 
are  usually  expressed  in  terms  of  flow  velocities  relative  to  the  flow  ahead 
or  behind  the  incident  and  reflected  shock  fronts.  These  equations  are  well 
known^,  but  difficult  to  apply.  Application  is  less  difficult  if  the  equa- 
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P1 

tions  are  expressed  in  terms  of  incident  shock  strength  %  «■  —  ,  incident 

P2  0 

angle  reflected  shock  .strength  ^ '  -  p~  ,  and  reflected  angled.'.  (See 
Figure  l)  The  term  PQ  is  the  absolute  pressure  in  the  undistrubed  region, 

P^  is  the  pressure  in  the  incident  region,  and  Pg  is  the  pressure  In  the  re¬ 
flected  region. 


For  incident  angles  or  to  oc.  extreme,  beyond  which  no  regular  re¬ 
flection  can  take  place,  an  observer  moving  with  the  point  of  contact  an  the 
wall  would  see  the  region  near  the  wall  split  into  the  zones  of  uniform  flow 
with  pressures  PQ,  P^  and  Pg.  Since  no  normal  velocity  component  can  exist 
on  the  wall,  oneob  tains 


■  l  Jf»  „  i 

— —  cos  c*  «*  -  cos  oC '  (l) 


The  velocity  component  parallel  with  the  wall  can  be  computed  from  the  con¬ 
ditions  ahead  of  the  reflected  shook,  R,  or  behind  the  incident  shock,  I,, 


l  +  6^f-5(£-l)s  ln^cr  _  l  +  6  -  ;(  4  '  -  1)  cos^of '  ^ 

^4(6  +  ?  )  sin  /T  +  6  4'  sinoc* 


This  set  of  two  simultaneous  transcendental  equations  expressing of' 
and  in  terms  of  given  %  and  of  can  be  solved  approximately,  for  oC 
less  than  ^0  degrees  by  postulating 


of '( ^  ,  a)  =  oc  +  A0f'(  ^ ,  of)  and 
4'(*,oc)  =  $(f,0) 

With  these  assumptions,  and  neglecting  second  order  terms,  one  obtains  from 
Eq.  1  and  2 


Ad.1 


Liz  1)  <5  cos^C  -  1)  tan  OC 
-M  -  5  + 


(3) 


-  ldO  - 


i 


and 


a*'=  jail*- 

^  (6  +4)  (3  *4<)*  ** 


(4) 


These  relatione  are  presented  in  Figure  2,  which  shews  that  up  to  39.2  degrees 
the  reflected  shock  strength  is  roughly  equal  to  the  head-on  value  (Of  -  0). 
Solutions  to  Eq.  1  and  2  in  the  form  of  0f'  and  ( %  K*  -  1) 

(  4  -  i) 

(reflection  coefficient)  as  a  function  of  ^  and  cx  are  given  in  the 

literature2* 3 , 

The  direction  and  magnitudes  of  absolute  flow  velocities  can  he  found 
from  Figure  1.  Consider  u,  the  particle  velocity  relative  to  fluid  ahead  of 
I  and  u  the  particle  velocity  in  the  undistrubed  region  as  known.  Then  one 
obtains  ror  the  flow  velocity  in  the  incidence  region,  u^ 


u.  -  (u  2  +  u2  *  2  u  u  sinof)^2 
1  o  o 


(5) 


which  is  deflected  from  the  component  of  u  ,  parallel  to  the  incident  shock 
front  by  the  angle  jS .  That  is 


0  a  tan 


■1  (  uQ  sin  of  +  u  \ 

V  uo  c°8  a  / 


(6) 


For  the  particle  velocity  in  the  reflected  region,  Ug 


u^  »  j^2  +  u'2  +  u^  u'  sin  (C<  +  0 ('  -/3).J 


1/2 


(7) 


which  is  of  course  parallel  to  the  vail  again.  The  term  u'  represents  the 
change  in  particle  velocity  when  passing  through  the  reflected  shock. 

For  relatively  weak  shocks  moving  into  an  'mdistrubed  region,  one  has 
approximately  a',  and  u  «  u'.  For  this  case,  where  u.  ■  u, 

0  -  ^/2f  Eq.  7  reduces  to 


10  c. 


"2 


2  u  sin  (X  - 


-  1 


JT  yl  *  6  4 


sin  of 


(7a) 


where  cq  is  the  velocity  of  sound  in  the  undistrubed  region. 
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Figure.  3  illustrates  infinite  re-entrant  comers.  Equations  1  and  2 
are  applicable  both  to  the  case  of  the  shock  moving  into  a  medium  at  rest, 
i.e.,  the  shocks  moving  into  the  re-entrant  comers,  as  well  as  to  shocks 
travelling  into  a  region  of  finite  flow  velocity,  i.e.,  shocks  emerging  from 
the  re-entrant  comers. 

For  example,  in  the  case  of  a  60  degree  re-entrant  comer  (Figure  3^>)#  » 

the  emerging  shock  is  still  represented  by  Equations  1  and  2,  but  the  meaning 
of  the  symbols  must  be  reinterpreted:  the  incident  shock  is  now  denoted  by 
§  '  and  the  incident  angle  is  <*'*  a  4’’  -  <*' ,  where  4>  is  the  angle  of  the 
corner.  The  quantities  %  '  and  tv'  ere  known  from  the  previous  solution. 

In  terms  of  the  reflected  shock  and  the  reflected  angle  oC* ' '  one  now 

can  write 

t  i  0  i  t  n 

, .  ,  ~  -  cos  a."  o  —  cos  <X"'  (is) 

^'(6  +  <*•)  \/\  *  6 


and 

1  +  6  g'  -  5(  -  1)  sin2cl"  a  1*6  -  5(4"  -1)  cos*  oc"  ' 

yj %r{6  +^*)  sin  <X"  /"IT 6  4' *  sin  <*"• 


In  the  case  of  the  U5  degree  re-entrant  comer  (Figure  3*)>  the 
shock  ij'  is  undergoing  a  normal  reflection  resulting  in  Then  Eqs.  1 

and  2  are  transformed  such  that  §  — 9  4'  ~ *  * ' ,  01  ■'*<*'  *,  and  <*.  - ’ ' . 

Hence,  in  a  60-degree  comer  the  maximum  pressure  in  the  re-entrant  comer 
corresponds  to  Prc  =  P3,  a  region  which  is  shocked  three  times,  while  in  the 
45-degree  comer  it  is  P4.  a  region  shocked  four  times.  In  the  90-degree  re¬ 
entrant  comer  (Figure  3C)  lb  is  simply  P?,  a  twice -shocked  region.  Evi¬ 
dently,  in  the  45-degree  comer  (Figure  3*7  the  upper  wall  is  equivalent  to 
the  symmetry  plane  for  the  shock  Incident  symmetrically  into  a  90-degree  cor¬ 
ner,  such  as  would  be  formed  by  the  front  wall  and  the  ground  in  the  regular 
reflection  region. 

The  peak  overpressure  in  the  re-entrant  comer,  Prc  -  P0,  can  be 
expressed  in  terms  of  the  Incident  overpressure,  P1  '  po,  by  meanB  of  the 
reflection  ratio 


t 
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Equation  9  represents  an  upperbound  for  1  .  A  lover  bound  Is  found 
by  substituting  ^  ■  1,  the  condition  for  a  sound  wave,  and  obtaining 

$  -  »  (11) 

Equation  9  can  probably  be  extended  to  the  limits  of  regular  reflection,  for 
example,  up  to  ^  •  1.1  for  <f>  =■  J  .  For  non- integral  values  of  m,  the 
upper  bound  for  the  reflection  6  coefficient  can  be  considered  as 

lying  in  the  range 


where  m+  and  m_  are  the  nearest  larger  and  smaller  integers  of  Eq.  10. 

For  the  45-degree  corner  (or  the  90"dcBree  with  the  f.hock  at  45-degree 
incidence)  the  reflection  coefficient  is  at  least  4  and  for  4>  becoming 
smaller  the  reflection  coefficient  becomes  very  large  so  long  as  <x  is  less 
than  ^extreme  (see  Reference  3»  page  60).  There  exist,  obviously,  incident 
angles  smaller  than  ^extreme1  Buch  ^  ^  M  larger  than  ^extreme*  In 

such  cases  the  above  analysis  doe3  not  apply  rigorously. 

For  moderately  strong  shocks  (^  <  2),  a  correction  term  can  be  found 
by  taking  advantage  of  the  flat  portion  of  Figure  2,  which  shows  that 


By  retaining  only  the  first  order  correction  terms,  one  obtains, 
instead  of  Eq.  9 
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where  m  is  again  defined  by  Eq.  10,  For  m  =»  2,  Eq.  9a  yields  the  familiar 
head-on  reflection  value 


6  + 

TT 


ill 


Thus,  for  a  shock  %  =  1,51  entering  a  60-degree  re-entrant  corner 
(Figure  3b)  the  reflection  coefficient  found  approximately  by  Eqs.  9  and  10 
is  !T  =  4,8,  by  the  improved  Eqs.  9a  '’’’'d  10#  3  4,3  while  the  rigorous 

Eq.  1,  2,  la,  2a,  and  8  give  J  =  4,2. 

The  only  experimental  data  at  present  applicable  to  this  problem 
consists  of  (1)  a  series  of  shadowgraphs?  which  show  plane  shocks  1,22  —  *j 
i.  1.79  impingent  on  45-degree  and  60-degree  re-entrant  comers,  and  (2)  in- 
terferograms  obtained  in  the  AHF  8-inch  shock  tube  for  ^  =  1.18  impingent  on 
60-degree  infinite  and  45-degree  finite  re-entrant  corners,  (see  Fig'ire  ")• 

Table  1  compares  computed  and  measured  shock  front  velocities  for  the 
shock  fronts  emerging  from  the  re-entrant  comers .  Table  2  compares  computed 
and  measured  re-entrant  comer  pressure  multiplication  obtained  br  irterfero- 
grams. 


Table  1 

VELOCITIES  OF  SHOCK  FRONTS  EMERGING  FROM  RE-ENTRANT  CORNERS 


Michigan  s 
Experiment5 
Fig.  No. 

Corner 

Angle. 

(deg) 

Inc  ider.t 
Angle,  a 
(deg) 

Inc Ident 
Shock  Strength, 

Emerging 

^  Shock  Frent 

Velocity  (fps) 

Measured 

Computed 

e- 

H 

VO 

H 

45 

45 

3.51 

1200 

1250 

25,  26 

45 

45 

1.51 

1200 

1250 

13,  14 

45 

45 

1.22 

1175 

1172 

37,  38 

60 

30 

1.79 

1260 

1290 

•  Physically  there  are  no  infinite  re-entrant  comers  and  hence  the 
constant  zone  behind  the  reflected  shock  is  bounded  by  a  rarefaction  wave 
originating  at  the  free  corner  edge,  and  by  the  wall  and  reflected  shock 
front.  Tills  rarefaction  wave  intersects  the  reflected  shock  front  at  S'" 
where 
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Table  2 


COMPARISON  OF  EXPERIMENT  AND  THEORY 
BY  FRINGE  SHIFT  MEASUREMENTS  IN  ARF  SHOCK  TUBE 


Incident 
Shock  Strength 

* 


Incident 
O'/  erp  res  sure 
(psi) 


Reflection 
Ratio,  J 


Re-entrant  Corner 
Overpressure 
(psi) 


Theory  Exnerlcent  Th  <_ ry  Experiment 


#0*  Infinite  1.18  2.6  +  .15  3.1*6  3.67  9.1  9.6  +  .5 

Re-entrant 
Comer 


45°  Finite  1.18  2.6  +  .15  ^.25  ^.04  11.1  10.6  +  .5 

Re-entrant 
Corner 


for  small  ''j/'  and  Of  Sr  ex* ,  For  T|r  *  0,  one  obtains  of  =*  oc  .  ^  OL 

sonic  extreme. 

In  this  section  the  corner  treated  is  as  illustrated  in  Figure  4  i.e. 
a  right-angle  corner  with  one  infinite  and  one  finite  leg.  Consider  a  shock  '* 
front  impingent  on  such  a  corner  with  an  incident  angle  of  45-degrees  to  the 
legs,  as  shown  in  Figure  4a. 

The  reflection  of  the  shock  front  from  the  corner  is  as  follows  t  The 
shock  front  enters  the  corner,  reflecting  in  a  regular  manner  from  the  finite 
leg  (Figures  4a  and  4b).  The  shock  front  reflected  from  the  finite  leg  is 
both  weaker  than  that  reflected  from  the  infinite  leg  and  curved  due  to  the 
rarefaction  initiating  at  the  corner  of  the  finite  leg.  The  shock  fronts  re¬ 
flected  from  the  finite  and  infinite  legs  then  interact  (Figure  4c)  resulting 
in  the  transmission  of  two  separate  shock  fronts  as  indicated  in  Figure  4d. 

The  region  between  the  two  transmitted  shocks  (region  P„  Figure  4d)  has  been 
swept  three  times  by  shock  fronts,  while  the  comer  region  (region  Pj,.  Figure  4d) 
has  been  swept  four  times  by  shock  fronts.  Finally,  the  entire  finite  leg  is 
immersed  in  a  region  that  has  been  four  times  (Figure  4e)  swept  by  shock  fronts. 

The  pressure  in  region  P4  (Figure  4e)  must  be  less  than  the  pressure  in 
an  infinite  comer  at  an  equivalent  time,  since  one  of  the  shock  fronts  enter¬ 
ing  into  the  interaction  was  necessarily  weaker  than  would  have  teen  the  case 
if  the  comer  were  infinite. 

It  is  not  possible  to  find,  rigorously,  the  strength  of  the  shock  ini¬ 
tially  reflected  from  the  finite  leg  of  the  comer  (^ '*  Figure  '-b) .  However, 
an  approximation  can  be  found  which  replaces  the  variable  shock  strength  X  ** 
and  non-constant  region  P2*  with  a  plane  shock  front  and  an  average  constant 
zone.  Hence,  the  regions  P3  and  P4  are  taken  as  approximately  constant  since 
the  shocks  * '  and  ^ ' '*  move  into  assumed  constant  regions  P?  and  P2*. 

Consider^  ..hock  with  strength  4  incident  on  the  finite  leg  with 
incident  angle  (g )  -  where  ot  is  nhe  angle  between  incident  shock  front 
and  infinite  leg.  The  strength  of  the  shock  reflected  from  this  leg  is 
known  for  the  limiting  cases  of  head-on  reflection,  (^)  -oc  ->o,  and  glancin* 
incidence,  (|)-d-*(|).  For  the  former  one  has  2 

<**  «  */2  -«-*>  (13) 

while  the  latter  gives 


4  '*  O  1,  */2  -  <*  -*ty2  (14) 
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Assume  a  linear  variation  with  «  between  these  limitirg  cases.  With 
the  approximation  of  the  previous  development,  £  Z  '  (i.e.,  the  shock  reflected 
from  the  infinite-  leg  is  approximately  equal  to  the  incident  shock  regardless 
of  angle  of  incidence)  and  n/2  -  01  «  the  angle  between  the  shock  V*  and  the 
finite  leg,  one  obtains  by  interpolating  betveea  Sq3.  13  and  14 


where  a.  <  c<extreme>  ensure  regular  reflection  ct  the  Infinite  leg. 

I 

The  computed  strength  of  the  weakened  shock,  Eq.  15  has  been  checked 
with  several  sets  of  experimental  data,  both  for  regular  and  for  Mach  reflec¬ 
tion  on  the  finite  leg®.  j 

The  head-on  collision  of  two  unequal  shocks  has  been  treated  in  the 
literature  by  solving  the  exact  velocity  eq-aticr.s  numerically  for  certain 
shock  strengths?.  In  the  present  paper  a  single  eq-attor  dealing  with  shock 
strengths  is  solved  only  approximately,  hut  thus  approximation  suffices  for 
all  practical  applications. 

The  shocks  and  (?**  (Figure  4)  move  into  pressure  region  with 
uniform  flow  normal  to  the  shock  motions.  Tnis  velocity  will  give  only  a 
tangential  component  to  all  shocks  and  will  lave  no  effect  c-n  the  interaction 
problem  before  or  after  collision.  Fr.-m  the  charge  In.  particle  velocity 
across  the  shock  fronts  before  collision 


jr,  . 

Ota..  a — "  -1-  . 

2  ST  /TTT?‘ 


and  after  collision 


0  ♦  .  *  * 


f7;  /i  +  eif  •* 


06) 


'"2  £  1  *-l 

-u,  .  +  U~  =* -  "* - •  t 

31  ♦  oj" 


5^*  .tVAxJL  ,17) 


+  m*  - -  - - 

"  nr  / TTT * 


« •  * 


where  u,  is  a  flow  component 
sion,  one^  obtains  by  adding  Eos 
an  equation  for  and 


in  the  region  between  the 
.  16  and  1”  ard  equating 

v  it. 


shocks  after  colli- 
thc-m  to  each  other. 


C1 


J^JL. 

/mrf 


+  lltjLj, 
/l  +  6{'* 


’  <"  - 1  1  .  -  - 

£’•*  -  l 

\/TTFT<'\  2 

yi  +~TF* 

06) 


where  the  c ' s 


refer  to  the  sound  velocities  in  regions  1,  2,  and  2*. 
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If  Since  the  pressure  between  the  shocks  after  collision  is  uniform,  one 
jObtal the  reciprocal  relationship 

II  v  V'  -  <•* 

I  (19) 

I  ?■-£.  <"•£.  r*-^i 

Eitfminating  the  sound  velocities,  one  can  combine  Eqs.  18  and  19  into  one 
ffjfuation  with  the  single  unknown 


-  1  £»»  -  1 


^'(6  +  ^*) 


«"*  - 1) 


/T+  6X*  / 1  +  6  <’*  I  x  + 


1  +  6 


For  =  <('*,  Eq.  20  reduces  to  the  head-on  collision  with  a  rigid  wall, 
which  is  Eq,  1  with  'X  *  ot'  =0,  i.e.}  head-on  collision  of  two  shocks  of 
equal  strength. 

As  zero-order  approximation  assume  that  the  two  shocks  pass  through 
each  other,  viz.,  and  '*  3s  .  The  first-order  approximation 

is  obtained  by  letting 

*"*  =  +  i  >  or  =  t,"  +  f  (21) 

and  finding  the  correction  term  from  Eq.  20  ’ 


where 


r  -  *.  ?(*'■  g'»)  +  .?(  *'♦,  s- 

* "  *  'H  v,  $•*, 


- 1)  /i  +  6  4'*  -/^'*(6+4'*)  (22) 


*  *•.  V)  - 
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The  correction  terra  &  is  equal  to  0  for  either  t*  or  equal  to  1, 

and  thus  for  a  shock  colliding  with  a  sound  wave.  For  two  equal  shocks 
colliding  the  explicit  solution  is  known: 

*”*  ”  4'  -  <23> 

which  with  Eq.  21  provides  an  upper  bound  on  £  .  Values  of  £  are  6hcwn  in 
Figure  5  for  shock  strengths  in  the  range 

1  <  4'  <  2 

1  <  $'*  £  2 

For  all  practical  purposes  is  negligible,  i.e.. 

Having  obtained  p,  by  Eqs.  22,  21,  19,  and  15,  P4  can  be  deduced  by 
realizing  that  the  particle  velocity  in  region  4  is  zero  and  that  the  compo¬ 
nent  particle  velocity  normal  to  the  emerging  shock  4,,r  is  roughly  the 
particle  velocity  behind  the  incident  shock.  The  sound  velocity  in  region  3 
corresponds,  however,  to  a  thrice-shocked  region.  Writing  the  expression  for 
the  change  in  particle  velocity  across  4  "*in  terms  of  4”  and  canceling 

the  factor  5c  7^7*  where  c  is  the  sound  velocity  of  the  undisturbed  re¬ 
gion,  one  obtains  a  quadratic  equation  for^'"i 


4-  1  . 

4(6  +  4) 

1/2 

4*(6  +  4') 

1/2 

4"(6  +  £•') 

1/2 

S'"  -  1 

v/l  +  6  £»' 

l  +  6< 

1  +  64' 

j 

1 

1  ♦  6<" 

/l  +  6$'" 

(24) 

whose  desired  solution  converges  toward  1  when  4  approaches  1.  The 
reflection  ratio  J  dor  the  finite  45-degree  re-entrant  corner  is  thus 

.  _<  ,25) 

d  - 1 

where 

<•  =■  f  4”  ■ 


from  Eqs.  21  and  15,  and  if1"  is  given  by  Eq.  24.  For  example,  a  shock, 
$  =  1.4l,  incident  into  a  finite  comer  will  lead,  according  to  Eq.  2’5,  to 
a  reflection  ratio  of  5.48,  while  the  solution  for  the  equivalent  infinite 
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corner  problem  gives  6.27 


For  shocks  Incident  under  an  angle  ot  different  from  1*5  degrees 
(Figure  6),  it  can  be  shown  that  the  emerging  shock  undergoes  a  shift  In 
angle  which  in  some  Instances  can  resemble  the  case  of  45-degree  incidence 
discussed  above  and  shown  in  Figure  4.*  The  peak  pressure  in  the  re-entrant 
corner  depends  again  on  the  approximate  head-on  collision  of  two  shocks  of 
unequal  strengths'.  However,  in  this  case  the  shocks  reflected  from  the  finite 
and  infinite  legs  move  now  into  a  region  where  there  exists  a  component  of 
steady  flow,  normal  to  these  shocks.  For  this  case,  Eq.  l6  would  be  modified  -to 


^  mw  • 


-U_  +  U~*  a  - 

^  a  {T 


(16a) 


-1 


Hence  by  adding  Eqs.  l6a  and  equating  them  to  Eq.  17,  one  obtains 
Eq.  18  again,  and  all  the  previous  results  are  applicable.  That  is,  according 
to  Eqs.  21  and  22,  the  two  shocks  cross  each  other  with  nearly  unchanged 
strength.  Also,  the  fact  that  the  shocks  do  not  collide  entirely  head-on  is 
not  too  important  since  according  to  Figure  2  the  reflection  ratio  is  nearly 
flat  and  equal  to  the  head-on  value  up  to  40  degrees.  It  is  believed  that 
according  to  Figure  5  the  reflection  ratio  for  equal  shocks  is  an  upperbound 
to  those  of  unequal  strength.  However,  the  velocity  vector  in  the  region 
between  the  shocks  is  deflected  towards  the  finite  leg  or  toward  the  weaker 
of  the  regions  2  or  2*,  Figure  6  as  can  be  shown  in  the  following  argument. 

For  the  case  of  the  infinite  comer  with  regular  reflection  on  each 
leg  (Figure  6a)  one  can  demonstrate  that  the  assumptions  5s  o<  and 

hence  c  «  c.  Cg,  lead  to  changes  in  the  particle  velocity  across  the  two 
shocks,  after  head-on  collision,  equal  in  magnitude  to  that  of  the  flow  behind 
the  incident  shock 


I  .!fo 


Hence,  in  region  3,  the  flow  vector  is  the  same  in  magnitude  and  direction  as 
in  region  1,  and  the  shock  emerges  parallel  to  the  incoming  shock. 


*  It  will  be  shovn  that  the  shift  also  occurs  when  £  =  45*  but  is  less  than 
the  shift  when  4  <45*. 
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In  ithe  case  of  the  finite  comer,  it  was  shown  that  the  two  unequal 
shocks  £| '  and  £'*,  pass  through  each  other  with  essentially  unchanged 
strength.  Consequently,  the  velocity  changes  across  these  shocks  also  remain 
unchanged,  (see  Figure  6h) 

The  quantity  |u.|  ,  the  velocity  change  across  ^ is  given  in  Eq.  26 
and  similarly  |u,#l  before  interac ting  with  is^given  as 


'C1  V*  -  1 

yr  j  i  *  64'* 


(27) 


Hence  after  interaction  of  shocks  and  the  velocity  vector  u^  asso¬ 

ciated  with  region  1  has  imposed  on  it  a  component  normal  to  the  shocks, 
and  -S'*,  which  is  given  by  the  difference  of  Eq.  2'6  and  27.  Hence,  the  flow 
vector  uj  is  deflected  from  by  j ,  where 


*,  t  .1^" 


.in.-.i.. 

/r+ 


1 


r1 


_sLn_2_0C; _ 

I  -  cos  2  C( 


sr*-*r 


(28) 


Thus,  Eq.  28  establishes  a  definite  tendency  for  the  emerging  shock  front  to 
shift  toward  the  weakened  region  2*.*  The  accuracy  of  Eq,  28  depends  upon  how 
precise  £j'*  is  giver.  >»y  Eq.  15  and  to  whan  degree  the  ass'xcption  of  parallel 
straight  shocks  holds  In  the:  neighborhood  of  the  comer. 


The  above  procedure  applies  approximately  ever,  if  there  Is  onset  of 
Mach  reflection  on  the  finite  leg,  since  the  strength  cf  the  weakened  shock 
Eq.  15  has  been  also  verified  for  seme  cf  these  case?.  The  manner  and  extent 
in  which  the  Mach  stem  grows  Is  treated  experimentally  by  D.  R.  White  in 
Reference  6.  While  it  is  obvious  tr.at  if  there  exists  regular  reflection  on 
both  legs,  the  contact  points  must  reach  the  comer  simultaneously,  a  Mach  stem 
travelling  down  the  finite  leg  will  reach  the  comer  before  the  regular  re¬ 
flection  contact  point  or.  the  Infinite  leg.  However-  the  ratio  of  time  inter¬ 
vals. 


At  infinite 
At  finite  leg 


(29) 


ordinarily  differs  little  from  unity.  The  value  cf  ^  represents  the  pressure 
ratio  across  the  Mach  stem  (as  a  f  unction  of  ^  and  ‘  ^/2  -  <*■)  feund  on  any 
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*  See  Figure  !*d  and  Jd 


standard  oblique  reflection  chart,  such  as  given  in  Reference  3« 

Thus  it  is  believed  that  even  for  finite  corners  with  regular  reflec¬ 
tion  on  both  legs,  or  for  a  weak  Mach  reflection  on  the  finite  leg  only,  Eq.  25 
in  conjunction  with  Eqs.  15,  21,  22,  and  24  holds  approximately  for  «■  alBO 
different  from  45  degrees  as  long  as  c*  £  *  extreme*  ^or  exanPle*  for 

%  =  1.62,  <X  a  20*,  the  reflection  ratio  J  for  the  corner  Figure  6b  is  5*37* 
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u  -  particle  flow  velocity  relative  to  fluid  ahead  of  I 
u*  -  Chance  in  particle  velocity  when  passing  through  R 
a)  Oblique  Reflection  Showing  Relative  Velocities 


■  Particle  flew  velocities 
b)  Oblique  Reflection  Showing  Absolute  Velocities 


Figure  1.  Oblique  Reflection  of  Plane  Shocks 
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Figure  5.  Heed-on  Colli j ion  of  T*o  Unequal  Shock* 
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a)  ShocK  Entering  60*  Infinite  b)  Shock  Emerging  from  60* 

Comer  (See  Figure  3b)  Infinite  Corner  (See  Figure  3b) 


c)  Shock  Entering  1*5*  Finite  d)  Shock  Emerging  from  1*5* 

Corner  (See  Figure  4b)  Finite  Corner  (See  Figure  Ud) 


Figure  7.  Interferograms  of  Shock  Entering  and  Emerging 

from  Re-entrant  Corners 
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